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We study smectic liquid crystals in random environments, e.g., aerogel. A low temperature 
analysis reveals that even arbitrarily weak quenched disorder (i.e., arbitrarily low aerogel density) 
destroys translational (smectic) order, in agreement with recent experimental results. A harmonic 
approximation to the elastic energy suggests that there may be no "smectic Bragg glass" phase in 
this system: even at zero temperature, it is riddled with dislocation loops induced by the quenched 
disorder. This result would imply the destruction of orientational (nematic) order as well, and that 
the thermodynamically sharp NA transition is destroyed by disorder. We show, however, that the 
anharmonic elastic terms neglected in the above approximate treatment are important (i.e., are 
"relevant" in the renormalization group sense), and may, indeed, stabilize the smectic Bragg glass 
and the sharp phase transition into it. However, they do not alter our conclusion that translational 
(smectic) order is always destroyed. In contrast, we expect that weak annealed disorder should have 
no qualitative effects on the smectic order. 
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I. INTRODUCTION 



A. Motivation and background 



The effect of quenched disorder on condensed matter 
systems is an important and challenging problem that 
continues to be actively investigated, because of its rele- 
vance to real systems, vi^hich always contain some amount 
of random inhomogeneity. Recently, much attention has 
focused on the random field XY model as a minimal 
model of a broad class oL systems such as disordered 
Josephson junction arraygj, roughening |Of crystal sur- 
faces growing on a disordered substratecl, the pinning 
of Wigner crystals, vortex lattices in superconductorsEl 
or charge density wavesu. Another system that is con- 
siderably less well understood is the superfluid transi- 
tion, e.g., He^'i-iin the random, "fractal-like" environ- 
ment of aerogelB. In apparent contradiction to the Harris 
criterionQ, the disorder modifies the critical properties of 
the transition (including the critical exponents). 

The central question addressed in all of the above stud- 
ies is: does a distinct low temperature "glass" phase exist 
in any of these disordered systems? And if so, what (if 
any) static property distinguishes it from the high tem- 
perature, conventional thermally disordered phase? The 
answer to the second question is clearly subtle: argu- 
ments dating back to LarkinQ show that it is impossible 
to have long-ranged translational order in three dimen- 
sions in a randomly pinned elastic medium. Hence, the 
glass cannot be distinguished from the high temperature 
"liquid" phase by long-ranged translational order. 



In the context of pinned vortex laltices in dirty type 
II superconductors, M. P. A. Fisher'scl original argument 
for the existence a vortex glass phase is based on a two 
(l-l-l)-dimensional model in which the random pinning 
is relevant at low tempyeratures but becomes irrelevant 
at higher temperatures.l3 The low temperature, pinning 
relevant phase is then identified as the glass phase, and 
the temperature at which the pinning becomes irrelevant 
is the thermodynamically sharp glass transition temper- 
ature. 

Unfortunately, this simple scenario cannot be directly 
carried over to three dimensions for vortex lattices 
(though, as we shall show later, it almost can for smec- 
tics in aerogel). In three-dimensional vortex lattices (and 
other pinned isotropic elastic media listed in the opening 
paragraph), the aforementioned Larkin argument shows 
that disorder is always relevant in d ~ 3. Hence, its rele- 
vance cannot be used as a criterion for distinguishing the 
glass and liquid phases. 

So wbaJ; can? One appealing proposal is the "Bragg 
glass"Lr1ij picture, in which the glass, while being elasti- 
cally disordered is topologically ordered, and is therefore 
distinguished from the liquid by being free of unbound 
dislocation loops, which proliferate in the liquid. The 
glass transition is then identified as an "unbinding" of 
the dislocation loops. 

This transition is then, qualitatively, very similar to 
the melting of the flux lattice in the absence of disor- 
der, which can also be thought of as an unbinding of 
dislocation loops. The only difference is that in the glass 
problem, the flux lattice is translationally disordered both 
above and below the transition. The absence of defects 
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below the transition, however, means that the low tem- 
perature solid phase still has a finite shear modulus, lead- 
ing to glassy behavior. 

Another related system that has been experimentally 
investigated is nematic liquid crystals in aerogel near and 
below the pure system's Nematic-Smectic-A (NA) transi- 
tion temperature.ta The mean-field theory of a bulk, pure 
smectic liquid crystal looks like that of an XY model, 
with a complex scalar order parameter ^ characteriz- 
ing the, sme ctic density wave. As was first noted by de 
GennesE^liJ a more precise model must include the cou- 
pling of ip to the nematic director fluctuations 6x1, jsihich 
takes the model out of the XY universality class .EtI13 Al- 
though the agreement between experiments and theory is 
far from perfect, the theoretical consensus is that the scal- 
ing near the NA transition should cross over from that of 
a three dimensional XY transition line, (as in a neutral 
superfiuid) to inverted XY-like behmdor with anisotropic 
X-ray correlation length exponents. E3 

The goal of this paper is to investigate the effects of dis- 
order on the liquid crystal phase diagram in the vicinity 
of and below the NA transition. Does the NA transition 
survive in the presence of even weak disorder? If so, what 
is the nature of the low-temperature phase? A condensed 
report of some of our results, the details and extensions 
of which are prese nted , in this paper, has appeared in 
recent publications .1111113. 

Our interest i n . thj s— problem was stimulated by re- 
cent experimentst3E3~t3 that aim to answer these very 
questions. In these experiments, a liquid crystal that ex- 
hibits a bulk Nematic-Smectic-A transition is introduced 
into an aerogel. X-ray scattering measurements show 
that this system never develops true smectic long-ranged 
(or even quasi-long-ranged) order: the translational cor- 
relation length S^-^ remains finite at all temperatures. 
This length smoothly increases as temperature is lowered 
across the bulk NA transition temperature, monotoni- 
cally and slowly rising to some finite asymptotic low tem- 
perature limit, which does not appear to be associated 
with any natural length scale of the aerogel itself. Our 
theory predicts this X-ray correlation length as the scale 
at which the pinning disorder energy begins to dominate 
over the smectic elastic energy, consistent with these ex- 
perimental observations; the smectic correlation length is 
not simply "pore" size of the aerogel.E3 This behavior is 
also in sharp contrast to the fast rise of the smectic corre- 
lation length to the nominal pore size at the bulk transi- 
tion temperature that one would expect in more regularly 
porous materials, which can be understood as a bulk-like 
sharp transition cut-off by the finite pore size. The spe- 
cific heat in aerogel was also measured and was found 
to exhibit a broadened but very well defined peak, near 
but slightly shifted down from the bulk transition tem- 
perature, Tna- These experimental observations sug- 
gest the destruction of the NA transition by the disorder 
imposed by the fractal aerogel pavironmcnt.Ej However, 
in contrast, recent experiments^, which study the NA 
transition in even lower density aerosils, appear to show 



a true resolution limited heat capacity singularity, while 
displaying a finite X-ray smectic correlation length in- 
dicative of short-range translational order. These latter 
experimental findings appear to support the idea that 
the Nematic-Smectic-A transition in pure systems, when 
confined in a low density quenched random structure, is 
replaced by a new phase transition into a novel thermo- 
dynamically stable phase with a finite smectic correlation 
length in both the high and low temperature phases. 



B. Summary, interpretation and consequences of the 
results 



The main conclusion of our work iS|-t|liat|— consistent 
with these and many other experiments J13'Ejl3 the three 
dimensional smectic phase, as defined by the existence 
of quasi-long-ranged translational order, is unstable to 
arbitrarily weak quenched disorder (i.e., arbitrarily low 
aerogel or aerosil density). Furthermore, we find that 
a new orientationally ordered low temperature "smectic 
Bragg glass" (SBG) phase replaces the smectic phase, 
and a thermodynamically sharp Nematic-to-SBG (N- 
SBG) transition, can survive in the presence of arbitrarily 
weak disorder, if and only if two universal positive defi- 
nite anomalous exponents rjs and rjK satisfy the bounds 



Vb + ^Vk > 4 . 



(1.1a) 
(1.1b) 
(1.1c) 



where the bounds in Eqs.l.la,b come from the require- 



ment of long-ranged orientational order and the condition 
for dislocations to remain confined, respectively. The re- 
gion in the 775, rjK plane that satisfies these three bounds 
is illustrated in Fig.l. 
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Fig.l. The region indicates those values of ry^ and rjx 
for which, in three dimensions, the long-range orienta- 
tionally ordered "smectic Bragg glass" phase is stable 
for sufficiently small disorder (sufficiently low density 
aerogel). 
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These exponents are universal; i.e., the same for all 
smectics in low density quenched disordered media. Un- 
fortunately, we have only been able to calculate them in a 
rather poor approximation: a d = 5 — e expansion. Since 
we are interested in d = 3, the ostensibly small param- 
eter e in this expansion is = 2, and so the expansion is 
expected to be rather poor. However, taking this expan- 
sion as our best estimate (since it is our only estimate) of 
the values of rjB and rjK, we obtain, in d = 3, ijb — 12/5 
and r]K = 2/5. These values, which are illustrated by the 
dot lab eled "e-expansion" in Fig.l, violate the bounds in 
Eq. |l.l| , and, hence, imply that the smectic Bragg glass 
phase does not exist, and that the thermodynamically 
sharp NA transition is destroyed by the presence of any 
disorder, no matter how weak. 

However, since the expansion parameter in this calcu- 
lation e = 5 — c? = 2 in the physical case of d = 3, this ar- 
gument that the transition and the SBG are destroyed is 
utterly uncompelling. Fortunately, experiments (which 
we will describe shortly) well below the pure transition 
temperature Tna can measure the exponents rjB and 
r]K- W ith these numbers in hand, the bounds given in 
Eqs.1.1 would then provide an unambiguous prediction 
as to whether or not the transition to the SBG is always 
destroyed by disorder. 

The physical significance of these exponents is quite 
intriguing: they reflect "anomalous elasticity" , which, in 
this context, refers to the fact that the smectic bend mod- 
ulus K and the layer compression modulus B are not 
constants in the randomly pinned smectic, but, rather, 
singular functions of the wavevector k under considera- 
tion. K{'k) and B(k) diverge and vanish respectively, as 
k ^ 0, according to the scaling laws: 



(1.2a) 



B(k) = B (A^^^D'" /^(fc.ef ■ (l-2b) 

Here kz and kj_ denote the projections of k along and 
perpendicular to the mean normal to the smectic lay- 
ers, respectively, and K and B without arguments k will 
here and hereafter denote the "bare" values of the smec- 
tic elastic moduli, i.e., their values in the pure (bulk) 
smectic. The unive rsal ex ponent ^ is determined by rjB 
and r]K through Eq. 1.12b| . 



The "non-linear" length scales S,^'" and are the 
distances in the _L and z directions respectively, beyond 
which the disorder -driv en anomalous elasticity is mani- 



fest. That is, Eqs.|L^ only hold for kj_^2 « 1 ^nd 



kzS,^^ << 1. If either of these limits is violated, both K 
and B ace-|k independent (up to logarithms), as in pure 
smectics.cj In three dimensions these nonlinear length 
scales are given by 



647r 



1/2 



^5/4 



1/2 
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(1.3a) 
(1.3b) 



Expressions for these lengths for general dimensional- 
ity d > 3 are given in Eqs. 16 and 6.8a of Sec.0. In 
the above is a measure of the strength of one of two 
types of disorder, (the other random "field" disorder Ay 
being less important in three dimensions) whose relation 
to various parameters of the aerogel is given in Sec.||. 
For now, it suffices to say that Ah is a monotonically 
increasing function of the aerogel density, and we expect 
it to be a smooth, analytic, non-singular, finite and non- 
vanishing function of temperature through the bulk NA 
transition temperature T/va- Since K is likewise well- 
behaved through Tna, Ci^^ is as well, a fact that we will 
make use of shortly. 

Although we have been unable to calculate the scal- 
ing functions fxix) and /b(x) exactly, we do know 
that they have the simple property of making Kih) 
and -B(k) independent of kz when the scaling argument 



opposite limit. This implies 



kz^z /if^±C± ) is << Ij and independent of fcj^ in the 



for fc.e"^ « {k±i±'')^ 
for fc,4^^ » (k^^^^)^ 



Brkloc^^'"' for fc.er « (fc^^i'^)C 
^Woc^ ^,B/C^ for fc.e'^ » (fc^^^^)C 



(1.4a) 
(1.4b) 



This strong wavevector dependence of K and B (driven 
by disorder) is caused by the same mechanism that leads 
to the much weaker (logarithmic) but still singular di- 
vergence and vanishing of K and B (driven by thermal 
fluctuations) that occurs in disorder-free smecticscfl: the 
anharmonic elasticity of large fluctuations in the smectic 
layers. 



As we will demonstrate in Sec. Vl, the effects are larger 
in the pinned smectic because the disorder induces layer 
roughness that is much larger than that due to thermal 
fluctuations, thereby leading to the stronger diverging 
and vanishing of K and B found here. 

The stability of the smectic Bragg glass phase depends 
on rjx and rys because the elastic moduli K and B deter- 
mine both the size of the orientational fluctuations and 
the stability of the phase against the unbinding of dis- 
locations. Requiring that real space orientational fluc- 
tuations remain flnite leads to Eq. l.la, while disloca- 
tions remain bound only if Eq.l.lb is satisfied. If the 
bounds Eq. 



1.1 



are satisfied for a real smectic liquid crys- 
tal, confined inside a low density aerogel, then the result- 
ing Pa — T phase diagram will be topologically identical 
to that displayed in Fig. 2. 

It is important to note that all of the results quoted 
here are equilibrium results. In contrast, in other pinned 
elastic media (e.g., Abrikosov flux lattices in dirty super- 
conductors and charge density waves in anisotropic met- 
als), it is often difficult to observe the true equilibrium 
behavior due to the extremely slow dynamics of those 
systems, which cause them to drop out|©|f equilibrium, as 
evidenced by strong hysteretic effects.0' 
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In this sense, smectics, and liquid crystals in general, 
appear to be far better systems for investigating the equi- 
librium effects of quenched disorder, since they typically 
do not exhibit hysteresis. This suggests that liquid crys- 
tals have intrinsically faster dynamics than, e.g., vortex 
lattices in superconductors, or charge density wave sys- 
tems in metals. 

That this should be so is hardly surprising: liquid crys- 
tals are, after all, "liquids" in the sense that, even in the 
translationally ordered phases like smectics, molecules 
are quite free to move around. In contrast, the atoms in 
the conventional crystalline solids are essentially locked 
into lattice sites, from which they can only escape by 
thermally activated hopping over substantial energy bar- 
riers. Even at room temperature, this hopping is ex- 
tremely slow; the fact that even "high-Tc" superconduc- 
tors are much colder than room temperature exacerbates 
the problem further. 

So, because of both their liquid-like microstructure and 
microdynamics, and the fact that they are at higher tem- 
perature than, e.g., superconductors, one might have an- 
ticipated that it would be far easier to observe the equi- 
librium effects of quenched pinning in liquid crystals than 
in other pinned elastic mediajStudied to date. This belief 
is supported by experiments,tZrEj and so we believe that 
the equilibrium results we obtain here should be directly 
testable in experiments. 
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Fig. 2. A schematic (aerogel density pA - temperature T) 
phase diagram for a smectic liquid crystal, confined in- 
side a low density aerogel, valid if and only if all of the 
bounds, Eq.1.1 are satisfied. As discussed in the text 
the "smectic Bragg glass" (SBG) and "nematic elastic 
glass" (NEG) phases are, respectively, translationally 
and orientationally disordered, but are topologically or- 
dered, and are therefore distinct from each other and 
from the fully disordered "isotropic" phase. The dotted 
line within the SBG phase represents a remnant of the 
disorder driven pinning transition of Fig^^ roun ded by 
the anharmonic elasticity studied in Secs. [v^[vill . 



The theoretical analysis that leads to these conclusions 



is quite interesting and novel. Our first-principles anal- 
ysis of smectics in aerogel demonstrates that the ran- 
dom pinning induced by the aerogel leads to only two 
potentially relevant types of random perturbations to the 
smectic: a random positional field disorder (hereafter re- 
ferred to simply as the "random field" and designated 
by Ay), which represents the aerogel's tendency to force 
the smectic layers to sit at particular positions, and a 
random orientational field disorder (hereafter called the 
"tilt" field, with designation A^), reflecting the aerogel's 
proclivity for particular orientations of the nematogens 
and smectic layers. If we ignore the anharmonic effects 
that lead to anomalous elasticity, we find that the re- 
sponse of the smectic-A phase to the random field dis- 
order in three dimensions is in very close mathematic 
analogy to that of the XY-model in two dimensions 
While this suggests that a nontrivial glassy phase might 
replace the smectic-A phase in analogy with the experi- 
ence with the 2d XY-model,El the presence of tilt disor- 
der (which is always generated by random field) actually 
destroys this phase. Even in the regime where the ran- 
dom field by itself has no effect at long scales, the tilt 
disorder leads to short-range smectic order parameter ip 
correlations, which fall off exponentially in the direction 
of the layer normal and as a Gaussian within the smec- 
tic layers (within a purely harmonic elastic model) . This 
absence of long-ranged order in the elastic model is a 
strong indication of its limitation and that dislocation 
defects are likely to proliferate. Focusing on the most 
important part of randomness, the tilt disorder, we find, 
in the approximation of ignoring elastic anharmonicities, 
that disorder always creates dislocations. The formalism 
we use to demonstrate this is new, powerful, and poten- 
tially applicable| to a wide variety of candidate "Bragg 
glass" systemsE^I. 

Once dislocations are present, the phase is best char- 
acterized as a nematic in a random tilt field. However, 
subsequent examination of orientational fluctuations in 
this nematic lead to the conclusion that tilt disorder de- 
stroys the orientational order of the smectic layers as 
well. Whether an orientationally and translationally dis- 
ordered low temperature phase distinct from the conven- 
tionally disordered high temperature isotropic phase is 
a "nematic elastic glass" (NEG), separated from it by 
a disclination unbinding traiwition, is the subject of an 
active current investigationHj 

Including anharmonic elastic effects considerably mod- 
ifies this picture. Although spatial dimensions d greater 
than three are clearly irrelevant to experiment, it is con- 
ceptually extremely useful (and fun!) to generalize our 
model to arbitrary spatial dimensions. We find that for 
d > 7, both the random field and the tilt disorder are ir- 
relevant, and a stable, translationally ordered, non-glassy 
smectic phase exists for sufficiently weak disorder. For 
5 < d < 7, tilt disorder remains irrelevant, but the ran- 
dom field and the anharmonic terms both become rele- 
vant, destroying smectic translational order and leading 
instead to a stable smectic Bragg glass phase. As in pre- 
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viouslsustudied "vortex glass" models of pinned elastic 
mediaE3iJ, here too we find that in all spatial dimensions 
5 < d < 7, real-space positional fluctuations diverge log- 
arithmically with system size: 



^(r)) =cdln [L^UXLJLD] 



(1.5) 



where i_L(z) is the linear spatial extent of the system in 
the _L (z) direction, 



A; 



(1.6) 



fu{x) is a universal scaling function and Cd is universal, 
dimension-dependent constant proportional to 7 — c? for 
d near 7. As a result, translational correlations through- 
out the range of spatial dimensions 5 < d < 7 decay 
algebraically, 



{ph {v)pG (0)) « r^^'^VdlAz/ri) 



(1.7) 

and 



with the exponent 77(d) a universal function of, 
fd{x) a universal (d-dependent) scaling function.c3 

Unlike the vortex glass case, however, here we find 
anomalous elasticity in this dimension range as well, ob- 
taining 



•A) 



if(k) cx |ln [k^fKih/Xk 
B{k) « |ln [k^fB{k^/Xkl)] 



|7K(rf) 

I 5 



(1.8a) 
(1.8b) 



where fK{x) and fsix) are universal scaling functions, 
and we have, rather remarkably, calculated the universal 
exponents jK{d) and 7_B(d) exactly for all d in this range, 
5 < d < 7, finding 



iKid) = - 
lB{d) = ^ 



1 f-d'^ + 12d- 23 



d2 -t- 6d - 
d^-l 
d2 + 6d - 13 
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(1.9a) 
(1.9b) 



Note that both the fluctuations of (^^(r)) c>c ln(L) 
and the anomalous behavior of -R'(k) and i?(k) are very 
weakly divergent functions of system size L and k respec- 
tively, depending only logarithmically on these quanti- 
ties. For d < 5, these weak logarithmic divergences are 
overwhelmed by power-law divergences caused by the tilt 
disorder, with K{\<i) and i?( k) d iverging and vanishing, 
respectively, according to Eq.|1.2|, with. 



^B = | + 0(e2) , 



where e = 5 — d, and 



{u\v)) = Ll^U 



(1.10a) 
(1.10b) 

(1.11a) 



with the universal roughness x ^-nd anisotropy C expo- 
nents given by 



X = 



C = 2- 



Vb + Vk 



(1.12a) 
(1.12b) 



Here fu{x) is another universal scaling function. 
This result Eg. 1.11 for (ii2(r)) leads to a quantitative. 



experimentally testable prediction for the X-ray correla- 
tion length f^, obtained by equating (^^(r)) = a^, where 
a is the smcctic layer spacing, and solving for L^, with 
Lj^ — !■ 00. This solution for is the X-ray correlation 
length that will be obtained by scattering off a powder 
sample (which probably means all samples, since, as dis- 
cussed earlier, there is probably no long-ranged orienta- 
tional order in d = 3, given our best estimates of rjB and 
riK)- The value of <^"^ so obtained is 



NL 



NL 



X «a 
X » a 



(1.13) 



This relation gives a way of experimentally measuring 
the exponents C and x- since the bare compression mod- 
ulus B(T) of the pure smectic is a strong function of tem- 
perature T near the bulk smectic NA transition tempera- 
ture TjvA, vanishing according to B{T) (x IT—TnaI"^, and 
since, furthermore, a, and K vary smoothly through 
Tna, a plot of In^"^ versus ln|T — Tna\ should yield a 
straight line for T near Tna, with a slope 0C/2x, provided 
that we stay far enough below T^a t hat A << a, so that 
the first expression for in Eq, 1.13 applies. In deriving 
this result, we have used the fact that A = KjB cx 

\T — Tn A l ^*^^^- This prediction for the X-ray correlation 



length, Eq. |1.13| , is schematically displayed in Fig. 3. 

A detailed renormalization group analysisS of the crit- 
ical behavior near the pure NA transition (detai ls of 
which will be presented elsewhere) shows that Eq. 1.13 
breaks down altogether at the temperature T^, which 
lies below (and for weak disorder very close to) the 
pure Tna- The corresponding reduced temperature <* = 
{Tna~T,)/Tna obeys 



^crit. ) ^ 



X 



(1.14) 



where is given by the second line of Eq. 1.13| , and we 
have defined a length scale ^crit.(^) derived entirely from 
properties of the pure system: 



^crit. (t) — ^'l 



2+0/21/^ 



(1.15) 



Lz^^^ ^ ior Lz/S,^'^ << {Lj_/(_^'")'' ,^ where £^,"^'^(t) is the X-ray correlation length within 



4^ fori,/e^»(i^/ci' 



the layers of the pure system at a reduced temperature 
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t = (T ~ Tna)/Tna above the pure transition tempera- 
ture Tna- Here z/_l is the X-ray correlation length expo- 
nent for the pure system in the _L direction, defined via 
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Fig. 3. Finite, temperature-dependent X-ray correlation 
length (T) for a smectic liquid crystal, confined inside 
a low density aerogel. The essential features of the X- 
ray correlation length that we predict are: (i) saturation 
as T ^ of S,^ (T) at a finite p^-dependent value for 
arbitrarily low aerogel density (diverging as pA 0), 
(ii) power-law scaling with t he disorder-free smectic bulk 
modulus B{T) (see Eq. |l.l3|) , (iii) crossover from the ex- 
ponent C/x to 2 in Eq. |1.13| as Tna is approached from 
below and anomalous elasticity becomes unimportant 
for , (iv) crossover to genuine critical behavior at the 
temperature at which ^crit.(7*) (indicated by the dot- 
ted curve below T^a) = ■ In the low density suerogels 
(e.g., PA = 0.08g/cm3) we find T^a - w 3°kM 

TheoreticallyEl, 4> is expectjed to be given by = 
vxY ~ 0.67. Experimentallyj23 the situation is more 
complicated, with (j) showing no universality, for reasons 
that are still unclear (at least to us). Hence, in extract- 
ing C/2x by the above analysis, the experimentally de- 
termined value of (j) of the particular bulk smectic that 
is confined in the aerogel should be used. One should 
also be careful, in such a fit, to treat Tjva as a fitting 
parameter, since the presence of even low density aerogel 
can presumable shift this temperature slightly. Alterna- 
tively, one could avoid these complications by studying 
the smectic outside the pure system's critical regime, but 
still reasonably close to the pure TAr^pSince in this regime 
mean- field theory applies and = l.Lil 

Once we know the ratio C/Xj we know C and x t hem 
selves, since the relation = 2 — x, implicit in Eqs.1.12 



provides us with a second equation for the two unknowns 
C and X- If the roughness exponent x so obtained is > 1, 
we predict that orientational fluctuations diverge, and 
hence, long-range orientational order is destroyed. As a 
result, there will be no "smectic Bragg glass" phase, and 
hence, no thermodynamically sharp transition associated 



with the pure NA transition that occurs i n the absence of 
disorder. In any event, the prediction Eq. |l.l3| for A >> a 
will inevitably apply sufficiently close to the pure T^a, 
since B{T Tna) 0, while K{T Tna) remains 
non-zero, and, hence X(T — > T^a) ^ oo » a. 

Should it prove that orientational order is possible in 
d = 3 (i.e., that the d = 3 values of the universal ex- 
ponents TjB and rjK do satisfy the inequalities in Eq 
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which would imply that the SBG phase is stable), then 
single orientational domain scattering zs, in principle, 
possib le. I n such scattering, given by as given 
by Eq. 1.13| , will be the correlation length along the mean 
normal to the smectic layers. The correlation length 
along the smectic layers is then given by 



ef = 



NL 



fNL 



1/x 



A << a 
A >> a 



(1.16) 



These X-ray correlation lengths should not, however, 
be interpreted as being the length scales beyond which 
smectic correlations cease. In fact, smectic behavior, by 
which we mean, e.g., the anomalous smectic elasticity, 
Eq.1.2, persists out to much longer lengths. Specifically, 
in the _L direction, it persists out to the shorter of the 
two lengths and S^^. £_o is the distance over which 
orientational order would be well-correlated in the ab- 
sence of dislocations, while ^j' is the distance (in the _L 
direction) below which dislocations would remain bound 
in the absence of large orientational fluctuations. In ac- 
tuality, only the prediction of the smaller of the two of 
these lengths is valid. We expect that when << 
the large orientational fluctuations also induce disloca- 
tion unbinding at about the same length ^o- On the 
other hand, in the opposite (^o >> £,±), limit the dislo- 
cation unbinding occurs first at and on longer scales 
our system is indistinguishable from a nematic in a ran- 
dom orientational field. Whether the orientational order 
is stable or not in this case must be re-evaluated withiu. 
an effective random-field nematic model using a LarkinQ 
type of analysis. c3 

The "orientational correlation length" is given by 



NL 



NL 



A 



(1.17) 



while the "dislocation length" scale in the _L direction, 
is given by 



_ cNL I '^^ 



■2/:NL \ 2{r,^-i) 



\cK^a^d^ 



(1.18) 



where d is a microscopic length of order the layer spac- 
ing a, and c{T) is a dimensionless constant. c{T) van- 
ishes like 6"^"=/-^ as T ^ and diverges like T/Ec{T) cx 
\Tna — T\~^^ , where Ec is the "core energy" of a dis- 
location line segment of length d, and 7i- is the "line 
tension" critical exponent describing how Ec{T) vanishes 
as T ^ Tna from below. We expect the resulting diver- 
gence of c(r) to overwhelm the corresponding divergence 
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of A(T), leading to a that gets smaller as T — > Tna 
from below. Clearly then, sufficiently close to Tna, 
gets to be less than the system enters the strong- 

disorder regime, and the weak-disorder theory we have 
presented here no longer applies. 

In the 2;-direction, the orientational correlation length 
is also given by while the dislocation length is given 
by 



NL £X tlD £ 
Z ^Z *^Z ''O 
1 1 1 ► 

Lz 



Fig. 4. One possible hierarchy of important length scales 
in the problem of a three-dimensional weakly disordered 
smectic, valid for A << a, x > 1, and riK > 1. If instead 
X < 1 and rjK < then the putative SBG phase is stable 
and and S,^ ^ are infinite. 

and for the case of A >> a in Fig. 5 



; NL 
)Z 



^Z ^0 



Note that the orientational correlation length given by 
Eq. 1.17 is always much greater than the non-linear length 
S^^-'^ in the weak disorder limit, in which ^^'^ oo. 
However, if we hold the disorder strength fixed (i.e., hold 
A/j fixed), and approach the pure NA transition from 
below, we will always, eventually at some temperature 
T* < TjsfA (T* Tma in weak disorder limit) leave this 
weak disorder regime (no matter how weak), since 



as given by Eq. l.3b| , does not change much as T — > Tjy^, 
while A(T ^ T"^) ^ oo. So the ratio Co/C±^ de- 
creases without bound as T ^ ^jvA' ^"^d, hence even- 
tually drops below 1, signaling our entry into the strong 
disorder regime for T > (see Fig. 3). This is not sur- 
prising: the smectic's resistance to the perturbing effects 
of the disorder is provided, in part, by B; so when i? — > 0, 
as it does as the pure NA transition is approached, any 
disorder will eventually look "strong" , and our theory will 
no longer apply. As long as we are at temperatures T far 
enough below Tjva, however, our weak disorder results 
will apply. As we weaken the disorder, by e.g., reducing 
aerogel density, we can apply our weak disorder theory 
closer to T^r^. 



The expression (1.17) for the orientational correlation 
length only holds if the condition for the stability 
of the long-ranged orientational order, Eq. 1.1a , which is 
equivalent to x < 1: is violated (i.e. Eq.1.17 only holds 
for Y > 1). If X < 1, Co is infinite. Likewise, Eqs.l.lS 
and 1.1£ for C;^^ only hold if ry^f > 1 (i.e., if disloca- 



tions are, in fact, unbound); otherwise, ^ are infinite. 
Clearly, if both bounds x < 1 a-nd -qx < 1 are satisfied, 
all three lengths are infinite, smectic behavior holds out 
to arbitrarily large length scales, and the smectic Bragg 
glass is a stable phase of weakly disordered smectic liquid 
crystals. 

A summary of the many length scales in our theory for 
three dimensions and their hierarchy, in the weak disor- 
der fimit, is illustrated for the case of A << a in Fig. 4 



Fig. 5. The hierarchy of important length scales in the 
problem of a three-dimensional weakly disordered smec- 
tic, valid for A >> a, X > 1, and t/k > I- If instead x < 1 
and rjK < 1, then the putative SBG phase is stable and 
Co and ^ are infinite. 

Just as the proposed "Bragg glass" phase of randcuxily. 
pinned elastic media is distinct from the liquid phase,Lrtij 
here, too, the absence of long-ranged orientational order 
does not preclude a "nematic elastic glass" phase that 
would be separated from the high temperature "isotropic 
fiuid" phase by a thermodynamically sharp equilibrium 
phase transition. That such a nematic elastic glass phase 
may indeed ^ist is suggested by dynamic light scattering 
experiment£3'E3 that show a dramatic slowing down of 
director fluctuation relaxations in liquid crystals in aero- 
gel below a temperature Tneg near the bulk NI transi- 
tion. 

We are currently investigating theoretically the possi- 
bility of such a nematic elastic glass phase, and initial 
results indicate that such possibility is, in fact, allowed 
theoretically. Our reailts on this subject will appear in 
a future publicationE3 

The low temperature analysis presented here is fur- 
ther supported by a complementary approach that in- 
vestigates the stability of the NA transition to disorder 
from the high-temperature nematic phase. l2I The results 
of this latter work substantiate our findings that in three 
dimensions the NA transition and the smectic-A phase 
are destroyed. The rounded remnant of this transition 
can, however, be studied utilizing self-consistent meth- 
ods previously successfully applied to pure systems below 
the lower critical dimension.L3 This method predicts a 
Lorentzian-squared structure function with a correlation 



7 



length that monotonically, slowly increases through the 
bulk NA transition, in good agreement with the X-ray 
measurements » The calculated specific heat-is also in 
qualitative agreement with the experiments ,EjIEI exhibit- 
ing a broad, well-defined peak at a temperature slightly 
lower than the bulk transition temperature Tjva- 

As we hope the introduction made clear, the phe- 
nomenology of this system is extremely rich. Our orga- 
nizational approach to presenting the derivation of these 
results is the following: starting from the full model, we 
first throw out all but the simplest effects; only after de- 
veloping a full understanding of this simplified model do 
we, one at a time, reintroduce the complicating effects 
we had initially thrown out, building gradually toward 
the complexity of the full physical model. 

Since many of these complicating effects prove to be 
important, this approach has the drawback that some of 
the results derived for the simplified models do not, in 
fact describe our physical system. We have attempted to 
alert the reader everywhere results of a simplified mod- 
els differ from our actual predictions for the full theory. 
Only the latter, of course, should be compared directly 
with experiments. 

The remainder of this paper is organized as follows: 
in Section ||, we introduce and motivate our model for 
smectics in disordered media, and discuss those aspects 
of the model specific to aerogel. 

In Section III, we derive an elastic low temperature 
model of a randomly pinned smectic, at first ignoring an- 
harmonic elasticity and dislocations. Using renormaliza- 
tion group methods, we study the stability of the smectic 
phase within this model and these approximations. 

we use the results of this RG treat- 



In Section IV 



ment to calculate various smectic correlation functions, 
still treating the elastic theory to harmonic order and 
ignoring dislocations. 

Section ^ incorporates the dislocations in the har- 
monic theory, and shows that, in that theory, in three 
dimensions, dislocations always proliferate, even in the 
presence of arbitrarily weak disorder, thereby destroy- 
ing, within the harmonic approximation, both the smec- 
tic Bragg glass phase and the thcrmodynamically sharp 
transition to it. 

The effects of the previously ignored anharmonic elas- 
tic terms are shown to be very important in Section 



VI, which develops a renormalization group treatment 
of them in a model that neglects both the random field 
and dislocations, but includes tilt disorder. The anoma- 
lous elasticity Eq 1.2 and the u — u correlation functions 
that follow from it are derived in this section. Given 
the importance of these elastic nonlincarities, we study 
their effects on dislocation unbinding and orientational 
order in Section VI] by incorporating anomalous elastic- 
ity into the duality theory de velop ed in Sec tion |v|. We 
thereby derive the bounds Eq.(l.la) and Eq.( |l.lb| ) on the 
anomalous exponents 773 and r]K for the stability of the 
smectic Bragg glass phase. The irrelevance of the random 
field in a full, anharmonic theory below = 5 is demon- 



strated by a fun ctional renormalization group treatment 
in Section VIII , in which we also calculate exactly the 
exponents 7b (rf) and "fxid) that govern the anomalous 



elasticity for 5 < d < 7. In Section IX, we demonstrate 
that in d — 3, dislocation unbinding is not induced by 
the random field alone. 

We conclude the main part of the manuscript with 
Sec.^ where we discuss interpretation of our results in 
terms of past experiments and their implications for the 
future experiments, as well as many remaining interesting 
and important theoretical problems, some under current 
investigation. 

The details of the analysis of the random field disorder 
in three dimensions, the functional renormalization group 
analysis in higher dimensions, an alternative derivation 
of the dislocation loop theory, and a fluctuation-corrected 
mean-field treatment of the dual model of randomly 
pinned smectic are presented in Appendices A, B, C and 
D, respectively. 



II. MODEL 

Our theory of the disordered NA transition is based 
on the de Gennes model. Near the mean-field transition 
from the nematic to the smectic-A phase, the center-of- 
mass nematogen molecular density p{r) (which is liquid- 
like in the nematic phase) begins to develop strong fluctu- 
ations dominated by Fourier components that are integer 
multiples of the smectic ordering wavevector qo = n27r/a 
(a is the layer spacing) parallel to the nematic director 
n. We take the dominant lowest Fourier component ip{r) 
as the local (complex scalar) order parameter which dis- 
tinguishes the smectic-A from the nematic phasalj. It is 
related to the density p(r) by 



p(r) = Re[po + e''i°-""^(r)] 



(2.1) 



where po is the mean density of the smectic. 

As first discussed by de Gennesllj, the effective Hamil- 
tonian functional HdoW'jf'] that describes the NA tran- 
sition at long length scales, in bulk, pure liquid crystals, 
is: 



c±|(V_L-%'5n)7/.|2 + c,||V||?M 



io|V'P+25ol^l' 



Hp[n] , 



(2.2) 



where oc [T — T^^) /T'^^ , ^jvA° transition 
temperature in the pure (bulk) system, 5n(r) = n(r) — no 
is the fluctuation of the local nematic director n(r) away 
from its average value fio, which we take to be z, sub- 
scripts II and _L denote the directions parallel and trans- 
verse to rio, and 7f_F[n] is the Frank effective Hamiltonian 
that describes the elasticity of the nematic order director: 



8 



+Kb{n X V X n)" 



(2.3) 



where Kg, Kt, and Ki, are the bare elastic moduh for 
splay, twist and bend of the nematic director field, re- 
spectively. 

The minimal coupling between n and ip is enforced 
by the requirement of global rotational invariancal3. It 
is important to emphasize, however, that although the 
de Gennes Hamiltonian Hdc is closely analogous to that 
of a superconductor, there are essential differences. The 
physical reality of the nematic Sn and the smectic ■0 or- 
der parameters, in contrast to the gauge ambiguity in the 
definition of the vector potential and the superconckujts. 
ing order parameter, selects the liquid crystal gaugeEflS 
i5n • rio = (since |nop = 1) as the preferred physical 
gauge. The strict gauge invariance of Hdo is already ex- 
plicitly broken by the splay term ifs(V • n)^ of the Frank 
Hamiltonian. This is one source of distinction between a 
smectic-A liquid crystal and a superconductor. 

In fact, since the only true symmetry of the smectic-A 
liquid crystal is invariance under global simultaneous ro- 
tation of the smectic layers (equivalent to V' ^ ipe^^''^^ ) 
and the nematic director Sn ^ 5n + 9 accomplished by 
a constant (in space) 6 — const., and not arbitrary local 
gauge transformations, the requirement of gauge invari- 
ance severely overrestricts the allowed effective Hamilto- 
nian. That is, the de Gennes model is more symmetric 
than is required by the physics of the NA transition, and 
a more general effective Hamiltonian will allow for terms 
that break nonlinear gauge invariance. All of the exten- 
sion terms that can be added to de Gennes model are, un- 
fortunately, irrelevant (which is probably why de Gennes 
did not include them in his original formulation). How- 
ever, in light of the above remarks, the bending stiffness 
Kg in the smectic phase, (i.e. the coefficient of (V_lu)^) 
can differ from the corresponding coefficient in the ne- 
matic phase (the coefficient of (V • n)^), due to the pres- 
ence of terms like, e.g., jV^-i/'p, in the smectic Hamilto- 
nian, which violates local gauge invariance, but not global 
rotation invariance. This is in contrast to what is usually 
assumed in all the standard treatments, of this problem, 
which take Kg and Kg to be the same.Ej The importance 
of this observation is that in mean-field theory it leads to 
a slope discontinuity in Ks{T) at the T^^, which is-likely 
to persist in the theory that includes fluctuationsEj. 

Guided by these observations, we now proceed to con- 
struct the disorder part of the effective Hamiltonian. 
Because ip{r) is related to the smectic density p{r) via 
Eq. p.l] , the randomness can couple to it directly. In par- 
ticular, symmetry allows a coupling of the form 



l6t{r){p-por 



U{v)p 



(2.4) 



where both the random Tc {5t{v)) and the random po- 
tential U (r) are proportional to the local aerogel density 



Pyi(r). We will take this aerogel density to be a quenched 
random variable, and denote averages over it by a hori- 
zontal over-bar. 

This means in particular that the autocorrelations of 
(5i(r) and U{y) will be proportional to those of tbe aero- 
gel, and, hence, long-ranged for fractal aerogel.E3 That 
is, we expect 



U{r)U{v') = Cu{r-v') 



Tupa{y)pa{v') , 



5t{Y)6t{v')=Ct{v^Y') , 



TtpA{v)pA{v') 



(2.5) 
(2.6) 



where Vu and Tt are coupling strengths that depend on 
the microscopic physics of the aerogel-smectic interac- 
tion, but not, at low aerogel density, on the aerogel den- 
sity itself, nor on the smectic density. 

The fractal structure of aerogel over a range of length 
scales can be determined from the aerogel density-density 
correlation function pa{y)pa{Q) measured in an X-ray 
scattering experiment. Writing the aerogel density pA 
in dimensionless (volume fraction) units, the correlation 
function p^(r)p^(r') measures the probability that, given 
that a point r is occupied, the point r' is also occupied. 
The probability that an arbitrary point r is occup ied is 
clearly proportional to the aerogel density pA^Lf), mea- 
sured at (averaged over) a scale beyond which the aerogel 
ceases to be fractal. We call this scale Lf. The aerogel 
is also not a fractal for lengths smaller than the micro- 
scopic length a/ of the order the aerogel strand diam- 
eter. The density of a structure which is a fractal on 
scales Of < r < Lf, with fractal dimension dp, is, by 
definition, given by the ratio of the number of occupied 
sites N ^ [Lf /of Y'' in a volume V = L'j to this volume. 

Hence, the average aerogel density pA = pA{Lf ) is given 

by 



PA{Lf) « ( 1^ 



(2.7) 



Given that a site r is occupied, the conditional proba- 
bility P(r,r') that a site r' is also occupied is, roughly, 
the typical "mass" (actually total volume) of material 
M(|r — r'l) contained within a sphere of radius |r — r'| 
centered on a point r that is also on the fractal, divided 
by the total volume |r — r'j'' of that sphere. By definition, 
M{\r — r'l) « a'j{\r — r'|/a/)'^^ (this is what we mean by 
the fractal dimension dp). Hence, the conditional proba- 
bility is P(r, r') = Af(|r-r^/|r-r'|'* - (a//|r-r'|)''-''^. 



Combining this with Eq.2.7, we obtain 



PA{r)pAir') « ( 



PA 



d—dp / \ d—dp 



r — r 



|r — r' 



(2.8a) 
(2.8b) 
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We will see in a moment that these long-ranged correla- 
tions have no effect on the long distance behavior of the 
smectic. 



Combi ning the random field energy Eq. 2.4 with the re- 
lation Eq,2.1 between the smectic order parameter ^ and 
the density p, we obtain: 



St{r)\il;\'^ + V{r)i; + V*{r)'ilj'' 



(2.9) 



where we have defined the complex random potential 
V{r), which acts on the NA order parameter ip like a 
random magnetic field acts on a spin, and is related to 
the potential disorder U (r) by 



V{r) EE C/(r)e 



iqoz 



(2.10) 



Note that, despite the long-ranged correlations of U{r), 
V{r) has only short-ranged correlations. To sec this, con- 
sider a double Fourier transform of the V{r)V*{r) corre- 
lation function 



y(k)y*(k') = J d'^rd'^r'e-'^''e-'^' '' Vir)V*ir') , 

= (27r)'*,5'^(k + k')Cy(k) , (2.11) 

where 

Cv{k) = J d'^ir - r')e'(k+9o2)-(r-r')[;(j.)[/(j.,) ^ 

= Ccr(k + qo2), (2.12) 

and Cjj ( k) i s the Fourier transform of the U —U correla- 
tion, Eq j2.5| . Even though we expect this Fourier trans- 
formed correlation function Ct/(k) to diverge as k ^ 0, 
due to the power-law spatial correlations in U (r) , there 
is no reason for it to diverge as k ^ q^i, since the aerogel 
itself has no particular spatial structure at the wav evec- 
tor of the smectic ordering. Hence, we see from Eq. 2.12 
that the V —V autocorrelation function remains finite as 
k ^ 0. Thus, the correlations of T^(r) are short-ranged, 
and hence we can accurately capture the long distance 
physics of the problem by taking those correlations to be 
zero-ranged, and writing 



V{v)V*{r') = W5'^{y-y') , 



where 



d—df 



afqo 



dp 



(2.13) 

(2.14a) 
(2.14b) 



as can be readily seen by Fourier transforming Eqj2. 

I n a ddition to the density coupling of Eq.2.4 and 
Eq.^.9|, we expect that the long nematogenic molecules 
that make up the smectic will tendi to line up with the 
randomly oriented aerogel strands ,E3 whose orientational 



correlations we expect to be short-ranged. This interac- 
tion leads to an additional orientational random coupling 
of the form: 



d^risii 



(2.15) 



where the quenched field g(r) is random and short-range 
correlated in direction, 



5,(r)g,(r')=A,<5'^(r-r')<5, 



(2.16) 



with strength Ag proportional to the local aerogel den- 
sity. 

As discussed in the Introduction, motivated by the ex- 
perimental observations, we assume that near the NA 
transition, in low-density aerogel samples, the nematic 
order is well-developed. In this regime, we can consider 
small disorder-driven deviations from perfect nematic or- 
der fio = z as small, writing n(r) as 



n(r) = z-\/l — |(5np + 6n{r) . 



(2.17) 



Substituting this representation inside Hdn, Eq. ^.l5 



and, given that \Sn\ 1, keeping only up to quadratic 
terms in 5n, and dropping (5n-independent terms, we ob- 
tain 



H, 



d7i 



2h(r) • Sn 



g,(r)2|5np + (g(r)-<5n)^ 



(2.18) 



where we have defined a quenched random tilt field 

h(r) = g,(r)g(r) . (2.19) 

It is easy to see that for isotropic disorder g(r), the terms 
quadratic in 5n(r), above, cancel each other on average 
and therefore only make unimportant contributions that 
are beyond quadratic order in (5n(r). As a consequence, 
for such isotropic disorder, the effecti ve or ientational dis- 
order is given by the last term in Eq.2.18 



Hdn ~ J d'^r h(r) • 6n 



(2.20) 



We expect that random field h(r) is related to the aerogel 
density structure by the following correlation function 



= np\{r)p\{r'%{rMr)t.{r%{r') , (2.21) 

where t(r) is the local tangent to the aerogel fibers, and 
F/i is a p/i-indepcndent coupling constant that depends 
on the microscopic aerogel-smectic interaction. 

Since we expect this tangent t(r) to have only short- 
ranged correlations (with range of order the orientational 
persistence length of the silica fibers), the above corre- 
lation function of the tilt disorder should also be short- 
ranged. Furthermore, it must be isotropic. These con- 
siderations, taken together with the fractal nature of the 
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aerogel, and Eq.2.21, lead to the following form for the 
correlatfon function of the random tilt-disorder h(r) 



= AhS^ir - r') , (2.22) 

which is short-ranged. This defines the tilt field disorder 
variance A;;. 

Thus our model of the NA transition in aerogel is char- 
acterized by the effective Hamiltonian functional H = 
HdG + Hdp + Hdn, with HdG, Hdp, and Hdn given by 
Eqs.(Ul), (H^l), and ( |2.20D , respectively. This total ef- 



fective Hamiltonian must be supplemented with corre- 
lation functions for the random Tc {6t(r)), random field 
{V{r)), an d random tilt ( h(r)) diso rders, which are given 
by Eqs. ^^ , |2.13| - ^.14i and ^.22i respectively. 



Finally, it is essential to keep in mind that, as discussed 
above, even for fractal aerogel, although the correlations 
of the random Tc field St{r) are long-ranged (power-law 
correlated), those of the random field V{r) and the tilt 
h(r) are short-ranged, with W « A^. 



III. SMECTIC PHASE AND ITS STABILITY 
WITHIN THE HARMONIC ELASTIC MODEL 

In this section we study the disordered NA model, 
defined in the preceding section, within the low- 
temperature phase. While it is tempting to directly 
analyze the model H — HdG + Hdp + Hdn written in 
terms of the smectic order parameter if), we will not do 
so here. Our motivation for this is two-fold: (i) As dis- 
cussed in the Introduction, even for the bulk (disorder- 
free) case the results obtain through such a direct (de 
Gennes model) approach, have, so far, failed to give pre- 
dictions that agree with experiments even on the hulk 
NA transition. Having not fully understood the diffi- 
culties with the bulk NA transition, we hesitate to use 
such a high temperature ( "soft-spin" ) analysis to study 
the significantly more complicated NA transition in the 
presence of disorder (i.e., confined inside the aerogel), (ii) 
While such a "soft-spin" approach is often superior in un- 
derstanding the critical properties of the transition and 
the high temperature phase, it is significantly less suc- 
cessful in the analysis of the low temperature phase. For 
example it is known to incorrectly predict the lower crit- 
ical dimension in, e.g., the random field Ising model, as 
well as completely missing the existence of the Kosterlitz- 
Thouless transition in the XY-model. Since one of the 
main goals of the current initial investigation of the dis- 
ordered smectics problem is to ascertain the stability and 
nature of the low temperature phase, we leave the direct 
de Gennes rnojdel approach to this problem to subsequent 
publications^. 

Instead, here, as a first analysis of the problem, we 
choose the elastic, low-temperature approach in terms of 
the Goldstone phonon mode u. The predictions of such 



an approach for the ordered pure smectic phase are free 
of controversies, and are in agreement with experiments 
on bulk smectics. 

We will not, in cither this or the next section (Section 
[V), treat the full elastic model in all its complexity, but 
rather, begin by simplifying the model by ignoring elas- 
tic anharmonicities (i.e., terms higher than quadratic in 
gradients of u). We will also, in both this and the next 
section, neglect the effects of dislocations, and focus on 
three spatial dimensions. All of these restrictions will be 
removed in later sections of the paper, as we build up to 
the full complexity of the disordered smectic system. 

We begin by assuming the existence of smectic order 
and investigate if/when this assumption is violated be- 
cause of the interaction of the smectic with the random 
environment of the aerogel. 

Within the ordered smectic phase, the fluctuations are 
conveniently described in terms of the fluctuations of the 
magnitude and phase of V'- It is easy to show that the 
fluctuations of the m agnitude of ip around the average 
value 1^0 1 = 1/^0/50 = const, are "massive", and can 
therefore be safely integrated out of the partition func- 
tion, leading to only finite., unimportant shifts in the ef- 
fective elastic moduli. In contrast, the phase of ■0 is a 
U(l) massless Goldstone mode, corresponding to sponta- 
neously broken translational symmetry. It is the essen- 
tial low energy phonon degree of freedom of the smectic 
phase, describing the local displacement of the smectic 
layers from perfect periodic order. In accord with this 
discussion, deep within the smectic phase, we can repre- 
sent the smectic order parameter as 



V'(r) = IV'ol 



(3.1) 



safely ignoring (actually integrating out the "massive" ) 
fluctuations in the magnitude lipol oiip. It is important to 
note that this can be done at any temperature below the 
transition, without any qualitative consequences for phe- 
nomena occuring on sufficiently long length scales, larger 
than a well-defined crossover length ^*(T). The elastic 
model is then rigorously valid on length scales larger than 
(T) and breaks down on shorter scales, and therefore, of 
course can only make predictions about phenomena (e.g., 
length scales like larger than C*(r). As T ^ T^^, 
^* (T) diverges and the range of length scales about which 
the elastic model is able to make predictions shrinks, be- 
ing pushed out to infinite scales. Consequently, for exam- 
ple, sufficiently close to T^^, i.e. for T > T* (see Fig. 3), 
such that ^* from above, our prediction for , 

based on the weak disorder elastic theory is no longer 
strictly valid. 

Using this low-temperature ansatz inside the tp- 
dependent part of the effective Hamiltonian, given by 



Eqs.2.2 and 2.9 and dropping constant terms, we find 



H[u,Sn] 



d^r 



-^(V • Snf + ^{z ■ V X Snf + ^{z x V x Snf 



2 
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+h(r) • Sn~ \ipo\U{r)cos[qo(u{r) + z)] 

(3.2) 

where Bj^ = cj^lV'oP'Zo ^^'^ ^ = C|||V'op9o- We observe 
that the fluctuation mode (V_lw — 6n) is "massive" and 
leads to the Anderson-Higg's mechanism, a hallmark of 
gauge theories. As a consequence, after a simple Gaus- 
sian integration over Sn, we find that at long length 
scales, dn fluctuations are constrained to follow 'V±u. 
The resulting effective elastic Hamiltonian is then ob- 
tained by the replacement 



Sn — * V±u , 



(3.3) 



everywhere in the Eq.3.2. This is valid in the long wave- 
length limit, to quadratic order in gradients of u, and 
provided dislocations are confined. We obtain. 



H\u] = / d 



K 
"2 



h(r) • V±u 



IV'o|C/(r) cos[(7o(u(r) + z)] 



(3.4) 



or, equivalently, in terms of V{t) 
H[u] = j d'^r 

-\i>o\ (y(r)e**"('") -I- y*(r)e-**"("")) 



(3.5) 



In the above K = Kg, although, as discussed in the be- 
ginning of Secjl|, in a model of the NA transition that is 
more general than the de Gennes model, K and Kg can 
differ by a singular function of the reduced temperature 
\T-Tna\- 

To compute self-averaging quantities, (e.g., the disoi.- 
der averaged free energy) we employ the replica "trick"E3, 
which allows us to work with a translationally invariant 
field theory at the expense of introducing n replica fields 
(with the n limit to be taken at the end of the calcu- 
lation). For the free energy this procedure relies on the 
identity for the ln(x) function 



-TlnZ 



-T lim 



Z" - 1 



(3.6) 



After replicating and integrating over the disorder using 
Eqs. |2. 13)^1^ , ^2|, we obtain 



/ [dua]e 



(3.7) 



The effective translationally invariant replicated Hamil- 
tonian H[ua\ is given by 

H[u„] ^ j d'r (f(Vi^a)' + |(9.«a)' 



Ay COs[go(wQ - Uf})] 



a,/3=l 



(3.8) 



where Ay = |-0opVr- The statistical symmetry under 
global rotation forces the disorder generated replica off- 
diagonal terms to be invariant under Ua{v) — > UnW) + 
6 ■r±. In the replicated effective Hamiltonian Eq. 3.S the 



nonlinearities only depend on the difference between dif- 
ferent replica fields and therefore do not depend on the 
"center of mass" field X]"=i """j which is therefore a non- 
interacting field. This implies an exact result that K and 
B are not renomialized by the disorder in this harmonic 
approximation.p A;i, of course, will be renormalized by 
the random-field nonlinearity. 

It is easy to see that the random tilt disorder term with 
coefficient A/j can be rewritten in Fourier space as 



Q,/3 



which, in the n - 
the Hamiltonian 



a,/3 



(3.9) 



limit leads to the quadratic part of 



1 r " 

Q,/3 



(3.10) 

from which the propagator Gap{(l) defined through 

{u^{q)up{qf)) = G„/3(q)<5'^(q + q/) (3.11) 
can be easily obtained, 

^"^(^^ - Kqi + Bql + [Kqi+BqlY^ ^^'^'^ 
using an identity for inverting matrices of the type 

Aap^a5ai3 + h, (3.13) 

namely: 



~ a a{a + bn) ' 

n^O —OaB ^ • 



(3.14) 



In three dimensions, for Ah = Ay = (disorder- 
free liquid crystal) at low temperatures (T < Tna), the 
smectic-A phase is described by a fixed plane, defined by 
the bare values of K and B. Our initial goal is to estab- 
lish how this fixed plane is destabilized by the disorder, 
i.e. by the random field and random tilt terms. The cal- 
culation ispargeneralization of that for the 2d random-field 
XY modelBH to the anisotropic elasticity of the smectic-A 
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in three dimensions. We employ the standard momentum 
shell renormalization group transformatiorCJ, by writing 
the displacement field as 



Ua{r) 



:(r) 



(3.15) 



integrating perturbatively in Ay the high wavevector 
part M^(r), nonvanishing inside a thin cylindrical mo- 
mentum shell 



Ae < |q_L|< A , 

— OO < Qz < CO , 



(3.16a) 
(3.16b) 



terms and topological defect loops in u; the latter be- 
come important at high temperatures, where they in- 
duce th£_lsLA transition, by driving the bulk modulus B 
to zero.EjH It is convenient to choose the anisotropy ex- 
ponent u! — 2, becau se such a choice keeps K{£) and 
B{e), Eqs. ^.l8b| , ^l8^, fixed under the RG. Ahhough it 
appears from Eq.3.18a that the glass transition temper- 
ature (relevance and irrelevance of Ay) depends on the 
arbitrary choice of it does not. This choice is actu- 
ally completely arbitrary and does not affect any physi- 
cal quantities. This can be seen by looking at the proper 
dimensionless coupling constant 



and rescaling the lengths and long wavelength part of 
the fields with 



z 

M<(r) = e 



V(r') 



(3.17a) 
(3.17b) 
(3.17c) 



so as to restore the ultraviolet cutoff back to A. Because 
the random-field nonlinearity is a periodic function, it is 
convenient (but not necessary) to take the arbitrary field 
dimension (j) — Q, thereby preserving the period 2sc./qa un- 
der the renormalization group transformation.eil Under 
this transformation the resulting effective Hamiltonian 
functional can be restored into its original form Eq.3.8 
with effective ^-dependent couplings. We relegate the 
details of these calculations to Appendix A, and focus 
here on the results. These can be succinctly summarized 
in the renormalization group flow equations 



M 
dKje) 

d£ 
dB{£) 

de 

de 



(2 + t^-77)Ay - AiA^ , 
(lo - 2)K , 
(2-w)S, 
LoAn + M^l , 



where we have defined 



SttA^ViOB ' 
egg 

47rAVA'3B ' 



(3.18a) 
(3.18b) 
(3.18c) 
(3.18d) 

(3.19a) 
(3.19b) 
(3.19c) 



A, 



AiA, 



(3.20) 



whose recursio n flow equat ion can be eas ily obtained by 
combining Eqs. 3.18a 3.18b, and Eq. 3.18c 



dAyje) 
di 



(4 - r;)Ay 

[4-77- Ay] Ay . 



(3.21) 



Obviously this flow equation is independent of the arbi- 
trary rescal ing ex ponent uj and has the same form as that 
for Ay, Eqj3.18a|, with uj = 2. 



From Eq.3.21, wc then find that for 77 < 4 (large elastic 
moduli and low temperature), or, equivalently, below the 
pinning transition temperature Tp, given by 



9^ 



(3.22) 



the smectic (Ay = 0) fixed plane is unstable to disorder. 
However the initial runaway of disorder is halted by the 
nonlinear terms in Ay, which terminate the flow at a 
new flnite disorder flxed line. 



Ay(T) 



(3.23) 



and c is a dimensionless number of order 1. Note that 
dT]{i)/d£ = 0, exactly. 

As discussed abov e, t he symmetry of the effective 
Hamiltonian H in Eq.3.8 guarantees that the flow equa- 
tions for K{£) and B{£) are exact (arising from simple 
length rescaling with no diagrammatic corrections), ig- 
noring (for now) the effects of both anharmonic elastic 



This new flxed line then describes a randomly pinned, 
glassy smectic-A phase, analogous to the 2d suaer-rough 
phase of crystal surface on a random substratea and the 
1-f 1 vortex glass phase of flux-lina-vortices (conflned to a 
plane) in type II superconductoraJ. This RG flow struc- 
ture is summarized in Fig. 6. 
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Fig. 6. A fixed line characterizing the pinned glassy 
phase of a three dimensional smectic in aerogel. Tp — 
Wiry/ KB/qq is the pinning transition temperature, tak- 
ing place within the 3d SBG phase. This transition only 
survives within the harmonic elastic approximation. 

It is important to note that the transition at Tp, il- 
lustrated in Fig. 6, is taking place within the SBG phase. 
That is, Tp is a transition from the (T > Tp) SBG phase, 
within which the translational random field disorder Ay 
is unimportant at long scales, to the (T < Tp) SBG phase, 
in which Ay is relevant and flows to a nontrivial fixed 
line. This Tp pinning transition must not be confused 
with the distinct higher temperature transition, occuring 
just below T/Vyi, in which the putative unpinned (T > Tp) 
SBG "melts" into the NEG by unbinding of dislocations. 
It is most certain that the transition at Tp, derived here 
within the harmonic elastic approximation, will be con- 
verted into a crossover by t he n eglected elastic nonlin- 
earities, analyzed in Sees, Vl- VIII , However, it is possible 
that a rounded "ghost" of the transition at Tp will per- 
sist in our full theory and will therefore be experimentally 
observable. 

It is enlightening to contrast these RG flows, with 
their perturbative fixed line, Eq,3.22, at which the rel- 
evant flows of Ay terminate, with those of disorder- 
free thermal SinfeGordon models (e.g. thermal rough- 
ening transition where the cosine coupling runs away 
to strong coupling at low temperatures. In these lat- 
ter systems the cosine potential becomes relevant upon 
lowering the temperature, with the system settling down 
in one of the minima of this periodic pinning potential. 
As the cosine coupling continues to grow, it further re- 
duces thermal fluctuations, suppressing their ability to 
average away its pinning effects, thereby further increas- 
ing its pinning influence. Mathematically this manifests 
itself in the positive nonlinear contribution to the flow 
equation for the cosine coupling, which results in the ab- 
sence of a stable finite coupling fixed line. In contrast, 
for disordered problems of the type considered here, as 
the disorder becomes relevant at low temperatures and 
begins to grow, it leads to fluctuations (roughening) that 



are larger than those from purely thermal fluctuationso. 
This disorder-enhanced roughening subsequently leads to 
a more effective averaging away of the pinning potential 
thereby suppressing its effects beyond a certain strength. 
Mathematically, this is captured by t he n egative sign of 
the nonlinear Ay term in the flow Eq.3.21, which results 
in the termination of the flow of Ay and the glassy flxed 
line described b y Eg . 3.2S . 

The flow Eq. 3.21 also implies that the random-field 
disorder is irrelevant for 77 > 4. Since the bulk modulus 
B vanishes, while K remains finite through Tna, rj di- 
verges as r — > T^^, and hence, sufficiently close to the 
NA transition temperature T^na, we are guaranteed to 
have a range of temperatures (within the smectic phase) 
over which the random-field disorder is irrelevant. How- 
ever, as we will see below, because tilt disorder Ah is 
a strongly relevant perturbation, the three-dimensional 
quasi-long-range smectic order for rj > A will be converted 
into short-range correlations, even when the random field 
disorder Ay is irrelevant. 

It is essential to stress that the renormalization group 
fiow diagram described above (i.e. relevance for 77 < 4 
and irrelevance for 77 > 4 of the random-field disorder) 
survives even despite the strong relevance and runaway 
of the random tilt coupling A;^. This occurs because the 
dimension of the cos [go(wct ~ Ufj)] operator (the random 
field disorder) is independent of Ah, to all orders in per- 
turbation theory in Ay. This can be easily seen to first 
order in Ay: 



-Hint) 



Ay 

T 



/ d.^i~^{(^oa[qo{ua - Up)]) 

a/3 

^ / dV^cosfeo^-^Dle-^"" 
^Jd\Y.'Os[qo{u<~u<)]e-^'(' 



-5„ 



where 



/a,3=go'[G'L(^ = 0)-G>^(r = 0)] , 
in G'qjj(O) there is no implied sum over a, and 



(3.24) 



(3.25) 



G>,(r) 
Using 



(2^)2 27r 



dqz 



Gaa{q) ~ Gaf}{q) = 



T 



Kql + Bq, 



GM^'^"" ■ (3.26) 



■(l-<5a^), (3.27) 



which is obviously independent of the tilt-disorder cou- 
pling Ah, and Fourier transforming, we obtain the Ah- 
independent 77 given in Eq.3.19a. Performing length 



14 



Eq.3.21 



IS 



rescalings, Eqs.3.17, to restore the new cutoff Ae ^ back 
to A give s, to first order in Ay, the flow equation for Ay, 
Eq. |3.184 

To see that the recursion relation for Ay 
independent of A^ to all orders in Ay, in this otherwise 
harmonic theory, on e can return to the non-replicated 
Hamiltonian, Eq.3.5 and completely eliminate the ran- 
dom tilt field h(r) via a change of variables 

uir) = u'ir) + /(r) , (3.28) 



with /(r) completely determined by the random tilt field 
h(r) via 



(-i?92 + xvl)/(r) = V^.h(r) 



(3.29) 



Inserting Eq. ^.28 



into Eq.3.5, and dropping terms that 
only depend on the random tilt field h(r), since, as dis- 
cussed earlier, these have no effect on any correlation 
J functions, we have 



H = 



B 



K 



{V\u'f + h(r) • + Bdjd.u' + KV\fV\u' - |i^o| (l^'(r)e*«""'("") + v'* (r)e-^*"'("")) 



(3.30) 



where we have defined V'(r) = e*''°-^^''^y(r). Note that 
the statistics of V'{v:) are exactly the same as the statis- 
tics of V{y), since all we have done is added a new ran- 
dom phase factor to V{y). But since the phase of V{y) 
was uniformly distributed between and 27r, the phase 
of V^'(r) will also be. H ence, its statistics are the same. 
Therefore, the model Eq. 3.30 is only affected by the pres- 
ence of the random tilt field h(r) through the terms 



AH 



h(r) • Vxu' + Bd,Jd,u' + K'^\fy\u' 



(3.31) 



which, after integration by parts, vanish by virtue of the 
choice Eq.3.29. Thus, the new Hamiltonian is unaffected 
by the random field h(r). Since the partition function is 
obviously invariant under a change of variables, it and, 
therefore the recursion relations for the parameters B, K, 
and most importantly Ay, cannot be affected by h(r), 
and hence are independent of Ah- 

The results described above for the disordered three- 
dimensional smectic-A phase are quite closely analogous 
to those of Cardy and Ostlund for random symmetry 
breaking fields in the two-dimensional XY-modelEi, ex- 
cept that for pinned smectics, the random-field disor- 
der (the cosine) is guaranteed to become irrelevant as 
T — > T^^ since ri{T Tma) oo. 

As can be seen from the detaile d reno rmalization group 
analysis of Appendix A and Eq. 3.18d , even if the bare 
A/j = (i.e., no initial tilt disorder), this type of disorder 
is generated by the random-field disorder (Ay) after the 
high-wavevector degrees of freedom are integrated out. 
In contrast to the 2d random-field XY-model, where the 
generated A/j disorder is only marginally relevant and 
only weakly affects the QLRO order found for Ah = 
(converting ln(r) phase correlations to In^(r)), for the 3d 
smectic-A phase the A;j tilt disorder is strongly relevant. 
As we will see in the next section, the dimensionless cou- 
pling that determines the effect of the tilt disorder is 



From the recursion relations Eqs j3.18a[ - p.l8d| we find 

(3.33) 



dl 



For ?7 > 4, Ay(^) — > (as we have seen), and so 
dAh/dl = 2A/i, which is trivially solved to give 



Ah(l) = Ahe"' , 



(3.34) 



where A^ is the non-universal dimensionless "bare" cou- 
pling, Eq. 3.3S. Thus, we see that the tilt disorder is 



strongly relevant, in contrast to the behavior of the two 
dimensional random field XY model,El for which the tilt 
disorder was marginal, in the RG sense, in the phase in 
which the cosine was irrelevant. 

For ?7 < 4, Ay ^ A^ > 0, with A^ given by Eq.E23 



Now, in the 2d random field XY model, the existence 
of a non-zero Ay in the low temperature phase implied 
completely different behavior for the tilt coupling Ah{l) 
than in the high temperature phase: a runaway of Ah{t} 
to infinity as ^ — > cxd in the low temperature phase, as 
opposed to^ constant Ah{t — *■ cxd) in the high tempera- 
ture phase.0 In our problem, however, Ah{t) runs away 
to infinity in both phases, and asymptotically in exactly 
the same way\ One can see this by solving Eq. |3.33 with 
Ay (i) replaced by its non-zero fixed point value Ay given 
by Eq.p.23|, finding 



A^(^) = [A, 



2 

/ A * ^2l 21 



£9o 
A2 



* \2 



(3.35) 



A,, = AiAi, 



(3.32) 



which has the same asymptotic behavior as £ — s- oo 
(namely Ahoe"^^) as in the low temperature phase. Non- 
universal constants (like A^o) can be different in the low 
temperature phase, but the scaling (e^^) is not. We will 
see in the next section that this implies that equal time 
correlation functions scale in exactly the same way in 
both the pinned (77 < 4) and non-pinned (77 > 4) phase. 
The only difference between these two phases is in their 
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dynamics, which are divergently slower in the randomly, 
pinned phase, as we shall show in a future publication.E£l 
Physically the strong relevance of A/i is expected, be- 
cause in the smectic-A phase the rotational invariance is 
spontaneously broken. Since each realization of the ran- 
dom tilt disorder explicitly breaks the rotational invari- 
ance (preserving it only statistically), the smectic layer 
orientation has a divergent response to tilt disorder. 
At long scales where Ah w is no longer valid, the 
quasi-long-range order, implied by the nontrivial fixed 
line, Eq,3.23, will undergo a rapid crossover to short- 



range correlations of tp{r), even in the regime (77 > 4) 
where the random-field disorder (which always generates 
random tilt Ah) is irrelevant. 



IV. RANDOM TILT-ONLY HARMONIC ELASTIC 
(TOPOLOGICALLY ORDERED) MODEL 

For finite Ah both the random field (Ay) and ran- 
dom tilt (A/i) disorders must be treated simultaneously. 
Since, as described above, our conclusions about the 
phase diagram for Ay are not affected by the tilt dis- 
order Ah, we can study the effect of A^i on smectic-A 
order within the regime ry > 4, where the random-field 
disorder is irrelevant. In this regime, at long enough 
scales the smectic-A phase is effectively subjected to tilt- 
disorder only. We can therefore safely set Ay — and 
analyze exactly the remaining quadratic theory for any 
value of Ah- 

As argued in the previous section, and to be shown in 
the next, the scaling of the results we find here will also 
apply in the pinned phase where 77 < 4 and the cosine 
disorder is relevant. 

Since, within the harmonic-elastic approximation of 
this section, the effective Hamiltonian is quadratic in the 
displacement fields u, when Ay = all the correlation 
functions can be computed exactly. In particular the 
quantity of interest is 



Fourier transforming the u-fields, we get 
(fq±dq. 



(4.1) 



C(r^,z) = 2 



(27r)3 



-a 



S^q + q') 



We can write the quenched and thermal averaged 
{u{q)u{q')) in terms of their replicated correlation func- 
tion: 



(u(q)u(q')) = (Ma(q)Ma(q')) , 

= 53(q + q')G„„(q) 



(4.3) 



where no sum on a is implied, and the replica propagator 
Gap is given by Eq.3.12. Using that equation, and Eq.4.3 
for the correlation function, gives 



C(r^,z) = CT(rj.,z) + CA(ri,z) 



(4.4) 



where we have separated C{r±, z) into "thermal" Ct and 
"disorder" (frozen) Ca parts, given by 



Ct = 2r 



(p-qA_dq^ 1 - 6*1"" 
(2^)3 Kql + Bql 



T 



2ttVKB 
T 



2ttVKB 



zX 




9 1 

rl 




a / 





4A|z| 



(4.5) 



and 



Ca = 2Ah 



A, 



327r52A3 



d\^dq, qj (1 - e'^ ") 
(2-)3 {Kqi + Bq^^f 

-rl/iX\z\, 



Ei 



{4A| 

-4AL, 
Ll 



z\e 



Ei 



(-ri 
W\z\ 



(4.6) 



where Ei{x) is the Exponential Integral function, and, 
as promise d, the effective dimensionless coupling Ah, de- 
fined in Eq.3.32, naturally appears above. In the above, 



we have considered a finite system whose shape is a rect- 
angular parallelepiped of linear dimensions Lj_ x Lj_ x L^, 
being the length of the system along the ordering (z) 
direction. Unless it has a huge aspect ratio, such that 
Lz ^ L\/\ » L±, any large system {L±,Lz » A) 
will have XLz << L\. In this limit, the asymptotic be- 
haviors of Ca are 



^^[^X\z\+rl\n\if\\ , A|z|»r^ 



Ca 



A. 



XL,^ « L\ , 
A|2;| « r\ , 
XLz « . 



(4.7) 



Note that, although in principle the second {r\ In jL^/zj) 
term in the A|z| >> r\ expansion of Ca dominates the 
X\z\ term in the thermodynamic limit ^ oo (tak- 
ing that limit at fixed and z), in practice, for any 
reasonable system size Lz, the first (A|z|) term actually 
dominates if Ajzj/r^^ gets appreciably bigger than 1. 

An unusual feature of this result is that not only do 
the mean squared fluctuations of m at a given point in 
space diverge as a function of system size, but even the 
relative displacement of two points with finite separations 
(rj^, z) diverge as the system sizes (i_L, Lz) go to infinity. 
This is because the mean squared real space orientational 
fluctuations (|(5n(r)|2) also diverge as L_\_^z oo: 
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(I^n(r) 



= A„ 



(|V^zi(r)|2) , 



\2 ' 



(2^)^ {Kq\ + Bqlf 
Ah r (fq± 



■ In (mm[\/XL^, L±] 



(4.8) 



Although not experimentally relevant, it is instructive 
to generalize this calculation to spatial dimensions d> 3. 
Keeping only the tilt disorder (we will return to justify 
this later), and repeating the calculation just presented, 
we find 



A„|z|(5-'^)/2, A|z| »rl, 



X\z\ << r\, 



(4.9) 



which diverges at large distances for all d < 5. This 
divergence signals the destruction by tilt disorder of 
the (quasi-) long-ranged smectic translational order for 
d < 5. Note that this divergence occurs even for ar- 
bitrarily weak disorder (i.e., arbitrarily small A/j). This 
agrees with the experimental observation that, even when 
the aerogel density becomes very low, smectic transla- 
tional order is still destroyed, as manifested in the non- 
zero width of the X-ray scattering peaks associated with 
the smectic layering. 



How are these results, Eq.4.1, for the equal time cor- 
relation functions modified by the presence of the cosine 
(random field) disorder term in Eqs, 3.4,3.8? Aside from 
a temperature-dependent modification of the prefactor 
A/i — > A/i -f const. (Tp — T)^, below d = 5, they are not 
modified at all. We will no w de monstrate this fact in 
d = 3, and defer to section VIIl the demonstration for 
3 < d < 5. 

For T] > A, the cosine clearly does not affect the long 
distance behavior, since its coefficient Ay is irrelevant 
(i.e., flows to zero under renormalization) . Furthermore, 
we need not worry that it might be dangerously irrele- 
vant, since we have just calculated correlation functions 
with the Ay set equal to zero, and encountered no diffi- 
culties. 

What happens with r] < 47 For i] near 4 (4 — 77 << 1), 
the fixed point value Ay = 4 — 77 of Ay is also << 1, and 
our perturbative (in Ay) renormalization group remains 
valid. We can therefore use this RG to calculate the equal 
time correlations for 77 < 4 and 4 — 77 < < 1 . This calcula- 
tion must itself be done perturbatively in Ay, since the 
full Hamiltonian with Ay ^ is not quadratic. However, 
such a perturbation theory diverges for small wavevectors 
when ?7 < 4. Indeed, the growth under the RG of small 
Ay for 7/ < 4 is a signal and a consequence of this diver- 
gence. Fortunately, we can use the renormalization group 
transformation to relate the correlation functions at small 
wavevector to those at large wavevector, where perturba- 
tion theory should be reliable. Indeed, simply repeating 



for correlation functions the rescalings done earlier for 
the partition function, we can show that: 



C.,(q;i3(0),if(0), A,(0), Ay(0)) ^ ^"jJ^^J'^g^^ 

= e^^+^^'Co^p{e'q^,e'^%;B{£), K {t) , An{l) , Av {I)) . 

(4.10) 

Now, in order to insure that the rescaled correlation func- 
tion on the right hand side of this expression can be safely 
evaluated using perturbation theory in Ay(i'), we will 
choose £ such that e^qi_ = A, the ultraviolet cutoff. With 
this choice, we have 



a/3(q;i?(0),X(0),A„(0),Ay(0)) = 

Ca.f,{A, {K/q^rq,-B{t),K{t),Ah{t), Ay (r)) , 

(4.11) 



qj-J 



where £* = ln(A/g_L), and B = B{0), K = K{0), 
Ah = Ah{0), Ay = Ay(0). Note that T -> 00 as 
q±^0. 



To calculate the right hand side of Eq.4.11 in pertur- 
bation theory in A y (^* ), we expand the cosine in the 
full Hamiltonian Eq.3.8 to quadratic order, obtaining, in 
Fourier space: 



1 /• r " 

H[u^]^- {Bql+Kqi)Y,\u^{ci}\' 
1 " 



a,/3=l 



(4.12) 



which we immediately reco gnize as identical to the tilt- 
only Hamiltonian Eq. 3.1C , except for the replacement 
Ah Ah + 2Avqo/q±- Thus, we can immediately cal- 
culate the correlation function on the right hand side of 
Eq.4.11 by simply making this replacement in the tilt 
only propagator Eq. [3.12 , obtaining 



Ca^A, (A/q^rq,;B{t), K{t),Ah{n,Ay{t)) 

_ TSaf3 

" K{£*)A^ + B{£*){A/q^)^^q^, 
Ah{£*)A^ + 2Avi£*)ql 



[K{£*)A^ + B{£*){A/qj 



(4.13) 



To finish writing this expression entirely in terms of q, 
we need to calculate the Hamiltonian parameters B{£), 
K{£), Ah{£) and Ay(^) from the recursion relations 
Eqs. |3.18a -3.18d and evaluate them at £* = \ii(A/qj_). 
We find 
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Bit) 
K{t) 

Ay(r) 



A, B 



A 



K 



A 

Av(r) 

A_ 



A,, 



' At, 



A^Ai 



(4.14a) 
(4.14b) 
(4.14c) 

(4.14d) 



where Ai = Ai{£ = 0), and in this last equation we have 
assumed that q± is sufficiently small that £* = ln(A/gj^) 
is sufficiently large that Av{i*) will have flown to very 
near its 77-dependent fixed point value Ay, as given by 
Eq. |3.23 . Likewise, in the expression for Ah{£), we have 



used the solution Eq. 3.3f; for A/i(£), the effective tilt cou- 
pling constant, and again assumed that £* — lii{A/q±) 
is ver y lar ge, so that the second (^-independent) term 
in Eq. 3.35| is negligible compared to the first, exponen- 
tially growing term. Both of these approximations be- 
come asymptotically exact as q± — > 0, and £* , as a result, 
00. Inserting these results for the renormalize d ela stic 
and coupling constants into the expression Eq.4.13 for 
the rescaled correlation function, and using that result in 
the matching formula Eq.4.11 for the original correlation 



function at small wavevector yields, after a bit of algebra 



Ca,p{q;B,K,Ah,Ax 



Kqi 



{Kqi + Bql 



(4.15) 



where all of the parameters {K,B,Ah,Ai) in this ex- 
pression are bare parameters (i.e., evaluated at ^ = 0). 
This expression is identical to the result in the Ay ir- 
relevant phase (ry > 4) except for: (i) the Ay((7j^/A)^ 
term, which is clearly negligible, (as q± ^0), relative to 
the {Ah + c q^A'y / Ai)q\ term, and (n) the enhance- 
ment of the strength of the tilt disorder according to 
Ah ^ Ah + cqlA*y / Ai. The heretofore neglected ef- 
fects of the anharmonic terms coming from the cosine on 
the rescaled correlation function (and hence on the origi- 
nal one, since they are related by the matching formula), 
should be even smaller, since A^ is << 1 for 4 — << 1. 
As discussed earlier, this argument cannot be invalidated 
by possible infra-red divergences in the perturbation the- 
ory, since we are calculating the rescaled correlation func- 
tion at large wavevector. The effects of these divergences 
on the short wavelength correlation functio ns are implic- 
itly included in the matching expression Eq. 4.11 through 



our use of the renormalized parameters on the right hand 
side, since these rcnormalizations include anharmonic ef- 
fects. 

Hence, at long wavelengths, aside from the modifica- 
tion of the prefactor Ah Ah -\- const. (Tp — T)^ (for 
T < Tp), we recover the same asymptotic form for the 
equal-time correlation functions in the randomly pinned 
(77 < 4) and non-pinned (77 > 4) phases. Although, 
strictly speaking, we have only derived this result near 
rj — A, where our perturbation theory is valid, it must ap- 
ply throughout the entire low temperature phase, since 
only a phase transition can change the asymptotic scal- 
ing of the correlation functions. Barring the existence 
of such a phase transition controlled by some strong cou- 
pling fixed point, which we of course cannot rule out with 
the perturbative analysis performed here, our results for 
the scaling behavior of the equal-time u — u correlation 
fimctions should persist all the way down to 77 = 0; i.e., 
throughout the smectic phase. 

In any case, in this elastically harmonic phonon (no- 
dislocations) model, there is guaranteed to be a pinned 
phase below rj = A whose static correlation functions 
have r (q) dependence identical, at long wavelengths, 
with those in the non-pinned phase above 77 = 4. The 
distinction between these phases lies in the temperature 
dependence of the parameters in these static correlation 
functions (as we discuss below) and in their dynamical 
properties. 

That is, it is important to note, that in spite of the 
above finding that static correlation functions scale (with 
position or wavevector) identically in both the 77 > 4 
and 77 < 4 phases, the existence of the transition at Tp 
to the low temperature translationally pinned (77 < 4) 
phase can still in principle be detected in a static experi- 
ment. This can be seen by noting that in the presence of 
such a transition, all A/j-dependent physical quantities, 
which would otherwise be smooth functions of temper- 
ature, are non-analytic in T at Tp. That is, because 
of the additional contribution to A^*^ in Eq.4.15, pro- 
portional to (Ay)^ cx 
Af (T) = A;,-Kconst.(Tp- 
ical quantities derived from Ah will have a discontinuous 
second derivative. Although this is a quite subtle effect, 
it is an unambiguous static experimental signature of the 
transition to the pinned smectic glass phase at T = Tp. 
However, as briefl y disc ussed just below Fig. 6, we demon- 
strate in Sees, Vl- VIII that this transition, derived here 



|Tp - Tp, that arises below Tp, 
T)^. As a result, other phys- 



within the harmonic elastic theory, will be rounded once 
nonlinear elastic effects are taken into account. It is, how- 
ever possible, that rounded remnants of the nonanalytic- 
ities at Tp, discussed above, will still be experimentally 
observable. 
The static 



4.15 



correlation function derived in Eq, 
shows that long-ranged smectic translational order is de- 
stroyed in d = 3. To make quantitative comparisons with 
experiments, it is useful to calculate translational corre- 
lation lengths and , whose inverses will give the 
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width of the broadened X-ray diffraction peaks. These 
correlation lengths are the distances r± and z at which 
the mean squared relative displacement correlation func- 
tion C{r±, z) is of order a? , where a — 27r/qo is a lattice 
constant. 

When there is no disorder (A/i = 0), the above crite- 
rion leads to 



(4.16a) 
(4.16b) 



As disjcussed in the extensive literature on the 
subjectjEJ however, the slow (logarithmic) divergence of 
the thermal part Ct of C implies that, for ry < 2, the 
X-ray scattering peaks do not, in fact, become broad 
for |q_L| < 1/^"^; instead, they become power-law diver- 
gences, rather than the Lorentzians one might otherwise 
expect. 

The more strongly divergent part Ca of C does lead 
to genuine broadening, for |q| < l/Cj^^, where ^jf^ are 
defined as those length scales (i.e., values of r± and z) at 
which Ca {r±_ ,z) = a? . Having computed Ca within the 
elastically harmonic theory, we thereby find 



X 



,8ttKB 



Aftln(2VAi;/ef ) 



1/2 



(4.17a) 
(4.17b) 



In fact, we powder average the X-ray scattering, since, 
as we will show in a moment, the smectic in aerogel lacks 
long-ranged orientational order, as well as translational 
order. The broad ring of X-ray scattering that results 
from this averaging will have its width determined en- 
tirely by : 



^powder — (^jg ) 



(4.18) 



a gratifyingiy simple prediction. 

Although it is more difficult to measure experimentally, 
it is nonetheless of interest to calculate the harmonic ori- 
entational correlation lengths S^j^ ^ , which are defined as 
the values of Lj_^z beyond which the mean squared ori- 
entational fluctuations (|^np) get to be of order 1. If we 
ignore dislocations and anharmonic elasticity, we can ob- 
tain these lengths by simply equating the "pure phonon" 
result Eq.4.8 for (|(5np) to 1. This gives: 



= ae4.^^A-VA, ^ ae^it/2a^ ^ (4.19a) 

A A 

In the next section, we will show that, unsurprisingly, 
these lengths give the distance beyond which disloca- 
tions unbind, invalidating the purely elastic phonon the- 
ory studied in this section. 



Thus, orientational order persists out to much larger 
distances than translational order, in the limit of weak 
disorder where all the correlation lengths get large. In- 
deed, in this limit the orientational correlation lengths 
grow exponentially with the translational ones. This 
qualitatively agrees with experimental determinations of 
the orientational correlation length, as indirectly inieiTed 
from specific heat data and related measurements. E3 

In Sec. VI, we will show that anharmonic elastic effects, 
which we have ignored up to now, change the relation be- 
tween the X-ray correlation length and the orienta- 
tional correlation lengths ^ from exponential to power 
law. However, the fact that ^ both remain >> 
continues to hold, validating our use of elastic theory 
(which was predicated on the assumption that orienta- 
tional fluctuations are small), to calculate the X-ray cor- 
relation lengths. In any case, it would be very interesting 
to measure the relation between translational and orien- 
tational correlation lengths as a function of temperature. 

All of the above results apply subject to our two initial 
assumptions: (1) that dislocations were not generated by 
the disorder, and (2) that anharmonic terms in the elastic 
Hamiltonian could be neglected. 

In the next section, we will show that, if we continue 
to assume (2), assumption (1) is wrong: in the harmonic 
elastic approximation, in 3d dislocations are created even 
by arbitrarily weak disorder. 

However, the effects of these dislocations turn out, in 
the weak disorder limit, to be felt only on length scales 
longer than the smaller of the orientational length ^ 
and dislocation unbinding length ^, both of which are 
much longer than the translational correlation lengths. 
Thus, our above calculations of the correlation lengths 
remain valid. However, in the presence of these unbound 
dislocation loops, the smectic glass phase is destroyed. 
It may be replaced by a nematic glass phase, however, 
as discussed in the Introduction. The static properties of 
our system in the presence of dislocations are analytically 
accessible, and we analyze them in the next section. 

Furthermore, the behavior of the correlation functions 
found above is modifled by clastic nonlinearities (which 
lead to anomalous elasticity), which also have a no ntriv - 
ial effect on dislocations, as we demonstrate in SeclVHl. 



V. DISLOCATIONS IN THE RANDOM 
TILT-ONLY HARMONIC MODEL 



In the previous sections, we have shown that there is 
no long-ranged smectic translational order in the pres- 
ence of the random pinning and tilting fields. However, 
this by itself is not sufficient to prove that there is no 
phase transition in this model. Indeed, we know of many 
examples of transitions between two phases which both 
lack long-ranged order. Tba Kosterlitz-Thouless transi- 
tion in the d = 2 XY modeluj is perhaps the most famous 
example. In that problem, the transition is associated not 
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with the disappearance of long-ranged order, but rather, 
with the unbinding of neutral pairs of topological defects 
(vortices) with increasing temperature. 

It is reasonable, therefore, to ask whether the same 
thing can happen in our model: are topological defects 
(i.e., smectic dislocation loops) still bound even in the 
translationally disordered phase we have discussed? If 
they were, then an equilibrium phase transition would 
be required to produce the expected high temperature 
phase, in which dislocations are unbound. 

Of course, the divergences that destroy long-ranged 
order are much stronger (power-law) in our model for 
d < 5 than in the 2d XY model, where they are logarith- 
mic. However, there are examples of phases much more 
strongly disordered than the Kosterlitz-Thouless phase-Ln 
which topological defects nonetheless remain bound.Cj. 
FurthermoEe, there has been considerable speculation 
recentlydild that a "Bragg glass" phase might exist in 
pinned superconducting flux- line lattices. This "Bragg 
glass" would be a phase in which the random pinning 
destroyed the translational order of the flux lattice, but 
did not induce dislocations in the lattice. It seems quite 
reasonable, therefore, to ask whether an analogous phase 
occurs in smectics. This section addresses this question, 
and shows analytically that smectic Bragg glass does not 
occur in d = 3 (within a model with harmonic elasticity) . 
We will show later that this result maybe invalidated by 
elastic anharmonic effects. 

The starting point of our analytic theory is the "tilt 
only" model, by which we mean Eq.3.5 with the ran- 
dom potential V{r) set to zero. As discussed earlier, 
this theory correctly reproduces all of the static corre- 
lation functions in both the pinned (T < Tp) and the 
non-pinned (T > Tp) regimes. However, in view of the 
very strong irrelevance of Ay to the static correlation 
function in both phases in the "phonon only" approxi- 
mation, it seems quite plausible that it is also irrelevant 
when dislocations are included. Indeed we will prove this 
in Sec.fxl. 



Setting V{r) =0 in the Hamiltonian Eq.3.5 reduces it 
to a quadratic theory 



H[u] 



B 



K 



(Viw)2 + h(r) • V^u 
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Fig. 7. An edge dislocation in the smectic liquid crystal, 
with the core coming out of the page at the position 
indicated by the lazy "T" symbol. 



This is the edge of an extra smectic layer inserted into 
the smectic. Clearly, a pure edge dislocation must form 
a closed loop lying in the smectic plane, as illustrated 
in Fig. 8. More generally, dislocations can tip out of the 
plane of the smectic layers. An extreme case is a screw 
dislocation, which runs perpendicular to the layers. 




(5.1) 



The effect of dislocations on such a theory can thetii be. 
treated using the techniques that many authorsac3c3 
have applied to a variety of disorder free systems, includ- 
ing smectics. 

Most aspects of this procedure are already described 
quite well in the literature; therefore, our explanation will 
be brief until we come to those points that are affected 
by the presence of the disorder. 

We begin by recalling the definition of dislocations in 
smectics.Ell The simplest type is an edge dislocation, il- 
lustrated in Fig. 7. 



Fig. 8. A three dimensional illustration of dislocation loop 
in the smectic liquid crystal. 

Mathematically a dislocation is a line, or, more gener- 
ally, a curve, with the property that when the gradient of 
the displacement u is integrated around any curve that 
encloses the dislocation, the result is not zero, as it would 
be in a dislocation free system, but rather an integral 
multiple N of the layer spacing a. Mathematically, this 
means 
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Vu- dl = aN , (5.2) 

or, in differential form: 

V X V = m , (5.3) 
where we have defined 

V = Vu , (5.4) 

and 

m(r) = XI / o-Niti{s^)5'^{Y - Vi{si))ds., , (5.5) 

i 

where Si parameterizes the i'th dislocation loop, ri(si) is 
the position of that loop, ti{si) is its local unit tangent, 
and Ni the "charge" or number of excess layers associ- 
ated with the dislocation. Note that Ni is independent 
of Si, since the charge of a give n lin e is constant along 
the line defect. Furthermore, Eq j5.3| implies that 



V • m(r) = , 



(5.6) 



which simply means that dislocation lines cannot end in 
the bulk of the sample; they must either form closed loops 
or extend entirely through the system. 

Now, our procedure for adding dislocation lines to our 
previous pure phonon model Eq. 
dard one. We separate the field v 



is the following stan- 
Vm into phonon 



(single-valued) and dislocation (singular) parts 



(5.7) 



where the dislo cati on part minimizesca the ela stic 
Hamiltonian Eq. ^.lj , subject to the constraint of Eq.5.3 
(V X V = m). (For an alternative derivation see Ap- 
pendix C.) This uniquely determines Vd(r) given the dis- 
location configuration m(r). We then insert the decom- 
position Eq.pJ back into the elastic Hamiltonian Eq. 5.1 



As a r esult of the construction that Vd(r) minimizes 
Eq.^.ll, all the cross couplings between Vd(r) and Vp(r) 
vanish. We are thus left with a decoupled elastic Hamil- 
tonian for m(r), which we can use as a basis for studying 
the statistical mechanics of dislocations. 

Let us now implement this procedure. T he E uler- 
Lagrange equation, obtained by minimizing Eq.BJ^, is 



(B9f - KVi)ud{r) + Vi_ ■ h(r) - . (5.8) 
Rewriting this in terms of Vd(r) = VM(;(r) gives 



9,«^-A2viV^.vi + -V^.h(r)=0, (5.9) 

where = K/B. 

In Fourier space, this becomes 

q^v'a + X'qlq^ ■ + iq^ • h(q) = . (5.10) 



The constraint Eq j5.3| becomes, in Fourier space, 

iq X v<j = m , (5.11) 
which has the general solution 
iqx m 



q0 



(5.12) 



where (/) is the smooth elastic distortion around the dis- 
location line , to be determined by the Euler-Lagrange 
equation Eq. 5.5. Inserting the above expression for v,; 
into Eq. 5.10| and solving for (j) gives 



iqzjl - \^q\)ezijqimj _ q± h 



(5.13) 



where we have defined the inverse of the smectic propa- 
gator 



(5.14) 



Inserting this into Eq.5.12 and then substituting this 
final expression for into the original elastic Hamilto- 
nian Eq.5.1 gives the defect interaction Hamiltonian 



-^P,j-TOi(q)TOj- 



-q) + m(q) • a(-q) 



(5.15) 



where Pjj"(q) = Sj^ — q^qj' /q\, a(q) (not to be confused 
with the lattice spacing a) is a Fourier transform of the 
quenched field related to the original random tilt field 
h(q) via 



a(q) 



q X h (z X q) q • h 

q^ r„g2 



-7.(1 



(5.16) 



and we have dropped unimportant terms that depend 
only on the quenched random variables (and not on m 



The first (h-independent) term in this expression Eq. 5.15 
is just the usual smectic dislocation energy for a pure 
system.c3 

To treat this model, we perform a duality transfor- 
mation. We begin by putting the model on a simple 
cubic lattice (to make the model well defined at short 
distances); now, m(r) is defined on the sites r of the 
lattice, and takes on values 



m(r) 



d? 



(n:,(r),7iy(r),n2(r)) , 



(5.17) 



where the n^'s are integers, and d is the cubic lattice con- 
stant used in the discretization. The partition function 
for this model is then 



Z[h(q)] 



{m(r)} 



Sfml 



(5.18) 



where 
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i/,[m]+i?,-^^|m(r)| = 



(5.19) 



and the sum is over all discrete configurations of m's 
given by Eq.5.17, satisfying the dislocation line continu- 
ity constraint 



V • m = 



(5.20) 



where the divergence now represents a lattice divergence, 
Hd is given by Eq. ^.lf: , and we have added a core energy 
term Ec |m(r)|'^, to account for energies near the core 
of the defect line that are not accurately treated by our 
continuum elastic theory. We call the reader's attention 
to the fact that the partition function still depends im- 
plicitly on the configuration of the random tilt-disorder 
fields {h} through a in Eq.5.15. 

To proceed, we enforce the constraint V • m = by 
introducing a new auxiliary field 0(r), rewriting the par- 
tition function Eq. [5.18| as 



Z = \{jde{v) g-S[m]+.^^fl(r)V.m(r)dVa^ 

(5.21) 



{m(r)} 



where the sum over {m} is now unconstrained. The con- 
straint V • m(r) = is enforced by integration over 0{y), 



(5(V •m(r)) = 



dO{r) .6l(r)V-m(r)dVa 

27r 



(5.22) 



where the (5 is a Kronecker delta, since mrf^/a, and, 
hence, V • md^/a, are integer-valued. 

Now we can "integrate" (actually sum) by parts, and 
rewrite 

6'(r) V • m(r) = — m(r) • V(?(r) + surface terms . 

r r 

(5.23) 

Our next step is to introduce a dummy gauge field A 
to mediate the long-ranged interaction between defect 
loops in the Hamiltonian Eq.5.15. This is accomplished 



by rewriting the partition function as 



Z = \{ j de{v)dK{v) e-^[™'^^^l(5(V-A)^(A,) , 

(5.24) 



{m(r)} 



J with 



Td- 



m(r) • [-i V6l(r) +d3(^A(r) +a(r))] + E^ — \ni{- 



2T ^ Kql ' 



(5.25) 



where a is the quenched gauge field defined in Eq. 5.16| . 

It is straightforward to check that, upon performing 
the Gaussian integral over A, subject to the indicated 
constraints V • A = = 0, we recover the original long- 



ranged interaction between dislocation lines in Eq. ^.l5 



The two goals of all of these manipulations have now ^ 
Z = Y[J deir)dA{r)6{V ■ A{r))5{A,) exp - ^ Vp[e{r 



been achieved: the sum on {m(r)} is now unconstrained, 
and the sum on each site over m(r) is now decoupled 
from that on every other site. Furthermore, this sum is 
readily recognized to be nothing mor(j_than the "peri- 
odic Gaussi an" m ade famous by VillaiiO. The partition 
function Eq.5.24 can thus be rewritten 



E 



Kqi' 



(5.26) 



where the well-known 27r-periodic Villain potential Vp{x), 
defined by 



Ea/T+ixn 



(5.27) 



n— — oo 



has the usual property that the smaller Ec/T is (i.e., 
the higher the temperature in the original random-tilt 
smectic model), the sharper the potential minima. Thus 
raising the temperature in the original model is lik e low- 
ering the temperature in the dual model Eq.5.26. It is 



precisely this familiar temperature inversion associated 
with duality that leads to an inverted XY transition for 
three-dimensipaal disorder-free superconductorst£l and 
bulk smecticsliB. It also plays an important role here, 
as we shall see in a moment. 

Standard universality arguments imply that replacing 
the periodic potential Vp{x) in Eqj5.26|by any other non- 
singular periodic function should not change the univer- 
sality class of the transition. In particular, we could re- 
place Vp{x) by cos(a;). The resultant model would be 
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precisely the "fixed length" version of the "soft spin" , or "action' 
Landau-Ginsburg- Wilson model, with the complex 



+ E^|A(q)P, (5.28) 



ad 



(iA + a)]V* 



ad 



(iA + a)]V + t|V'l' + u|V'l' 



where "0(1") is a complex "disorder" parameter field whose 
phase is 9{v); the reduced temperature t, quartic coupling 
u, and c(T) are parameters of the model with 



c(T) = 



<PVp{x) 



dh 



X 0(1) , 

x=0 

\ T/2E„ T » 



(5.29a) 
(5.29b) 



Because of the duality transformation's inversion of the 
temperature axis, the reduced temperature t is a mono- 
tonically decreasing function of the temperature T (of 
the original dislocation loop model), which vanishes at 
the mean- field transition temperature Tmf of the fixed 
length model Eq. ^.26 . 



Universality also implies that this "soft-spin" model 
should be in the same universality class as the fixed 
length model Eq. 5.26 . We shall, therefore, henceforth 
work with model Eq. 5.28 , because it is more straightfor- 
ward to analyze perturbatively. 

As we undertake that analysis, it is important to keep 
in mind that, as a consequence of the duality inversion 
of the te mpera ture axis, the ordered phase of the dual 
model Eq. ^.28| corresponds to the disordered (i.e., dislo- 
cation loops unbound) phase of the original dislocation 
loop gas model. That is, the low dual-temperature phase 
described by 



(V'(r)) ^ , 



(5.30) 



corresponds to the disordered dislocation-unbound phase 
of the smectic liquid crystal. 

In the absence of disorder (a — 0) this model is ex- 
actly the dual version of the NA [dislocation loop model 
derived and studied by one of us .113 



Disorder is included in Eq. 5.28 through the quenched 
gauge-field a(r), which is related to the random tilt field 



h(r) by Eq.5.16. The partition function 



Z[h\ = j [d^P][dA]e-^^'l''^'^^{V ■ A)5{A,) , (5.31) 



with S given by Eq. 5.28 , is thus an implicit function of 
the random tilt field configuration h(r). 

As in the previous section, we will cope with this de- 
pendence of Z on the quenched field h(r) using the replica 
trick. Doing so leads us to calculate 



[da] n MA„] 



-Sp[i/J„,A„,a] 



xP[a]<5(V-A„)<5(A^) , 



J with 



St- 



E 



E 



r„ 



2TKqi 



|A„(q) 



(5.32) 



(5.33) 



The p robability distribution P[a] of the field a in 
Eq. ^.32 is Gaussian, since a is linear in h and the distri- 
bution of h is Gaussian, defined by Eq.2.22. Thus, the 
distribution P[a] is completely specified by the average 
Qi( q)a-,-( — q). This is easily evaluated, using the relation 
Eq. |3.16 between a(q) and h(q). We find 



where eijk is the usual fully antisymmetric third rank 
unit tensor. To obtain the second equality we have used 
/ii(q)/i„(— q) = AhSin, and defined the standard trans- 
verse projection operators 



aj(q)aj(-q) = /i;(q)/i„(~q) 



qiq^il^X^ql) 



(5.34a) 



Pi J = Si J - q^qj/q , 

/^| = (l-<5..)(l-<5,.)(<5. 



(5.35a) 
(5.35b) 



We now need merely to consid er the sta tistical me- 
chanics of the model defined by Eqs j5.32| , |5.33 , in the limit 

From previous classic worksli3 we have consider- 
able experience analyzing the critical properties of gauge 



theories of the type defined by the action in Eq. 5.33 



As we know from work |©ii the zero field Normal-to- 
Superconductor transitionlla, and the analogous analysis 
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of the bulk, clean Nematic-to-Smectic-A transitionJl^, 
there are two possible regimes: (i) the extreme type I 
regime, in which u « 1, and, as a result, the tp fluc- 
tuations are subdominant to those of the gauge field A, 
allowing a mean- field treatment oi ip. In this case the re- 
maining gauge field fluctuations can be treated exactly. 
(ii) the type II regime, in which m > 1 and both ip and A 
fluctuations must be treated within a considerably more 
involved RG analysis. 

A complete analysis of the critical properties of this 
smectic dislocation loop unbinding transition, described 
by the actiour-ip Eq.5.33, is beyond the scope of the 
present paper .Ea Here our goal is to determine the stabil- 
ity of the topologically ordered smectic Bragg glass phase 
to disorder-induced dislocation loop unbinding. To an- 
swer this question, it is sufficient to study the effect of 
diagrammatic corrections on the reduced dual tempera- 
ture t, as we demonstrate in detail in Appendix D for the 
type I limit, delasdng the analysis of type II limit to a fu- 
ture publication.Ej By computing the disorder-averaged 
free energy, we find that in either regime the lowest or- 
der contribution to the renormalized dual temperature t 
comes from the average of the "diamagnetic" terms 



r 

graphically illustrated in Fig. 9. 



(5.36) 
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Fig. 9. Feynman diagrams that dominate the renormal- 
ization of the reduced dual temperature t, leading to an 
upward shift in the dual Tc- In d < 3 these drive the 
transition temperature in the original smectic model to 
0. 

Generalizing to d dimensions these give 



{d - 2)ca^d'^ 



tR — ^0 



d'^q 
(2^ 



(5.37) 



where the first term in the square brackets comes from 
the first graph in Fig. 9, with the internal wiggly line rep- 
resenting Aq, fluctuations, while the second comes from 
the second graph with the internal dotted line represent- 
ing the quenched gauge field a. The negative sign of the 
disorder contribution leads to an increase in the dual T^, 
an d can be traced back to the fact that the action Sr 
Eq. ^.33 , is complex. 



The second, disorder, term in this integral dominates 
the first as q ^ 0. Indeed, this integral diverges in the 
infra-red for d < 3 



d'^q 



9 9 



q^q^ + X^qiy 



oc 



d q±--^ 



9i 



(5.38) 



where the first proportionality follows from the fact that 
the dominant regime of the integral over qz is qz ^ Xqj_, 
as can be easily seen by contour integration. This diver- 
gence implies that the renormalized dual temperature tji 
is driven to — oo (note the sign) by the disorder in d < 3. 
Indeed, we find in d = 3 



tR — to — 



ca'^d'^Ah 
AT2 



ln(L/a) X 0(1) 



(5.39) 



where L is an infra-red cutoff (e.g., the lateral extent of 
the smectic layers) and a is the ultra-violet cutoff (e.g., 
the size of the liquid crystal molecules ~ lOA). 

The unusu al negative (d = 3) divergent correction to 
t in Eq. |S.39 implies that the dual order parameter is al- 
ways in its ordered phase, which, in turn, implies that 
the dislocation loops of the original smectic model are 
always unbound, thereby destroying the topological or- 
der of the smectic Bragg glass, even at T = 0, for any 
amount of disorder, no matter how small. This conclu- 
sion, of course, only holds within the harmonic elastic 
approximation made in this section, with the treatment 
that approximate ly i ncludes these nonlinear elastic ef- 
fects given in Sec. VII. 

It is important to note that although the dual dis- 
ordered gauge model, Eq. 5.28 , can be extended to any 
dimension, it rigorously only describes the original dis- 
ordered smectic in d = 3, the dimension in which the 
duality correspondence was established between the two 
models. However, if one can entertain the idea that the 
relation between the two models extends to d > 3, then 
we have found a quite nontrivial result for this higher 
dimensional disordered smectic. For 3 < d < 5 the disor- 
dered smectic will then have strongly divergent displace- 
ment fluctuations Eq.4.9, but will maintain topological 
order, i.e. no condensation of dislocation loops takes 
place in this range of dimension (for d > 3). From this 
point of view this higher dimesskmal disordered smectic 
is in fact a true "Bragg" glass Bilil 



VI. ANOMALOUS ELASTICITY IN THE 
TILT-ONLY MODEL OF DISORDERED 
SMECTICS 

Up to now we have analyzed disordered smectics ig- 
noring elastic anharmonic effects. These effects are im- 
portant even in pure smectics, as was disccpered some 
time ago by Grinstein asd Pelcovits (GP)cil. As re- 
ported in a recent Letterll3, we have discovered that in 
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disordered smectics these elastic anharmonicities play an 
even more important role; specifically, we find disorder- 
driven power law divergences, in contrast to the weak, 
thermally-drivea logarithmic divergences that occur in 
pure smecticsBJ Here we provide details of these cal- 
culations, describe their effect on the analysis of previ- 
ous sections r-ati4r4f^'^'^^^ them in the context of recent 
experimentsOEj Ej j—. 

The anomalous elasticity of pure (bulk) smecticsEa 
is characterized by layer compressional and tilt moduli 
B(k) and K(k.) which vanish and diverge, respectively, at 
long wavelengths (k — > 0). This beautiful result is actu- 
ally a general property of all of one-dimensional crystals, 
in which the direction of the Id ordering wavevector is 
chosen spontaneously. As a consequence of this sponta- 
neous breaking of rotational symmetry (a property pos- 
sessed by smectics but not, e.g., charge density waves), in 
the presence of fluctuations a compression can be relieved 
by smoothing out these fluctuations, thereby leading to 
an effective compressional modulus i?(k) that vanishes 
at long wavelengths. This is very similar to the ability 
of free dislocations in a solid to relieve an imposed shear. 
Similarly, in the presence of fluctuations, a bending of 
smectic layers necessarily leads to a compression, which 
implies an effective tilt modulus -R'(k) which diverges at 
long wavelengths. This phenomenon is also analogous 
to the thermally-driven anomalous elasticity of polymer- 
ized membranes, in which the bending rigidity modulus 
diverges, ancLihe-shear and bulk moduli vanish at long 
length scales. E3til 

Although quite fascinating, the anomalous elasticity of 
GP for pure 3d smectics is somewhat academic, because 
for d = 3 the nonlinearities responsible for the anomalous 
elasticity are marginally irrelevant, and therefore lead 
to weakly (logarithmically) fc-dependent elastic moduli. 
These, in practice, are quite difficult to detect experi- 
mentally, and have never been observed in bulk (pure) 
smectics. 

As mentioned above, the main ingredients necessary 
for the existence of anomalous elasticity are sponta- 
neously broken rotational invariance, and the presence 
of fluctuations-i.e. wrinkles in the smectic layers-both 
of which exist (even at zero temperature) in smectics 
in random environments. In this section, we demon- 
strate the existence of anomalous elasticity that is signif- 
icantly stronger, in all spatial dimensions d < 5, than the 
marginal anomalous elasticity of thermal smectics. This 
new, stronger anomalous elasticity is driven by quenched 
disorder, and controlled by a new, zero-temperature fixed 
point that is perturbatively accessible in d = 5 — e. One 
can physically appreciate why this elastic anomaly is so 
much stronger in smectics with quenched disorder by 
realizing that the quenched disordered fluctuations are 
more divergent (at long wavelengths) than their thermal 
counterparts (see e.g., Eq.4.5 and Eql 



As discussed at the end of Sec. [II, in dimensions d < 5 



for the properties of smectics with quenched disorder. 
We therefore expect the effective elastic Hamiltonian in 
Eq.3.4 with U{r) = to be a good starting point for the 
description of smectics confined inside low-density aero- 
gels. However, such a model misses an important ingre- 
dient: nonlinear elasticity which takes into account tha 
underlying rotational invariance of the smectic phase,Eil 
hidden by the spontaneous choice of the layers to stack 
along the z-direction. A careful analysis, starting with 
the de Gennes model,— that keeps track of such elastic 
nonlinearities, leads tocJ 



H[u] 



+ h(r) • V_ 



(6.1) 



as the proper rotationally invariant elastic Hamiltonian 
for the disordered smectic. This form guarantees that a 
uniform rotation of smectic layers costs zero energy in 
the absence of the random tilt field.E3 The field h(r) = 
ff^(r)g(r) is the quenched random tilt disorder, defined 
in Sec|l|. As discussed in that section, it has a vanishing 
average and Gaussian statistics, completely characterized 
by the two-point correlation function 



(6.2) 



which we take to be short-ranged for the reasons dis- 
cussed in detail in Sec.|l[ 

To compute self-averaging quantities, (e.g. the disor- 
der averaged free energy) and to assess the importance 
of the nonlinearities at long scales it is coipenient (but 
not necessary) to employ the replica "trick"Ej that relies 
on the identity InZ = lim„^o(^" ~ I)/"- discussed 
in Sec. [II, this allows us to work with a translationally 
invariant field theory at the expense of introducing n 
replica fields (with the n ^ limit to be taken at the 
end of the calculation). After replicati ng a nd integrating 
over the disorder h(r), and using Eq.6.2, we obtain an 
effective translationally invariant replicated Hamiltonian 



1 f 

•'^ Q = l 



(6.3) 



Q,/3=l 



from which the noninteracting propagator GQ^(q) 
^^^('"a(q)u/3(— q))o (where V is the system volume) can 
be easily obtained (see Sec.pl 



(6.4) 



= TG(q)^„^ + A^giG(q)2 



we expect that the random field disorder to be sig- 
nificantly less important than the random tilt disorder 



with G(q) = l/(ifgi-f Sg2). 

We first attempt to assess the effects of the anhar- 
monicities, disorder and thermal fluctuations by perform- 
ing a simple perturbation expansion in the nonlinearities 
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of H[ua]- The lowest order correction SB to the bare 
elastic compressional modulus B comes from a part of 
the diagram in Fig. 10. A standard analysis gives 



SB = p [TG(q)2 + 2A„giG(q)3] qi , 



(6.5a) 



-B^/^h 



^1 



_^ 2^ J (2^)-^-! {Kqi + BqlY 
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l6 5^^^Mif^ 
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L 



5-d 



(6.5b) 



(6.5c) 



to restore the ultraviolet cutoff back to A. The un- 
derlying rotational invariance insures that the graphical 
corrections preserve the rotationally invariant operator 
[dzU — ^(Vj^w)^), renormalizing it as a whole. It is 
therefore convenient (but not necessary) to choose the 
dimensional rescaling that also preserves this operator. 
It is easy to see that this choice leads to 



2-w 



(6.9) 



where we have dropped the subdominant thermal part, 
focused on d < 5, which allows us to drop the uv-cutoff 
(A) dependent part which vanishes for A ^ oo, and 
cutoff the divergent contribution of the long wavelength 
modes via the infra-red cutoff restriction q± > 1/L, 
where L is the linear extent of the system. The con- 
stant Cd = 27r''/2/[(27r)''r(d/2)] is the surface area of 
a d-dimensional sphere divided by {2n)'^. Clearly the 
anhamonicities become important when the fluctuation 
corrections to the elastic constants (e.g., SB above) be- 
come comparable to the bare values. The divergence of 
this correction as L —^ oo signals the breakdown of con- 
ventional harmonic elastic theory on length scales longer 
than a crossover scale which we define as the value 
of L at which \SB{£_^'")\ = B. In d dimensions, this 
definition gives: 



fl6{5-d)K^/ 



2\ 1/(5-^) 



which for physical 3d smectics is given by 



647r K^/^ 



3 B^^Ah 



1/2 



for d = 3 



(6.6) 



(6.7) 



We can also obtain a non-linear crossover length in the 
z-direction 



for d = 3 



(6.8a) 
(6.8b) 



where A = {K/B f^, by imposing the infra-red cutoff in 
the ^-direction. 

To understand the physics beyond this crossover scale, 
- i.e., to ma ke sen se of the apparent infra-red divergences 
found in Eq. 3.5c - we turn to the renormalization group. 
As discussed in Sec. ID, we employ the standard rriomen- 
tum shell renormalization group transformatiorcd, by 
writing the displacement field as Ua{v) — u<(r) -f u>(r), 
integrating perturbatively in the nonlinear coupling B 
the high wavevector part (r) with support in Fourier 
space in an infinitesimal cylindrical shell (Ae^^ < q^^ < K 
and —CXI < q^ < oo), and rescaling the lengths and 
long wavelength part of the fields as in Eqs.3.17 with 
r± = r^e^, z = z'e"^ and M^(r) — e"^^u^(r'), so as 



This rescaling then leads to the zeroth order RG flow of 
the effective couplings 



K{£) = Ke'^'^-'^"^^^ 



T 



T 



(6.10a) 
(6.10b) 

(6.10c) 



From these dimensional couplings one can construct two 
other effective dimensional couplings. 



1/2 



52 ^ A. ( 



1/2 



(6.11a) 
(6.11b) 



whose flows are independent of the arbitrary rescaling 
exponents uj and 0, and are given by 



5l(^)=5le(3-'^)^ 



(6.12a) 
(6.12b) 



gi is just the coupling that becomes relevant in d < 3 
and was discovered by Grinstein and Pelcovits to lead to 
anomalous elasticity in pure thermal smectics. It is, how- 
ever, only marginally irrelevant in d = 3 and therefore 
only leads, to weak anomalous elasticity in physical 3d 
smectics.cJ In contrast, the upper critical dimension due 
below which g2 becomes relevant is d^c — 5, and we there- 
fore expect a significantly stronger anomalous elasticity, 
that should be experimentally observable in disordered 
3d smectics. These observations imply that temperature 
is a strongly irrelevant variable near the disorder domi- 
nated fixed point and does not feed back in a dangerous 
way. We will therefore set T in all subsequent calcu- 
lations. 

We now turn to the RG computation of the one-loop 
graphical corrections to the flow of the couplings defined 
above. The required integration over the high wavevector 
components of Ua can only be accomplished perturba- 
tively in nonlinearities of H[u\. Since the most relevant 
coupling g2 becomes important for d < 5, we will control 
the infrared divergences by performing an expansion in 
e = 5 - d. 

The change in the Hamiltonian due to integrating out 
these short length modes is 
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6H[u<] - {H,[t 



-^{Hl[u<+ul])\ 



(6.13) 



where Hi[u'^ + u^ ] is the nonlinear part of the Hamil- 
tonian H, Eq. p.Sl , which contains three- and four-point 
vertices. The averages ab ove are performed with the 
quadratic part of H, Eq j3.8| , using the propagator 
Ga/3(q), Eq.B.4, with only an infinitesimal cylindrical 
shell (i.e. no cutoff on q^) of modes integrated out. 
The superscript c denotes a cumulant average. It is easy 
to verify that the first term in SH (first order in Hi) does 
not lead to corrections of the elastic constants. Aside 
from correcting the free energy, it generates an opera- 
tor linear in ((?zU — ^(Vj^it)^) which corresponds to a 
renormalization of the smectic wavevector qq. These cor- 
rections turn out to be finite (irrelevant) near d = 5, and 
we will therefore not keep track of them here. 

The renormalizat ion o f K, B and comes from the 
second term in Eq.6.13, for example from parts which 



are second order in the three-point vertex, with 4 of the 
6 fields contracted. The generic diagram that corrects 
the elastic moduli and the disorder variance is illustrated 
in Fig. 10, with the part diagonal in the replica indices 
a,/? (i.e. part proportional to Sap) renormalizing K and 
B, and the part independent of a, /3 correcting A/i. 




Fig. 10. Feynman graph that renormalizes the elastic 
moduli B and the tilt-disorder variance A^. 

For the calculation of 6B^ the loop integrals can be 
performed at vanishing external momentum, and after a 
short calculation one finds 



SB 



B 



1 r> 



dqz 



271 J {2nY-i {Kq\ + Bqlf 
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(6.14) 



where going from the first line to second line above we 
dropped the irrelevant finite temperature part, and de- 
fined a dimensionless coupling constant 



92 ^ A, ( — j C.-iA-^-^ . 



(6.15) 



Similar, but more involved calculations for K and A/j 
give 



5K « ^g2Kdl 
A, 



A^ 

T 
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T 



dl 



which lead to the following RG flow equations 



dB{£) 
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dK{t) 

dl 
dA^il) 
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92{£))B{£) , 



1 
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g2{i))K{i) , 



(6.16a) 
(6.16b) 

(6.17a) 
(6.17b) 
(6.17c) 



These, together, lead to t he RG flow of the dimensionless 
coupling g2 defined in Eg. 6. 15 



dff2(^) 5 2 

^ - ^92ie) - ^g2ii) 



(6.18) 



where we remind the reader that in this section e — b — d. 
As required, the flow of g2 is independent of the arbitrary 
choice of the anisotropy rescaling exponent lo. The RG 
flow Eq. 3.18 shows that the Gaussian 52 = fixed point 
becomes unstable for d < 5, and the low temperature 
phase is controlled by a stable, nontrivial, glassy T = 
fixed point at 
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(6.19) 



The existence of this nontrivial fixed point leads to 
anomalous elasticity. To see this it is convenient to 
use our RG results to evaluate the connected disor- 
dered averaged 2-point M(k) correlation function G'(k) (x 
(u(k))(u(— k)). As discussed (and utilized) 



(|z.(k)P 

in Sec, III, the power of the renormalization group is that 



it establishes a connection between a correlation function 
at a small wavevector (which is impossible to calculate in 
perturbation theory due to the infra-red divergences) to 
the same correlation function at large wavevectors, which 
can be easily calculated in a controlled perturbation the- 
ory. This relation for G(k) is 

G(ki,A:„if,i?,52) = e(3+'i-)^x 

G(k^e^fc,e"^if(^),i3(^),52(^)), 

(6.20) 

where the prefactor on the right hand side comes from 
the dimensional rescaling (remembering the momentum 
conserving ^- function), after using the exact rotational 
Ward identity (j) = 2 — lo, and we have "traded-in" the 
disorder variable A^ for the dimensionless coupling 172. 
To establish the anomalous behavior of K , we look at 
fcz = 0. We then choose the rescaling variable £* such 
that 



k^e^ =A 



(6.21) 
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We also choose k± sufBciently small such that 52 (^*) has 
reached our nontrivial fixed point gj. Eliminating £* in 
favor of k±, we then obtain 



3-\-d—LU 



(6.22) 

Since the right hand side is evaluated with the transverse 
wavevector at the Brillouin zone boundary, it is easily 
evaluated perturbatively for a small fixed point coupling 
§2- To lowest order we obtain 



(A/fci) 



1 



, (6.23a) 



(6.23b) 



where we integrated Eq. |6.17b| to obtain K{i*), and de- 
fined the anomalous tilt modulus which diverges at long 
length scales 



(6.25a) 
(6.25b) 
(6.25c) 



K{k^) = K{k^/A)-"'' , 
with an anomalous exponent 

1 

= r' 

— - , for d = 3 . 
5 

Similar calculations for the other coupling constants 
and other directions of k show that, in general, 

K{k) = K (fc^^i'^)""" fK{k.^^'^/{k^^^'')<) , (6.26a) 
Bik) = B (fc^e^^^)"" /B(A:.ef V(^±e^''')'^) , (6.26b) 

A,(k) = A, (fc^^i'^)""" fAik.^'^/ik^^^'-'f) , (6.26c) 

with the anisotropy exponent 

C = 2-(ryB+r;K)/2, (6.27) 

which would = 2 in the absence of anharmonic effects, 
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— — , for d = 3 
5 



and 
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1 

10^' 



— - , for d = 3 
5 



(6.28a) 
(6.28b) 
(6.28c) 

(6.29a) 
(6.29b) 
(6.29c) 



Of course, we do not completely trust the extrapola- 
tion of these small e results down to e = 2 (d = 3). 
However, since by definition dg2/d£ = at the nontriv- 
ial fixed point, this condition implies an exact relation 
between the anomalous exponents 



(6.30) 



whi ch is obviously sati sfied by the anomalous exponents. 



Eqs. |6.25b| , |6.28b| , |6.29b 

e 



computed here to first order in 
This Ward identity between the anomalous exponents 
can be equally easily obtained from a self-consistent in- 
tegral equation for the u — u correlations functions, using 
renormalized wavevector dependent elastic moduli and 
disorder variance. l3 

At length scales beyond and the elasticity 

and fiuctuations of the disordered smectic are controlled 
by our new glassy fixed point. One of the important con- 
sequences can be seen in the layer fluctuations, which 
can be inferred from X-ray scattering experiments. For 
instance, layer displacement fluctuations along z are de- 



(6.24) scribed by 



C(z) 



((u(0±,z)-u(0±,0))") , 

d'^k 2(1 -cos(/c^z))A/,(k)fc 



(271-)'* [if(k)fc4 +s(k)fc2]' 



(6.31) 



One can then naturally define the X-ray translational 
correlation length as the length along z at which 
C{z = = a^, where a is t he smectic la yer spacing. A 
simple calculation, using Eqs.6.26a-6.26c leads in 3d to 



where 



= (a/A)^/'^e^ 



X= ivB +Vk)/'2. , 



(6.32a) 
(6.32b) 



(6.33) 



and this result should be contrasted with the expression 
for the X-ray correlation length Eq. 4.17a, calculated 
within the harmonic theory, i.e., ignoring the anomalous 
elasticity. Note that this X-ray correlation length is fi- 
nite even as T ^ 0, as illustrated in Fig. 3. Tius-i-esult 
is consistent with the experimental observationEj'CJ that 
the X-ray correlation length for smectics in aerogel satu- 
rates at some finite value at low temperatures. Note also 
that this length should be different for different smectics 
in the same aerogel, since B, K, and A^, will change from 
smectic to smectic. Since, as discussed in Sec.|l| we ex- 
pect Ah to be a monotonically increasing function of the 
aerogel densi ty pA , with a simplest aerogel model giving 
Ah oc PA, Eq.2.8b, the aerogel density dependence of 
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could test the prediction of Eq. 6.32| . Likewise, the tem- 
perature dependence of could be used to determine 
the ratio C/x> since the bulk K{T) and B{T) that im- 
plicitly appear in Eq. 3.32 have temperature dependences 
that can be extracted from measurements on bulk mate- 
rials. 

Note also that this correlation length is longer than the 
non-linear crossover length for A < a (i.e., for large B). 
For A >> a (small B), C{z) reaches a? before z reaches 
and hence anharmonic effects are unimportant. In 
this case, the correlation length can be determined in 
the harmonic theory (which amounts to evaluating the 
integral in Eq.S.31 with if(k), -B(k), and A?i(k) replaced 
by their constant (bare) values K, B, and Ah). As shown 
in Sec.0, this gives = a^BK/Ah = {a/Xf£,^^, 
which is, reassuringly, much less than in the limit 
a << A in which we have asserted it applies. 

In this section, in our treatment of elastic nonlineari- 
ties we have ignored dislocations, which in Scc^we have 
shown to unbind in a harmonic theory, for arbitrarily 
weak disorder. However, as we have shown in that sec- 
tion, their effects are only felt on length scales greater 
than ~ (a^/A)e'^^i , which is much longer, in the 
weak disorder (A;i — > 0, — > oo) limit, than the trans- 
lational correlation length found here. However, it 
is, nevertheless important to study the simultaneous ef- 
fects of a noma lous elasticity and dislocations. We do 
this in Sec. VII by reconsidering the criterion for disloca- 
tion unbinding, studied in Sec.^, now in the presence of 
anomalous elasticity. We will find that anomalous elastic- 
ity leads to shortening of , which nonetheless remains 
much longer than Hence, for low density aerogels, 
our predictions for and the anomalous elasticity stud- 
ied here remain valid. 

Finally, our entire discussion so far has focused only 
on the effects of orientational tilt disorder. In Sec. VIII 
we will consider the translational disorder (i.e., random 
pinning of the positions of the layers) and will demon- 
strate that it is less important, at long wavelengths (in 
d < 5), than the orientational disorder. Thus, the results 
described herein are directly applicable to real smectics, 
where both kinds of disorder are present. 



VII. EFFECTS OF ANOMALOUS ELASTICITY 
ON DISLOCATION UNBINDING AND SMECTIC 
CORRELATIONS 



In Section ^ we established that, at least at the 
"phonon-only" level, the anharmonic terms in the elas- 
tic energy are extremely important below c? = 5. We 
therefore need to also consider the effects of anharmonic 
elasticity (and the resulting anomalous elasticity) on dis- 
locations and smectic correlation functions. 



A. Dislocations and orientational order 

In this subsection wc incorporate anomalous elasticity 
into our duality model and use it to derive a criterion, 
in terms of the anomalous exponents r]B and r]K, for the 
dislocation unbinding to take place in the presence of 
anomalous elasticity. Combining these results with an 
analysis of orientational order we establish bounds on r]B 
and ijK for the stability of the smectic Bragg glass phase. 
However, since we do not know the numerical values of 
1]B and r]K in c? = 3, we cannot say whether or not, in 
the full, anharmonic theory, dislocations are, in fact, un- 
bound. 

Unfortunately, a full anharmonic theory of dislocations 
is simply intractable. In particular, the fact that, in an 
anharmonic theory, the interaction energy between dis- 
locations cannot be written as a sum of pairwise inter- 
actions (since their fields do not simply add) makes it 
impossible to even write down a general expression for 
the energy of an arbitrary dislocation configuration. At 
best, one might hope to be able to write down the energy 
for a few simple, high symmetry configurations (e.g., a 
single, straight dislocation line). Such specialized results 
would be of no use in a full statistical theory of defect 
unbinding, which requires consideration of very compli- 
cated, tangled configurations of dislocations, which, for 
entropic reasons, dominate the free energy near the dis- 
location binding transition. 

Furthermore, even if one could write down the anhar- 
monic energy for an arbitrary dislocation configuration, 
it would presumably be anharmonic in the dislocation 
fields m(r), and hence, those fields could not be decou- 
pled by a simple Hubbard-Stratonovich trick, as they can 
in the harmonic approximation. 

For all of these reasons, a completely honest treatment 
of dislocations in the full, anharmonic model is impos- 
sible, at least for mere mortals like us. Rather than 
despair, however, we will use a slightly dishonest (i.e., 
uncontrolled) approximation: we will simply replace the 
elastic constants K and _B, and tilt disorder variance 



in our expression Eq.5.37 for the renormalized dislocation 



unbinding transition temperature, with the renormalized, 
wavevector dependent moduli K{c\), B(q), and A;i(q), 
derived from the anharmonic, tilt-only elastic theory of 
Seclvi. 



Doing so, we obtain for the downward renormalization 
of the dual transition temperature 



St = -0(1) X 



2^2 



A,M)qidU 



2 ' 



(7.1) 



where A;i(q) is given by Eq. |6.26c , c is given by Eq. 5.29 
and A^(q) = i ^(q)/-B (q), with K{(\) and B{q) given by 



Eqs. 6.26ci| and 6.26b , respectively. In writing Eq.7.1 



have anticipated that the integral i s do minated by region 
^2 << gj^, as can be seen from Eq j7.3| , below. 

Imposing an infr a-re d cutoff qi_ > on the wavevec- 
tor integral in Eq,7.1, where L is the lateral (_L) real- 
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space linear extent of the system, and changing variables 
in that integral as follows: 

L ' 



A 



LC (CD 



(7.2b) 



(5i « - 



L 
Tnl 



(7.2a) with x = {vb + ?7/-s')/2, wc find that 



AT2 



LA>Qi>l 



n 2 



/B(Q./Qi) 



(7.3) 



with 



7 = C + ?/A + 1 



(7.4) 



In deriving Eq.7.4, we have used the scaling relation 
C = 2 - (rys + 77^)72, and Eq.|!3^ defined in Sec.0. 

Now, it is equally straightforward to show, by 
anisotropic power counting, that the integral over Q in 
Eq. 7.3 is convergent in the ultra-violet when 7 > 0, and is 
therefore independent of its upper cutoff LA, for L — > 00. 
As a result, for 7 > 0, the correction to the reduced dual 
temperature t depends on L only through the overall prcf- 
actor L'^ , i.e. 



5t : 



L 
Tnl 



ca'^cPAh 

Ar2 



xO(I) 



(7.5) 



and therefore, for 7 > 0, diverges with system size. 

As discussed in our original treatment of duality, such 
an infinite downward renormalization of t in the dual 
model implies that dislocations are always unbound, even 
at zero temperature, even for arbitrarily weak disorder. 
Thus, our criterion for determining whether or not dis- 
locations are always unbound in d = 3 is very simple: 
if 7, as given by Eq.7.4, is > 0, dislocations are always 
unbound, and no stable Bragg glass phase exists, while 
if 7 < 0, dislocations will be bound, and a stable Bragg 
glass phase will exist, at least at sufficiently low temper- 
atures and weak disorder. 

So it behooves us to determine 7 in = 3. Using 
C ==_2 — {riB + riK)/2, and the exact scaling relation 
Eq. |6.30| between rjA^rjB, and rjK, we can rewrite Eq.7.4 
in terms of rjA and tjk as: 



7 = 2 [t^k - 1) 



(7.6) 



Using our e-expansion result Eq. |1.10a for rjK , we obtain 



7 = 2(^-1) +0(^2) 



<0, 



(7.7a) 
(7.7b) 



the last equality holding in d = 3 (e = 2), if we drop 
the unknown 0{e^) term. This result appears to imply 
that dislocations are bound, and the smectic Bragg glass 
is stable in three dimensions. Of course, it is always risky 
to quantitatively extrapolate an expansion for small e out 



to e = 2, but this is the best we can do in d = 3. And it 
seems to suggest that the smectic Bragg glass is stable in 
d = 3. 

Independent of the e-expansion, our result implies the 
strict necessary, but not sufficient condition for the sta- 
bility of the smectic Bragg glass phase 



r]K <1 ■ 



(7.8) 



If, on the other hand, this condition is violated, then 
dislocations will unbind on length sca les longer than , 
where is the value of L in Eq.7.5 such that the cor- 



rection St becomes of 0{t). This gives the value of£ 



quoted in the Introduction, with c(T) given by Eq. 5.29 
We note that, although A(T) ^ cxd as T ^ Tna from be- 
low, so does c(r), since the renormalized dislocation core 
energy vanishes at T^a- We expect the divergence of 
c(r) to overwhelm that of A(T), leading to a dislocation 
length that gets ever shorter as T — > T^a from below. 
Event ually, as we approach T^a from below, as given 



by Eq.1.18, will get smaller than , at this point, we 



are in the strong coupling regime, and our weak-disorder 
results no longer apply. 

There is another criterion in addition to Eq j7.§| , that 
must be satisfied for the smectic Bragg glass to be sta- 
ble in d = 3: long-ranged orientational order must exist. 
This is clear, since we have implicitly assumed the exis- 
tence of such order throughout our calculation (e.g., in 
writing n — z + Sn and assuming \dn\ 1). Surely, if the 
layer orientations are not ordered, the layers themselves 
cannot retain their integrity (i.e., remain undislocated). 
So long-ranged orientational order must be preserved for 
the smectic Bragg glass phase to be stable. 

An obvious corol lary of thi s conclusion is that orienta- 
tional fluctuations (|(5n(r)p) must remain finite (as the 
system size L — > 00) for the smectic Bragg glass to be 
stable. 

We will now show, that, within a 5 — d = e-expansion, 
this condition is definitely not satisfied: orientational 
fluctuations do diverge in d = 3. Hence, the e-expan sion 
implies that (notwithstanding the result 7 < 0, Eq. 7.7) 
the smectic Bragg glass phase is not stable in three di- 
mensions. 

To show this, we use the fact that Sn{r) ~ Vj^u(r) to 
write 
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(|5n(r)P) 



(2^) 



(2^) 
d'^q 



^9i(|u(q)P), 

A/.(q)9l 



(2^)'^ [i3(q)g2+X(q)qi] 



(7.9a) 
(7.9b) 
(7.9c) 



where we have kept onl y the dominant, disorder in- 
duced term in (|u(q)p), Eq.3.12,in writing the last equal- 
ity. Imposing an infrared cutoff q± > L^^, and power- 
counting on this integral exactly as we just did in evalu- 
ating Eq.7.3, we find that 



|(5n(r) 



cx L" 



with 



-d-3 + i]A + 2i]B + 3C 



(7.10) 



(7.11) 



Usi ng C = 2 — {rjB + tik) /2 and the exact scaling relation 
Eq. |6.30| , we can rewrite this as 



= 2(X-1). 



(7.12a) 
(7.12b) 



Thus, stability of long-ranged orientational order, and 
hence, stability of the smectic Bragg glass phase itself, 
requires (5 < 0, or 



i]B +riK <2 



(7.13) 



which, not coincidentally, is equivalent to the condition 
that the anisotropy exponent C > 1. The reason that 
these two conditions are equivalent is clear. As we go 
down in dimension, rjB and rjK grow, and the anisotropy 
exponent ( decreases below 2, with the system becom- 
ing more isotropic. The dimension at which the orien- 
tational order disappears must be exactly the dimension 
at which the system becomes completely isotropic, i.e. 
C = 1, since, without orientational order, there can be 
no distinction between different directions. Thus, the 
lower critical dimension for orientational order is given 
by C = 2 — {r]B +r]K)/'2 = 1, which impl ies tjb +'I]k = 2, 



the borderline of the i nequa lity Eq.7.13 



If the condition Eq. |7.13| is violated, then for length 
scales beyond the orientational correlation length, ^o, 



e — 



NL \ 2(x~i) 



A 



(7.14) 



obtained by equating (|(5n(r)P) to 1, the orientational 
order is destroyed, our smectic description breaks down, 
no stable SBG exists, and the s ystem is isotropic. On the 
other hand, if condition in Eq. 7.13 holds, i.e., X < 1; Co 
is infinite and SBG is stable. 

The length scales ^l'" and ^o, derived in Sec.^ and 
in this section, respectively, lead to nontrivial, and in 
principle experimentally measurable crossover behavior. 



At the shortest length scales, q± >> 1/^^^, nonlinear 
elasticity is unimportant and the dominant modes corre- 
spond to (?z ^ M\- At longer scales, for gj_ < this 
quadratic behavior crosses over to qzS,^^ ~ ('ZJ-Ci^)^; 
controlled by the nontrivial anomalous elasticity fixed 
point studied in Sec. VI . Th e anisotropy exponent ( sat- 
isfies 1 < C < 2 if Eq.7.13 is satisfied and orientational 
order is stable, but is less than 1 if this condition is vio- 
lated and Co is finite. In this latter case, the anomalous 
elasticity scaling qz^^^ {q±^±^)'' crosses over to the 
isotropic scaling q^ ^ q± for lengths scales longer than 
^o- This orientational crossover length scale is, reassur- 
ingly, precisely the orientational correlation length given 
in Eq.7.14, obtained from the condition of (|(5n(r)p) = 1. 
It is important to note, that even for ^ < 1, in the domi- 
nant regime of wavevectors for |q| > 1/Co: the condition 
qz << q± is always satisfied. Hence, even in the absence 
of long ranged orientational order beyond Coj this justi- 
fies our use, throughout the smectic regime (|q| > l/Co)j 
of the approximation qz << q_\_, which is indeed the usual 
approximation one makes in treating smectic elasticity. 

These dominant wavevector regimes in q±,qz plane, 
various important length scales, and the corresponding 
crossovers are illustrated in Fig. 11. 




Fig. 1 1 . Crossover in the dominant wavevector regimes in 
q±i Qz plane, for C < 1- 



Combining the stability conditions Eq. 7.8 and Eq.7.13 
with the rigorous bound 



rjB + 5r]K > 4 



(7.15) 



that follows from requiring 77^ > in the exact scal- 
ing relation Eq. 3.3C| in c? = 3, we find that the three- 
dimensional smectic Bragg glass phase can only be stable 
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when rjB and rjK lie in the shaded q uadril atera l in F ig.!. 

Our e-expansion results from Eqs. |6.25"c and 3.28c , ex- 
trapolated to d = 3 



12 

VB = — , 

5 



(7.16a) 
(7.16b) 



lie well outside this stable region. Whether this extrap- 
olation can be trusted in d = 3 is, of course, an open 
question. We are therefore left, despite our best efforts 
as displayed in this manuscript, unable to decide whether 
or not the smectic Bragg glass phase actually exists in 
d = 3. This question can therefore only be answered by 
experiments. 



B. Smectic correlation functions 

In addition to having important effects on dislocation 
unbinding, as demonstrated in the previous subsection, 
anharmonic elasticity also changes the smectic correla- 
tion functions from those derived for the elastically har- 

We study these effects in this 



monic model of Sec. IV 
subsection. 

We focus on the smectic lay ers displacement correla- 
tion function defined by Eq. 



(7.17) 



To compute C(rj_,2:), we proceed as in Sec, IV, except 
that instead of using the bare values of elastic mod- 
uli K, B and tilt-disorder variance A;i, we will use 
their wavevector-dependent counterparts, -ft'(k), -B(k) 
and Aft,(k), given in Eqs.6.26. We expect this proce- 



dure to take into account the long length scale effects of 
anharmonic elasticity. We obtain: 



C(ri,z) 



d'^k 2(l-cos(k-r))Afe(k)fcj 
(27r)'' [X(k)fc4 -t-B(k)fc2]^ 



(7.18) 



where we have only kept the dominant disorder-induced 
contribution. Performing an asymptotic analysis simi- 
lar to that of the previous subsection, we find that the 
harmonic elasticity result, Eq.4.9 is replaced by 



C(ri,z)« 



(7.19) 



valid for d < 5 and on scales beyond the crossover scales 
' when the smectic elasticity is controlled by the non- 
trivial anomalous fixed point with nonvanishing rjB,riK, 
and the anisotropy exponent ^ = 2 — (77^ + Vk)/"^ < 2. 

The above result for smectic layer undulations has ob- 
vious imphcations for the smectic order parameter {ip{r j) 



correlations. Within the Debye- Waller approximation, 
which amounts to ignoring higher order cumulants, the 
smectic correlation function is given by 



(7.20a) 
(7.20b) 



It is short-ranged, confirming our earlier claims about 
the destruction (in d < 5) of the long-range smectic order 
by an arbitrarily small amount of tilt disorder, even in 
the presence of anharmonic elasticity. 

Although the anomalous elasticity does not prevent 
the destruction of the long-ranged smectic order, it does 
lead to some physically appealing features in the above 
expression for 5'(r_L, z), and C(r), Eq, 7.1£ , which are ab- 
sent when these elastic anharmonic effects are not taken 
into account. In particular, we observe that at the lower- 
critical dimension in which the orientational order be- 
comes short-ranged, (i.e. rjs + rjK ^ "2 and C 1), the 
smectic correlation function S(r) is isotropic and decays 
as a Gaussian along both z and r± . In this case of short- 
ranged orientational order this faster than exponential 
decay of smectic correlations of course only persists up 
to the orientational correlation length, beyond which we 
expect the order to decay simply exponentially. 



VIII. ANOMALOUS ELASTICITY AND 
TRANSLATIONAL PERIODIC DISORDER IN 

3 < D < 7 

The analysis of nonlinear elasticity of the preceding 
two sections was confined to the tilt-only model in which 
the random field disorder is neglected. In this section 
we study the effect of the periodic field disorder in the 
presence of anharmonic elasticity and show that, for all 
d < 5, and, in particular in d = 3, such random-field 
disorder 



1 
T 



d''r\tl;o\U{r)cos[qo{u{r) + z)] (8.1) 



has no effect on the long- wavelength "pure phonon" tilt- 
only disorder theory; that is, as long as dislocations are 
bound, the anomalous elastic theory of the previous sec- 
tion is correct in d < 5. In the next section (Sec. IX), we 



show that periodic disorder does not induce dislocation 
unbinding. 

We begin by showing that the periodic disorder it- 
self induces anomalous elasticity in all spatial dimensions 
d < 7. Fortunately, the induced corrections to B and K, 
and to the w-fluctuations themselves, are so weak (only 
logarithmic) in all spatial dimensions d < 7 that, as soon 
as d < 5, they are completely overwhelmed by the power- 
law divergent fluctuations and renormalizations induced 
by the tilt disorder that we studied in Sec. VI. These 
latter fluctuations and renormalizations therefore domi- 
nate, and are unaffected by the presence of the periodic 
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random field disorder. Therefore, once again, the static 
correlations that we derived in Sec. V] from the 5 — e ex- 



pansion on the tilt-only disorder model are correct, in the 
long- wavelength limit, even for the full model. 

To derive these results, we consider the full disordered 
smectic phonon model. In addition to the two types of 
disorder (random tilt and periodic random field disor- 
ders), described by the effective Hamiltonian in Eq. |3.8t 
this model also includes the anharmonic elasticity (con- 
sidered in the previous section in the tilt-only phonon 
model), but still ignores unbound dislocations. General- 
izing Eq j3.8| to anharmonic elasticity, as we did in Sec. VI 
for the tilt-only model, we obtain the effective Hamilto- 
nian for such a phonon model of disordered smectic: 

1 " /A 

— ^ ( — ^|V_L(Ua - U/3)P - AvCOs[qo{Ua - Ufj)] 



T ^ \ A 

Q,/3 = l 



To study the long distance behavior of this model in 
general dimensions d, for reasons that we now explain, 
we must first generalize it further. 

First, recall that in the Sec, III derivation of the model. 



the periodicity of the smectic phase required the transla- 
tional disorder (i.e., the random field Ay) to be a func- 
tion that is invariant under discrete translations by smec- 
tic lattice constant a (i.e., periodic with period a). In 
writing down Eq.3.2, we took this function to be the 
lowest cosine harmonic. Our implicit justification for this 
was that in d = 3, in harmonic clastic theory, the eigen- 
value A„ for the nth harmonic (with wavevector nqo) 



d'^v\i>o\U{r) cos[n(7o(w(i-) + z)] (8.3) 



is given by a simple generalization of 
Ai =4-7? , 



(see Eq.3.21) to be 



A„ = 4 - rjn 



(8.4) 
(8.5) 



H 



E 



f(vi^.«)^ + f(a. 



As long as 77 is fixed under the RG fiow, this shows that 
the higher harmonics of the random pinning potential are 
less relevant, i.e., they are irrelevant at the point at which 
the lowest n = \ harmonic becomes relevant. However, 
here we argue that in d dimensions rj{t) flows according 
to 



(8.2) 



de 



(8.6) 



with the "thermal eigenvalue" 6 > for d > dc with 
dc > 3. We can demonstrate this explicitly in three lim- 
iting cases described below. 



The first is the truncated model Eq.3.8 analyzed in 
SecEn, which ignores anharmonic elasticity. Using the 



exact expressions, Eqs.3.18b 



B (sec discussion in Sec. HI ), an d the constancy of qo 



under the RG, inside r/, Eq.3.19a, we obtain 



3.18c for the flow of K and 



^harmonic 



= rf- 3 



(8.7) 



exactly, in this elastically harmonic model. 

We can also compute near d = 5 in a model defined 
by the Hamiltonian in Eq.B.3 that includes anharmonic 
elasticity but leaves out the random field disorder (Ay), 
which is, as we show bel ow, subdom inant for d < 5. In 
this case, combining Eqs.6.17a, 3.17b| with the flow of qo 
(governed by the eigenvalue (f) = 2 ~ lo) and the fixed 



point value for g2, Eq.6.19, inside Eq.3.19a, we obtain 



= d-3-^ 
2 



(8.8) 



exact to order — 5 — d. 

Finally, as we will show in this section, (with details 
given in Appendix B) for 5 < d < 7, 



W-W = rf-3 + 0(l/^) 



(8.9) 



In these three regimes for d > dc ^ 3 , 6 > 0, which 
leads to ri{£ 00) 0. Based on Eq.S.5, the implication 



of this is that the spectrum of random field harmoni cs is 
degenerate, i.e., an infinite class of operators H^^ , Eq. 3.3 



are all equally relevant and must be treated simultane- 
ously. Since these higher harmonics are generated by the 
RG, we are forced to consider a more general model 

1 " /I 

(8.10) 



T ^ V4 

a,/3=l 



where Ah{ua — up) and Ay(itQ — up) are arbitrary pe- 
riodic functions (with period a, where a is the smectic 
layer spacing), describing the field-dependent variances 
of the tilt and translational disorders. Rather than fol- 
lowing the fiow for an infinite number of couplings (one 



for each Fourier component of the random potentials) we 
apply functional RG (FRG) methods developed by D.S. 
Fisher in the context of the random field Ising and XY 
modelsE3 to our study of disordered smectics. 
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Functional RG for 5 < d < 7: irrelevance of 
random-tilt disorder 



We integrate out the short scale phonon modes within 
an infinitesimal momentum shell near the uv cutoff and 
study how the elastic moduli B{£), K{i) and the disorder 
variance function Av{u,£) evolve under this RG trans- 
formation. As discussed in detail in Sec. VI, the tilt dis- 
order Afi (c haracterized by the dimensionless coupling 
52, Eg. 6. 15 ) is irrelevant by simple power-counting for 
d > 5. We will show at the end of Appendix B that this 



conclusion persists, even in the presence of random-field 



disorder, although the argument is subtler. Hence for 
5 < c? < 7 we can focus on the random-field-only model 
of a randomly pinned smectic. We relegate the details of 
the FRG procedure to Appendix B. For 5 < d < 7, the 
results are summarized by the FRG flow equations 



dK 

m 

dB 

m 



d-3 



d-i 



(d2 - I2(i -f 23) 

3A^(0^" 
2 A2(^) 



K, (8.11) 



.12) 



J and 



dtKv{u,t)^eKv{u,l)- 



3(c!2-6d+ll).„ 



2(d2 _ i)y{t) 



A^(0,^)Ay(u,^) 



u,i) 



(8.13) 



where 

Av{u,i) 



AviuJ) , 



(8.14) 



is the appropriate measure of random field disorder, 
e = 7-d, 



ry(£) 



' K{e) 
B{i) ' 



4:TTy^K{e)B{£) 



.15a) 
.15b) 



and primes in Eq. 8.13 indicate a partial derivative with 
respect to u. 

Note that we have written the flow equations for K{£) 
and B(£) in general d, rather than specializing to d = 
7 — e. The reason for this is that for all d > 5 (i.e., 
wherever we can ignore tilt disorder), the recursion rela- 
tions Eqs. ^.ll" and B.12 become asymptotically exact as 
A"(0, £)/X^JT) 0, which, as we will show in a moment, 
it always does for 5 < c? < 7. 

To demonstrate the above assertion we first assume 
(and verif y a posteri ori) that near d — 7, the flow equa- 
tions, Eqs .|8.1iy8l^ lead to a perturbative fixed point 



which, assuming that Ay{Q,£ oo) 
the solution 



with 



X{e) = Xo^l + £/£o , 



Ao = \i£ = 0) , 



A',^,(0) < 0, has 



.17) 



.18a) 



£o 



|A^,(0)|d2 + 6d-13 



(8.18b) 



The main result of this analysis is that at large length 
scales, £ ^ £o, A^(^) grows as a universal function of £ 



X'^i£^oo) 



6d- 



d2 - 1 



(8.19) 



Hence, as asserted above, as X'^{£) grows, the graphi- 
cal correction s to the flow of K{£) and B{£) diminish, 
Eqs. S. 11 , ^.1^ become asymptotically exact, and the third 
term in Eq. B.13 vanishes as £ ^ oo. All this is still sub- 



ject to the assumption that Ay^(O) 0. 

We can begi n demonst rating this latter assumption by 
combining Eqs.B.11,3.12 into a recursion relation for ?;(£), 
given by Eq. ^.l5b[ y Comparing the resulti ng equation 



with the generic flow equation for ri{£), Eq,8.(;, we flnd 



at which Av{u,£ oo) = Ay{u) is a universal func- that 9{£) (for £ > ^o) is given by a universal function 
tion of order e with a flnite (negative) second derivative 

ilj8.12| into 



A',^,(0) « -g§0(e). Then, combining Eqs.p.ll 
a recursion relation for X{£), as defined by Eq 
obtain 



B.15a, we 



e{£) = d - 3 



d^ -3d+5\ 1 
d2 + 6fi - 13 / 1 



.20) 



dX 

d£ 



A'^(0,£) /d2 + 6d- 13 



X{£) 



2(d2 _ 1) 



(8.16) 



We note that besides being universal, 9{£) is furthermore 
independent of the fixed poin t valu e Ay^(O), as long as 



it is finite. As a result of Eq.8.20 we conclude that for 
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d > 3, r]{i) is strongly irrelevant. Physically, this re- 
flects the strong irrelevance of thermal fluctuations for 
d > 3, with the properties of the pinned smectic de- 
termined by a competition between elastic and random 
pinning energies (both infinitely larger than the thermal 
energy ksT). As a result, the second term on the right 
hand side of Eq. pTT^ flows to zero exponentially fast, and 
can therefore be dropped. Likewise, the 1/A(^)^ term in 
Eq. ^.l3 also vanishes as £ ^ oo (albeit more slowly) , and 
can also be neglected. This leads to a closed recursion 
relation for Av'(u,^), which has been extensively studied 
in the context of the raajiDm field XY modelEa and other 
pinned periodic mediaEliO. 

The fixed point solution A* (u) is readily found by set- 
ting 9f Av(m, t) = 0, which leads to an ordinary differen- 
tial equation for Av'*(m): 

eKv,{u) + i [KAu)f - KMKM = . (8.21) 

The exact form of Av*{u) is not important for our 
conclusions (except for one exponent); all that matters is 
that a fixed point solution Av*{u) exists (which it does 
for 5 < d < 7). However, for completeness we summarize 
the fixed-point results for Ay*(u). Focusing, for conve- 
nience, on the random field correlator A {u) = Ay ^ (u) , 
which satisfies 

cAp, (u) + ( A^, iu)y+ A';, (u) [Af* (u) - Ap, (0)] = , 

(8.22) 



we obtain 

Af^{u) 



nez Qq [6 



ell X n2 

mm — <; -((jou - (2n + l)7r) -- 



(8.23) 



which has cupes at u = n27r/go = na, n integer. It is 
easy to showES that at T = (77 = 0) the cusp non- 
analyticity in Af*{u) develops at finite "time" £, while 
for finite T (77 ^ 0) this singularity is cutoff within 
a boundary layer of thickness ~ ri{£)/£qo at the value 
{0,i) ~ e^ir'^ /9ri{£). We also see from the solution for 



A" 



Af*{u) = A'^,(u), Eq.m, that A^,(0) = -eTr^gq^ < 



0, thereby, together with Eq,8.1£, demonstrating that 
X^{i 00) — > 00, as assumed earlier. 

Having established the existence of the zero tempera- 
ture fixed point for 5 < d < 7, characterized by the uni- 



we now 



versal form of the disorder correlator, Eq. 
turn to the implications of these results for the elastic 
properties of the randoml y pinned smectic. 

Using Eqs. |8.15a - 8.18b insi de the flow equations for 
K{e) and B{1), Eqs. |8.11| , |8.l"2| , and solving for these elas- 
tic moduli, we obtain 



;.24a) 
.24b) 



where the physically measurable moduli Bp {£) and Kp {£) 

(8.25a) 

(8.25b) 
with 



Kpi£) ^Ko[l + - 

£ ^ ^'^-^ 
Bpi£) = Bo(l + - 



, „ 1 /-d2 + i2d-23 
^ 2 I d2 + 6d-13 



and 



3 / d^-1 
2 V d^ + 6d- 13 7 



(8.26) 



.27) 



In the range of interest the exponents (d) and 7b (d) 
are monotonically decreasing and increasing functions 
of d, respectively, and do not change significantly for 
5 < d < 7: 7^(5) = 1/7, 7k(7) = 1/13, 7s(5) = 6/7, 
7b(7) = 12/13. 

As advertised, we therefore find (for 5 < d < 7) 
that the randomly pinned smectic exhibits logarithmi- 
cally weak anomalous elasticity, with the bend modulus 
Kp diverging and the compression modulus Bp vanishing 
logarithmically at large length scales ^ 00. 

With these results in hand, we can calculate the cor- 
responding smectic time-persistent correlation function 



^ (u(k))(u(kO) 

,5'*(k-fk') 



(8.28) 



which determines the long-range spatial smectic correla- 
tions. Following the RG matching method, described in 
detail in Sec. VI, we relate C(k) at the small wavevectors 
of interest to us, to the same quantity computed in the 
effective theory with rescaled ^-dependent couplings and 
evaluated at larger wavevectors e^A:_L, e'^^fczi 

C(k) = e('^-i+'")^C7(e^fc^,e'^^fc„B(£),if(^), Av(m,£)) . 

(8.29) 

Here let us consider the limit kz << Afc^. Now, we set 
(for convenience) ui = 2 and choose = A/k±. The lat- 
ter condition guarantees that no infra-red divergences ap- 
pear i n the perturbative calculation of the right hand side 
of Eq, 8.29| (since it is done at large wavevector k± = A 
at the Brillouin zone boundary). At these short length 
scales (A""'^), phonon displacements u are genuinely small 
and the random potential variance Ay (uq— U/3) (Eq. ^.lC ) 



can be safely expanded in a convergent pow er series in 
Ua — up. Expanding the Hamiltonian Eq. |8.10| to the low- 
est (quadratic) nontrivial order maps the full nonlinear 
FRG problem at these high wavevectors onto the ran- 
dom tilt-only model Eq.5.3 with effective fc-dependent 
couphngs K{£*), B{£*) and A'f{£*) given by 
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Af(r 



-2A'^(0,r)/A2 



(8.30) 



(Ay(0,£*) < 0, as we show below). Using this effec- 
tive theor y and Eqs. 3.11 3.12| to compute the right hand 
side of Eq.8.29, and taking k± sufficiently small, so that 
£* is sufficiently large to guarantee that the fixed point 
Eq. ^.23 is reached, we readily obtain 



C(k) 



-2e(''+i)^*A'^(0,r 



(8.31a) 



[Kp{i*)ki + Bpi£*)k^f Cd-i 



(8.31b) 



To obtain Eg. 8.31b above, we used Eq. 
to Ay and Eq j8.25| to express our resu 



8.14 to relate A 



physical elastic moduli Kp and Bp 



V 

ts in terms of the 
Comparing our final 
expression for C (k) with the expression for this quantity 
in Gaussian theory Eq.3.12 justifies our calling Kp(£*) 
and Bp{£*) the "physical" K and B of the nonlinear the- 
ory of the randomly pinned smectic. Using Eqs. B. 25 for 
these, in the limit £* » £o, with £* = ln(A/fcj_), and 
kz « \k\ we find the anomalous moduli 



if(k) oc |ln(fci/A)r-('') , 
B(k) cx |ln(fc_L/A)r''''^''^ 



1.32a) 
.32b) 



Similar analysis f or a n arbitrary direction in fcj^ — kz 
space leads to Eqs. 1.8, advertised in the Introduction, 
Sec.| 

The Fourier transform of C(k) measures the spatial 
phonon correlations 



C(y)^{[u{ruz)-u[G^M ) , 
d'^k 



(27r)' 



^2(1 cos(k • r))C(k) 



(8.33) 



Using Eq.S.31t, but ignoring the subdominant (In(lnr)) 
corrections arising from weak anomalous elasticity, we 
obtain for z — Q 



C(r^,0) 



-8A^^^,(0) 



k±<r 



(27r)'^ [X(k)fci + S(k)fc2]2 



-4A^y.(0) 



k^<r-' (27r)'i-i t 



d-1 



(8.34) 



We therefore find that throughout the range of dimen- 
sions 5 < d < 7, the smectic displacements u grow 
logarithmically (up to In In corrections due to anoma- 
lous el astic ity), with a universal amplitude determined 
by Eq. 8.23| . This super-universal logarithmic roughness 
is analogou&,to that previously found for the random field 
XY modeled and foij-.the Abrikosov vortex lattice in a 
dirty superconductoia. 

The super-universal smectic Bragg glass phase found 
in this section is quite exotic. However, since it only ex- 
ists for 5 < c? < 7, it unfortunately has little experimental 
relevance for physical 3-dimensional disordered smectics. 
Furthermore, because the periodic disorder-driven layer 
fluctuations and the related anomalous elasticity diverge 
only logarithmically^ while the tilt disorder, which be- 
comes relevant for d < 5, leads to power-law roughness, 
Eq.4.9, the former is completely swamped by the latter 



in d < 5. 



B. Functional RG for 3 < ci < 5: irrelevance of 
random-field disorder 



The above claim can readily be demonstrated by ex- 
tending the functional renormalization group analysis 
presented above and in Appendix B to c? < 5. The re- 
sulting set of flow equations for K{£), B{£), and Ah{u,£), 
val id for d < 5 (assuming the validity of the elastic model, 
Eq.B.lO, i.e., that dislocations are bound) is 



dK{£) 
d£ 



dB{£) 
dt 



1 , , (d^ - l2d + 23) A" (0,^) 
-^-^+32^("'^)+ ^ 2(d^-l)W) 



B 



K, 

;.35) 

i.36) 



deAh{u,£) =(^d-l + ^9{u,i)^ A,,{u,£) + + D{u)^ A'/^{u,£) + c (A'{^,{u)f ^ K{£YB{£)K^-'' /Ca-i 



..37) 



where we have defined 

1 /2 

g{u, £) = Au{u, £) ( Cd-iA'"' 



^K^{£) 

the "position-dependent diffusivity" 



D{u)^C2 (a'^,(u)-A'^JO) 



(8.39) 



■'^^) and c and C2 are dimensionless constants of order unity. 
Note that these recursion relations reduce to those of 



the tilt-only model considered in Sec.Vl, if we take the 
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tilt disorder Ah{u,£) to be a constant independent of u, 
and drop the random field disorder. In these relations, 
for convenience, we have chosen the anisotropy exponent 
w = 2, so that, as discussed earlier, qo does not renor- 
malize. 

The recursion relation for the dimensionless random 



field disorder Ay is given by Eq, 8.13 . Note that since we 
are now interested in d < 5, Ay may grow quite large 
(since its linear eigenvalue e = 7 — (i is no longer small), so 
we cannot neglect the 0(Ay) non-linearities in Eq. ^.l3[ 
Because we do not know these, we cannot actually calcu- 
late the fixed point function Ay*(u) that Av{u,£) flows 
to upon renormalization. However, as in d > 5, we do 
not need to. Rather, we only need to argue that it ex- 
ists, is finite, and that its second derivative at the origin, 
Ay^(O) is finite and negative. 

To make this argument, we first note that if we drop 
the seco nd and third terms on the right hand side of 
Eq. ^.13 and ignore higher order nonlinear terms, the 
resulting equation certainly has a fix ed point solution; 
namely, our explicit solution Eq. 8.23 . Now let us con- 
sider the effects of the other terms. The ri{£) term can 
only reduce Ay, being a "diffusive" type of term. So, 
while it may drive Av{u,£) to zero as ^ — > oo, it cer- 
tainly cannot have the opposite effect of driving Ay to 
infinity. Hence Ay*(u), and therefore, Ay^(O), remain 
finite in the presence of this term. 

If we assume that Ay*(u) is finite, we can show a pos- 
teriori that the third term renormalizes exponentially (in 
renormalization group "time" £) to zero, since, as we will 
show in a moment, A(^) — > oo exponentially in £. 

This should be contrasted with the behavior of this 
term in 5 < < 7, where it vanished much more slowly 
(like l/£) as £ oo. The more rapid vanishing of this 
term in o? < 5 is a consequence of the stronger renormal- 
ization of B and K due to the tilt disorder field (A/j), 



which becomes relevant for d < 5. 

The crucial point here is that this much stronger renor- 
malization of B and K due to this relevance of the tilt 
field (A/i) does not destabilize the Ay fixed point. Since 
this relevance of the tilt field is the only qualitative differ- 
ence between d > 5 and d < 5, showing that it does not 
destabilize the Ay fixed point by itself shows that there 
is no more reason to doubt the stability of this fixed point 
for d < 5 than for 5 < d < 7. Indeed, if anything, the 
argument is strengthened, since the 1/X{£) term vanishes 
even more rapidly for d < 5. 

Hence, the only possible problem can come from the 
unknown higher (than second) order terms in Ay, in 



Eq. S.13 . Not knowing these terms, we cannot make a 
very compelling argument that they do not destabilize 
the fixed point. However, the same objection could have 
been raised for d > 5 but well below 7. Indeed, an anal- 
ogous objection can be raised for any fixed point that 
is found in an e-expansion once one goes well below the 
upper-critical dimension. We simply assume (as is done 
in all standard e-expansions) that the stability of the 
FRG fixed point, rigorously derived for e — > 0, persists 
well below the upper-critical dimension, i.e., for d < 5. 

So it seems plausible to assume that a finite fixed point 
function Av*{u), with finite Ay^{0) < 0, exists for d < 5. 
Let us now investigate the consequences of this assump- 
tion. We will see in a moment that they are virtually 
non-existent: the asymptotic RG flows of B{£), K{£), and 
Ah{£), and long distance behavior of the static smectic 
correlation functions that they imply, are exactly those 
we found in Sec.0, where we neglected the periodic po- 
tential altogether. 

To see this, let us take Ay(u,£) to have flown to its 
fixed point, and combine the recursion relations for B(£), 
K(£), and Ah{u,£) into recursion relations for th e di- 



mensionless coupling function g(u,£), defined in Eq. ^.38 



J which evolves under the RG transformation according to 



5-d - -^giuj) 



C3- 



A" 



(0) 



XH£) 



g{uJ)+['^ + Diu) 

% 



'{u,£) 



KM 



(8.40) 



where we have defined C3 = 3/ 4 -I- (d^ - 12d + 23)/(2(d2 - 
1)). Similarly, from Eqs. |8!35| , p^ , we obtain the recur- 
sion equation for l/\^{£) = B(£)/K{£) 



dX~^{£) 
dl 



--g{(),£)X-\£)+a,A'UO)\"\£), 



..41) 



where ai is a constant whose numerical values is unim- 
portant for our argument. 

We will now prove that for d < 5, at long scales (large 
£), g{u,£) flows to 52 = 32(5 — d)/5, the constant, u- 
indcpcndcnt fixed point value of the dimensionless cou- 
pling g2, Eq.6.15, that we found in our earlier, tilt-only 
treatment of the disordered smectic, described in SccjV^. 



To see this, we write 

g{u,£) = g;+5giu,£) , 
32 

= -{5-d)+6g{u,£) 



and expand the recursion relation Eq. B.40 
in 5g{u,£), obtaining: 

di5g{u, £)=[d-b + D{u)dl] Sg{u, £) + 

where we have defined the "source term" 



(8.42a) 
(8.42b) 

to linear order 



X^£) 



(8.43) 



Siu) ^ c,g*A'{,, (0) -I- c (A- (u)) ' , (8.44) 
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and have dropped the 77(f) term which decays away ex- 
ponentially to as £ ^ 00. 

To demonstrate the stability of the tilt-only u- 
independent fixed point, we now only need to show that 
Eq. ^.43 implies that Sg{u, £) decays to as f — *■ 00. To do 
this we expand Sg(u, €} in eigenfunctions of the operator 
D{u)d^; that is, we write 



ri=0 



gn{£)4>n{u) 



(8.45) 



where (pniu) are periodic functions of u with period a 
satisfying the eigenvalue equation 



Diu) 



= r„0„(M) 



(8.46) 



It is straightforward to see that all of the eigenvalues 
r„ are < 0, provided that D{u) > for all u, as it is 
in 7 — e dimensions, where we can explicitly calculate it. 
We will assume that D{u) continues to be > for d < 5, 
and, in fact, all the way down to the physical dimension 
d = 3. 

Using the orthogonality of the </)„(«) basis, with 



° (j)niu)4)n' (u) 

D{u) 



and 



Jo D(u) 



..47) 



(8.48) 



standard projection analysis leads to a set of decoupled 
recursion relations for the ^-dependent expansion coeffi- 
cients 5„(£) in Eq.B.45 



dgn{i) 
d£ 



(d-5 + r„)g„(i?) + 



where the source term is given by 

InD[u) 



(8.49) 



.50) 



Now, since all of the r„ are < 0, and we are considering 
d < 5, and since, as we show momentarily, X{£) — > 00 as 
£ ^ 00 , we see that all of the gn {€) flo w to zero as f — > 00 , 
provided only that all of the Sn, Eq.8.5Q|-are finite. We 
leave the demonstration of the latter factO as a straight- 
forward homework exercise. Hence as asserted, g{u,£) 
flows to a stable fixed point g*{u) = ^2 = 32(5 — d)/5, 



identical to that found in the tilt-only model of Sec.Vl 



It is now easy to complete the proof, by demonstrating 
via Eq.8.41 that near this fixed point A~^(€) flows rapidly 
to zero: 



X-\£)^X-^{0)e 



(0)e-^' 



;.5Ia) 
.51b) 



Although this result holds only near d = 5, it is straight- 
forward to generalize the above analysis to arbitrary 
d < 5, by relating the graphical corrections to B and 
K to the exact anomalous exponents rjB and rjK of the 
tilt-only model. We find that 172 (^ 00) again flows to 
the same fixed point value it would have in the absence 
of the periodic potential, while A~^(£) now flows to zero 
like 



X-^{£) = A-2(0)e-("«+"'^)^ 



.52) 



So in all dimensions d < 5, X^'^{£) flows exponentially 
to zero, in contrast to 5 < d < 7, where X~'^{£) vanishes 
slowly, like l/£. This completes the a posteriori argu- 
ment in our earlier demonstration that Av{u,£) flows to 
a finite fixed point and that the tilt disorder g(u) flows 
to a stable, M-independent fixed point identical to that 
found for the tilt-only model. 

All that remains to complete our argument that the 
random periodic potential has no effect on the static 
properties of disordered smectics in d < 5 is to consider 
the correlation functions themselves. Using yet again the 
trajectory integral matching formalism, we can show that 



(^(q)^(qO) 

d+l 



.53) 



-)■ 



K{£*)A'^ + B{£*){A/qjLYq'^, 

Afe(r)A2-2A'^(0,r) 
[K{£*)A^ + B{£*){A/q^)^qlf 



(8.54) 

with £* = hi{A/qj_). This only differs from our expres- 
sio n for this correlation function in the tilt-only model, 
Eq.6.31, by the presence of the constant term Ay{0,£*) 
in the numerator of the second term. Using 



^h{£) = g2{£)—, 

-52 



K'^/^£) 



By^(£)Cd-iA'i-5 
K^l'^{£) 



B^/^{£)Cd-iAd-5 



1.55a) 
.55b) 



and 
Av{u,£) ^e- 



4 A 

.00 K^^^i£)B'^\£)Av.{u) , (8.56) 



Cd- 



it is easy to show that the ratio of this Ay (0, £*) term to 
A^(4)A2 is 



A^,(0) Bjl,) _ const. 



9*2 Ki£, 



A (4 



.57) 



which vanishes as (7_l — s- and — > 00, since A~^(^*) 
does. 

To summarize: every possible static effect of the peri- 
odic potential Ay vanishes, as ^ 00, like X^^{£), which 
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vanishes exponentially as i —^ oo. Hence, the periodic po- 
tential (the random field disorder) has no effect on any 
static properties of the disordered smectic in d < 5. Our 
tilt-only model is sufficient to treat those. 

All of the above arguments, of course, only apply if we 
ignore dislocations. In the next section we include these 
in our analysis. 



IX. FAILURE OF PERIODIC POTENTIAL TO 
INDUCE DISLOCATION UNBINDING 



In Sections III-IV and VIII we have demonstrated that 



in three-dimensional (or more generally for d < 5) ran- 
domly pinned smectics, the tilt disorder (A/i) asymp- 
totically leads to significantly larger layer disordering 
than the random field displacement disorder (Ay). We 
thereby argued that static correlation functions for a 
general disordered smectic can be accurately computed 
within a tilt-disorder-only model, studied in Sections ^v[ - 
VII. However, the validity of these arguments and the 



conclusions that follow from them critically rely on the 
stability of the elastic model within which these assess- 
ments are made. That is, it is in principle possible that 
although the periodic disorder is less important than the 
random tilt pinning in disordering smectic layers in the 
elastic model, it could nevertheless drive dislocation un- 
binding, thereby invalidating the elastic model and the 
conclusions drawn from it. 

Therefore, to make the arguments of the previous 
sections complete, it is essential to assess whether the 
random-field disorder (Ay) alone can drive proliferation 
of dislocation loops, and this is the subject of current 
section. |— ■ 

To this end we adapt D.S. Fisher'aSi real space renor- 
malization group analysis to disordered smectics and 
show that the random periodic disorder alone does 
not induce a proliferation of dislocation loops in three- 
dimensional smectics, thereby concluding our justifica- 
tion for the tilt-only model of randomly pinned smectics. 

We first review Fisher's analysis as it applies to the 
three-dimensional random-field XY model. Recall that 
within a defect-free model the elastic energy is balanced 
against the pinning energy, both much larger than the 
thermal energy ksT. The compromise in this competi- 
tion leads to phase variations that grow logarithmically 
with distance and typical ground-state energies that vary 
as , where L is system size and 6 is (minus) of the 
thermal eigenvalue exponent, which, due to statistical 
symmetry, is given exactly by 



d-2 



(9.1) 



We wish to know whether it energetically pays for the 
system to insert a vortex loop. Confining our discus- 
sion to d = 3, a vortex loop of length L will cost the 
system elastic energy E^i ~ JLlnL, where the multi- 
plicative logarithm is associated with the long-range XY 



phase deformation around the vortex line and J is the 
spin-wave stiffness. For such a randomly placed vortex 
loop, which relieves strains by allowing spin rotations of 
order 27r, the system can typically gain disorder energy 
-^disorder ~ ~ L. Hcncc, this crude argument would 
suggest that, because of the additional multiplicative log- 
arithmic factor in E^i^ it never pays to insert arbitrarily 
large vortex loops into the random-field XY model. How- 
ever, as argued by Fisher, a judicious optimization of the 
vortex loop location and conformation can do better en- 
ergetically, and therefore it is essential to include in the 
stability analysis the energy minimization over the vortex 
position and conformation. 

To this end we start out with a circular vortex loop 
and, as illustrated in Fig. 12, try to lower its energy by 
moving various segments of length A^, of it first by a dis- 
tance w, then by 2w, etc. We then calculate the renor- 
malization of the line tension, e(£), by pinning at scale 




Fig. 12. Illustration of the optimization of a topological 
defect loop on length scale A„ . The full and dashed large 
loops represent the topological defect line before and af- 
ter optimization on scale A^,, respectively, and the small 
circles are the isotropic correlated regions corresponding 
to displacements w. 

The idea is that each roughly spherical region of size w 
is statistically independent of all the others; each gains a 
pinning energy of the order — w'^~^ = w (in d = 3). 
So the total pinning energy gained from optimizing the 
conformation of the vortex loop on scale w is given by 

in ^ — Ngpfieres j (9.2a) 

(9.2b) 
(9.2c) 
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Balancing this energy against the stretching energy nec- 
essary to deform the vortex loop 



etch 



gives 



^3/2 



A„ 



„3/2 



(9.3) 



(9.4) 



which predicts the aspect ratio of the energetically most 
favorable deformation of the vortex loop 



A. 



(9.5) 



Using this in the expression for the pinning energy, we 
get the best pinning energy gain of 



(9.6) 



Associating this energetic reduction with the renormal- 
ization of the vortex line tension e{£) gives 



J..^^i±i^.Jln2. 
A^ 



.1/3, 



,2/3 



W 



(9.7) 



where the first term arises from the elastic energy in- 
crease due to the factor of 2 increase in length scale. 
This result implies the RG recursion relation 



de{£) 
di 



J 



1 



;(€)i/3 



(9.8) 



Clearly, if the bare vortex line tension e(0) = e is suffi- 
ciently large, the line tension e{£) at scale grows like £ 
(w InL), and therefore it is energetically unfavorable for 
vortex loops to unbind. This therefore argues for the sta- 
bility of the topologically ordered "Bragg" glass phase in 
a weakly disordered three-dimensional random-field XY 
model. 

We now extend this analysis to smectics pinned by a 
random periodic potential, with the main difference from 
the random-field XY model being the anisotropy of the 
smectic state. We begin by calculating the thermal eigen- 
value exponents —9. As was the case for the random- 
field XY model, statistical symmetry guarantees that the 
smectic elastic moduli B and K are not renormalized by 
disorder .e3 This therefore implies that the 6 exponent can 
be computed exactly by estimating the elastic energy due 
to layer displacement u of order a lattice constant, 0(a), 
inside an anisotropic volume L± x L± x (L^ ~ ^j^)- 



E ^ B 







dzd^r± {dzu) 
1 



LI 



(9.9a) 

(9.9b) 
(9.9c) 



In this estimate Eq j9.9| we have used the fact that, in 
the harmonic smectic elastic theory, distances in the z- 
direction scale like the square of those in the _L directions. 



In three dimensions, we therefore find that the ground 
state of a randomly pinned, elastically harmonic smectic 
is described by 







(9.10) 



We now need to estimate the pinning energy gain from 
introducing a dislocation loop and optimizing its con- 
figuration, analogously to the analysis for a vortex loop 
above. The main difference here, however, is that, be- 
cause of the anisotropy of the smectic phase, there are 
three dislocation line optimization cases to consider: 

1. Vertical (along z) w displacement of a horizontal 
(perpendicular to z) segment (i.e., an edge disloca- 
tion). 

2. Horizontal w displacement of a horizontal segment. 

3. Horizontal w displacement of a vertical segment 
(i.e., a screw dislocation). 



A. Case 1: Vertical displacement of a horizontal 
dislocation (w \\ z) 



Because of the elastic anisotropy of the smectic state, 
in contrast to the spherical correlated volumes of the 
random-field XY model, here correlated regions are cigar- 
shaped (directed along z) characterized by dimensions 
L±(w) X Lj_{w) X w, with Lj_(w) ^ \/w, as shown in 
Fig.l3. 



w 



w 



1/2 



w 



Fig. 13. Illustration of smectic anisotropic cigar-shaped 
correlation regions inside the dislocation optimization 
length Ayj for a displacement w along z of a dislocation 
running perpendicular to z (case 1). 
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The corresponding pinning energy change due to opti- 
mization of a dislocation segment of length is given 
by 




(9.11c 



Balancing this against the dislocation stretching energy, 
which, as always, is 




SE. 



stretch 



we get 



(9.12) 



(9.13) 



We note that A^ >> Lj_{w) ^ ^Jw, consistent with the 
assumption of many "bubbles" per segment , on which 



this calculation of (5i?pin was based. This leads to 



<5i?pin 



A 



1/2 



(9.14a) 
(9.14b) 



and predicts a change of the line tension at scale w given 

by 



5e 



A„ 



pl/3 



w 



(9.15a) 
(9.15b) 



Note that in the above we have ignored the elastic en- 
ergy, because we do not need it to prove the stability 
of smectic "Bragg" glass! Even wit hout the stabilizing 
contribution of the elastic energy, Eq,9.15 implies 



d£ 



1 



£(^)l/3 



(9.16) 



Clearly, if the bare €{£ = 0) = e is big enough, e{£ oo) 
will be non-zero. Hence periodic disorder, by itself, can- 
not unbind horizontal loops by displacing them vertically. 



B. Case 2: Horizontal displacement of a horizontal 
dislocation [w _L z) 

As illustrated in Fig. 14, here, the volume of the cor- 
related regions is w x w x w'^ , with the long axis along 
z. 



Fig. 14. Illustration of smectic anisotropic cigar-shaped 
correlation regions inside the dislocation optimization 
length A„ for both the displacement w and the dislo- 
cation running perpendicular to z, i.e. along r± (e.g. 
along X and y) (case 2). 

The width of the correlated volume along A^, is now 
w, so we have 



SEpi^ 



(9.17) 



Equating this energy to the dislocation stretching energy 
eui^/Aiu, we find 



A„ 



,2/3^5/3 



which, upon inserting into SEpi^, Eq.9.17, gives 
and implies 



5e 



This leads to 

de{t) 
d£ 



A„ 



1 



.l/3,,,4/3 ■ 



e(£)1/3 



= -«/3 



(9.18) 
(9.19) 

(9.20a) 
(9.20b) 

(9.21) 



which also demonstrates that there is no divergent down- 
ward renormalization of e, implying that in the presence 
of only periodic disorder, horizontal dislocation loops will 
not unbind by displacing horizontally. 



C. Case 3: Horizontal displacement of a vertical 
dislocation {w _L z) 

Finally we consider optimizing a z-directed screw dis- 
location, by displacing it horizontally. This corresponds 
to a correlated region of volume w x w x w'^, but, in con- 
trast to Case 2, with the axis along A„, being the long 
dimension w^, as illustrated in Fig. 15. 
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Fig. 15. Illustration of smectic anisotropic cigar-shaped 
correlation regions inside the dislocation optimization 
length for displacement w along r± and dislocation 
nmning along z, i.e. along rj_ (case 3). 



As a result, we have 



SE, 



pin 



(9.22) 



which when balanced against the dislocation stretching 
energy gives 



A. 



(9.23) 



We note that for I this A^, is much greater tha n 
consistent with the assumption that goes into Eq. ^.22 
This gives 



SEp 
A,„ 



which implies 



deje) 
di 



-21 



(9.24a) 
(9.24b) 

(9.25) 



whose solution for e(£), again, is stable as I ^ oo. 

Since in all three cases we find that the effective dis- 
location line tension is finite at long scales, we conclude 
that periodic random-field disorder alone cannot unbind 
dislocation loops in such a randomly pinned smectic. 



EXPERIMENTAL PREDICTIONS AND 
FUTURE RESEARCH DIRECTIONS 



Our results imply many unequivocal dramatic predic- 
tions for experiments. The first and most unequivocal 
of these is that long-ranged, and even quasi-long-ranged, 
smectic translational order is destroyed by the presence 
of even arbitrarily weak quenched disorder. Furthermore, 
the decay of smectic translational order induced by such 
disorder is far more cajjiid—than in previously studied 
"Bragg glass" systemsES'LTtiJ, being exponential, rather 
than algebraic. The experimental signature of this is 
equally unequivocal: broad X-ray scattering peaks with 
finite peak height, as opposed to the infinitely sharp, di- 
vergent peaks associated with quasi-long ranged order in 
thermal smectics. 

This qualitative prediction is born out by all experi- 
ments on smectics in aerogel to date: all show only broad, 
finite height X-ray scattering peaks, which-can be fit by 
assuming a finite X-ray correlation length.t3 

In strong contrast to quenched (frozen) disorder of this 
type, it is straightforward to show analytically that weak 
short-ranged ^^annealed^ disorder, in which the random 
environment can rearrange to accomodate the smectic's 
preference for an ordered state, has no qualitative effect 
on the long scale properties of liquid crystal phases. Such 
weak annealed disorder only leads to finite renormaliza- 
tions of the effective parameters in our model. c3 Conse- 
quently, the smectic phase, its quasi-long-ranged smectic 
translational order, and the transition into it are stable to 
weak annealed disorder, which might arise, for example, 
due to a low concentration of microscopic inclusions. 

Of course, strictly speaking, all physical random con- 
fining structure are elastic, even aerogel, and therefore 
are able to deform to some extent in response to the 
smectic, which resists distortions. It is easy to see that 
this additional "elastic compliance" ingredient can be de- 
scribed as an annealed component of the disorder, in ad- 
dition to the quenched part that is the main subject of 
our work. In light of the qualitative unimportance of 
weak annealed disorder, discussed above, our theory for 
purely quenched random structures immediately extends 
to deformable structures, which are quenched for large 
deformations, i.e. described by a finite shear modulus. 

A notable example of liquid crystals confined in 
random structure in which these elastic annealed fea- 
tures a ppea r are the recent experiments on smectics in 
aerosilPWA Since typical weak hydrogen bonds that link 
the aerosil network together have bond-breaking ener- 
gies of order fc^T, a large fraction of the aerosil network 
can rearrange itself in response to stresses imposed on 
it by the smectic and thereby come ta_thermal equilib- 
rium with it. However, measurementsLj show that even 
these tenuous, weakly-bonded structures are character- 
ized by a small but a finite shear modulus. Hence even 
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aerosils are unable to perform arbitrary rearrangements 
to equilibrate with the ordered smectic. Therefore, while 
they have a large annealed component, aerosils contain a 
small quenched component, which dominates the qualita- 
tive physics of the confined smectic. We therefore expect 
a destruction of the quasi-long-ranged translational or- 
der for smcctics confined in aerosils, with the X-ray cor- 
relation length obeying our predictions for low density 
aerogels. 

However, there are three important consequences 
of the large annealed disorder component present in 
aerosils. Firstly, while weak annealed disorder should 
only lead to a finite renormalization of effective param- 
eters, strong annealed disorder can destroy the quasi- 
long-ranged smectic order by driving some of the effec- 
tive elastic moduli negative. Secondly, given the em- 
pirical fact that the bulk N A, transition appears to dis- 
play non- universal behavior]^ we would expect this non- 
universality to also manifest itself for smectics "confined" 
in low-density aerosils, where the effective model param- 
eters are functions of aerosil density due to the presence 
of the annealed disorder component. Thirdly, because of 
this large annealed disorder component in aerosils and 
a correspondingly small quenched component, we expect 
aerosils to be described by significantly smaller values 
of our effective quenched disorder parameters Ay and 
A/j than for aerogels of the same solid volume fraction. 
Consequently, experiments in aerosils are more appropri- 
ate for exploring our weak quenched disorder theory and 
to search for the delicate and "fragile" SBG and NEG 
phases. 

These nontrivial qualitative predictions are clearly in 
agreement with the experimental observations of smec- 
tics confined in aerosilsEHl: (i) for low aerosil densities 
these experiments observe significantly longer (than in 
equivalent volume fractions aerogels), finite smectic X- 
ray correlation lengths, (ii) they measure heat capacities 
with aerosil density dependent (and therefore nonuniver- 
sal) exponents, and (iii) in aerosils beyond a critical den- 
sity they find a behavior that is qualitatively similar to 
smectics confined in aerogels. 

It is important to note that in aerogels and aerosils, in 
contrast to what one might have naively expected, the 
smectic correlation length is not limited by some sort of 
"finite size" effect associated with the "pore size" of the 
aerogel. That is, the smectic confined in aerogel does not 
act as if it has simply been broken up into many small 
volumes which do not interact with each other. 

If this "finite size" scenario were the case, one would 
expect to see the smectic translational correlation length 
grow with decreasing temperature as T approached 
the pure NA transition temperature Tf^J° from above 
exactly as it did in the pure (bulk) smectic. Indeed, one 
would expect it to track its temperature dependent value 
?purc(^) i'^ the pure system until T got close enough to 



that e^,„(T) 



aerogel 



(T) 



urate, and remain constant at Lp for all lower tempera- 
tures. 

This is ernphatically not what is seen in experiments 
on aerogels. E£lE£fCJ Rather the growth of continues 
right through Tf^^^, and only saturates at some tempera- 
ture well below T^"^°. This is completely consistent with 
our picture, in which (^-^ is determined by a competition 
between the smectic elasticity, which tries to keep the 
system ordered, and the random forces exerted by the 
aerogel strands, which disorder it. In this picture, the 
saturation of with decreasing temperature occurs not 
above T^^"^, but, rather, well below T^^°; specifically, at 
the temperature at which the smectic elastic compression 
modulus B{T) saturates. 

Liquid crystals have also been confined in other ran- 
dom but significantlv_more ordered-structures. Some ex- 
amples are AnoporeEj and VycoiO, which in contrast 
to the tenuous aerogel strand structure present on all 
scales, consist of a random distribution of empty pores, 
with fairly regular wall structure and much narrower pore 
size distribution. Even these systems can in principle 
be understood within the general model presented here. 
However, they fall into the strong disorder regime of our 
model, and we are therefore unable to utilize our weak 
disorder approach to make quantitative predictions, aside 
from the observation that the quasi-long-ranged smectic 
order should be destroyed in these systems as well, con- 
sistent with experiments. 

The qualitative agreement between our theory and ex- 
periment can be made quantitative through our expres- 
sion Eq.1.13, which relates to the smectic elastic con- 
stant B{T). As discussed in the Introduction, this ex- 
pression implies that sufficiently far below 7)^'^° (where 
our weak disorder theory applies), a simple power-law re- 
lation between (T) and the "bare" layer compressional 
modulus i?barc(7') of the pure system holds: 



e^(T) « [i?barc(r)]^/''^ , 



OC [i?bare(r)]^/'^-^/^ 



(10. fa) 
(10.1b) 



the 



pore size 



The bare Share (T) can be determined in a number of 
ways. One way is simply to measure the layer compres- 
sional modulus i?purc(T') of a pure (bulk, aerogel- free) 
sample of the same smectic liquid crystal material, and 
pray that this is not changed when the smectic is placed 
in aerogel. However, since prayers are not always an- 
swered, a direct, simultaneous in-situ determination of 
Share (r) in the smectic inside the aerogel is clearly de- 
sirable. 

This could be done by determining the squared mag- 
nitude of the smectic order parameter 1-01^ from, e.g., the 
integrated intensity under the (broadened) smectic X-ray 
peak. In mean-field theory, which should hold for bare 
quantities sufficiently far below Tjva, Share (T') cx jV'P- 
Thus, we should, in this mean-field regime, observe a 
universal power-law relation between the X-ray correla- 



of the aerogel. At this point ^ (T) would abruptly sat- tion length and \Tp\'^: £, 



iX 



P/^-i. Note that 
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the abse nce of anomalous elasticity, the second regime of 
Eq.|l.l3| (^^ = i^^{a/Xf) would apply, and we would 
get^-^oc B cx IV'P, the last proportionality holding only 
in the mean-field regime. Thus, any departure of the 
S^^-B relation from linearity is evidence for anomalous 
elasticity. Furthermore, observing such a power-law re- 
lation determines %. i— . 

Preliminary analysis of versus B dataEj, as deter- 



mined above, supports the relation Eq.10.1, with x = 
0.8 ±0.1. Note that x < 1, favoring stability of the smec- 
tic Bragg glass phase, although we have, as yet, no exper- 
imental way of determining if the other necessary relation 
for the stability of the SBG, namely tik < 1, is satisfied. 

While our prediction for (T) is not as quantitative as 
we would like, since we do not have quantitatively reliable 
predictions for C, and x in d = 3 (our 5 — d = e-expansion 
results not being quantitatively trustworthy when e = 2), 
it nonetheless makes it possible, by experimentally deter- 
mining Q/x for one smectic in one aerogel, to predict the 
temperature dependence of for any smectic in any low 
density aerogel. Furthermore, once C/x is known, C and 
X separately are also known, since ^ = 2 — x- Knowledge 
of these exponents then makes other predictions possible. 

In particular, it leads to our second prediction which 
concerns whether or not the transition (which replaces 
the NA transition in bulk smectics) from the nematic to 
the low-temperature phase (which replaces the ordered 
smectic) is destroyed by disorder. Surprisingly, this ques- 
tion can, in principle, have a different answer than the 
question of whether or not smectic translational order is 
destroyed. This is because our analysis shows that, even 
though smectic translational order is always destroyed, 
a distinct "smectic Bragg glass" phase (SBG), in which 
dislocations remain bound and long-ranged orientational 
(nematic) order persists, maybe possible in the presence 
of sufficiently weak disorder. If it is, then a thermody- 
namically sharp transition into it (with decreasing tem- 
perature) replaces the NA transition of the pure system, 
as this phase itself replaces the smectic-A phase. Conse- 
quently, an exotic and quite unique feature of this tran- 
sition is that, while thermodynamic quantities (e.g., heat 
capacity) display singularities at the transition to the pu- 
tative smectic Bragg glass, the X-ray scattering peak 
remains broad through the transition. Unfortunately, 
whether or not this (SBG) phase exists, in principle, de- 
pends on whether or not the imprecisely known anoma- 
lous exponents rjB and tjk satisiy the bounds Eq.^. 

Fortunately, the above-described experimental deter- 
mination of X could answer this question, since one of 
the bounds Eq.Lla is simply x < 1- Thus, if the X-ray 
measurements, described above, find x < Ij a SBG phase 
can exist, and a thermodynamically sharp remnant of the 
pure system's NA transition could be seen in sufficiently 
low density aerogels. On the other hand, if x proves to 
be > 1, no sharp transition should be seen. This is a 
definite, falsifiable prediction that could be tested exper- 
imentally: if X > Ij no sharp transition should be seen. 
The converse conclusion, that if x < Ij a sharp transi- 



tion should be seen, need not hold, since it is possible 
that X < Ij but the bound tjk < 1 is violated. 

For the same reason, even if the current experimental 
evidence that there is no transition even for arbitrarily 
low density aerogel holds up, we can not conclude that 
X>1- 

It is also important to note, that despite of the absence 
of the long-ranged ordered smectic phase, it is possible, 
as always, to have a thermodynamically sharp first-order 
transition within the disordered (nematic) phase. At such 
transition, which would be from a nematic with a short 
smectic correlation length to a nematic with longer (but 
still finite) smectic correlation length, thermodynamic 
quantities (like the heat capacity) could displaying sharp 
non-analytic features (like latent heat). We believe that 
the recent calorimetry experiments on 70.4/aerogel sam- 
ples by Haga aiuL Garland are examples of such a first 
order transition.E^ As a consequence of our main result 
that only short range smectic order is possible for smec- 
tics confined in quenched random environments, we un- 
ambiguously predict a broad, finite X-ray scattering peak 
for the system studied by Haga and Garland. 

If the SBG phase does prove to be possible, then we can 
make far more detailed tests of our predictions. In partic- 
ular, the simple fact that the SBG phase has long-ranged 
orientational order is, in itself, a striking prediction that 
could be tested very easily by looking at optical birefrin- 
gence in smectics in aerogel. The existence of such long- 
ranged order would make it possible even more quanti- 
tative tests for our predictions of anomalous elasticity, 
by making it possible to measure director fiuctuations 
{6ni{q)6nj{—q)) , by light scattering. Since these are re- 
lated to smectic layer fiuctuations by the usual smectic 
relation Sn{r) Vu(r), we can immediately write 



((5ri,(q)(5nj(-q)) 



q,gj(|u(q)|2) , 

[B{q)ql+K{q)qiy 
ksTqiqj V 
B{ci)ql+K{ci)qi 



(10.2a) 



(10.2b) 



where i and j run over all indices but z, and V is the 
volume of the system . The first, quenched, term in this 
expression, Eq. l0.2b , dominates as q ^ 0; thus fitting 
this piece to the observed light scattering would enable 
a direct experimental measurement of K{c(), and 
Aft (c[], and, hence, a direct test of the scaling predictions 
Eq.6.26 for these quantities, as well as a determination of 
the anomalous exponents rjx, rjB, a nd 77 A . For instance, 
if we take q in the x — z plane, Eq.10.2 implies 



ilSn^iciW) ^ q- 



4-r;A+2)) 



K f 



NL 



(10.3) 



Thus, a log-log plot of {\5nx{q)\'^) for = has slope 
— 4 — 77A + 277if . The simple observation that this slope 



44 



is not —4 would confirm the existence of anomalous elas- 
ticity, by showing that tja and rjK were 7^ 0. The deter- 
minatio n of the ani sotropy exponent ^ by collapsing the 
plots of (|(5na;(q)p) versus qx for a variety of g^'s would 
then fix rjB + fjK ■ Ou r exact scaling relation between 77A , 
riB and rjx, Eq.6.3C would then provide a third equa- 
tion to fix all 3 exponents. That the results satisfied our 
bounds Eq.1.1 for the stability of the SBG phase would 
then provide a non-trivial check on the consistency of our 
theory. 

A more stringent test could be p rovided by measuring 
the second, ksT term in Eq. 10.2 . Although, as men- 
tioned earlier, this term is subdominant to the first, dis- 
order term, it could nonetheless be resolved by time- 
dependent light scattering. The reason for this is that 



the first, disorder term of Eq.10.2 represents the static, 
time-persistent response of the smectic to the random tilt 
field of the aerogel, while the ksT term represents ther- 
mal fluctuations about this disorder-determined smectic 
layer configuration. That is, if we consider the unequal 
time correlation function {6ni{q, t)dnj(~q, 0)), we expect 

'^hicOqtQjqi V 



(5ni(q, t)Snj{-q, 0)) 



kBTq,qj f{q,t}V 
+ i?(q)<z2+A'(q)q4 



(10.4) 



where /(q, t) is a function we do not know, but which 
presumably will have some sort of slow, "glassy" decay, 
vanishing as t 00, and going to 1, as t 0, thereby 



recovering the equal-time correlation function, Eq.10.2 
Using these two limits immediately implies 



limt^oo {Snt{q,t)Snj {-q,0)) - {Sn^{q,0)Snj{-q,0)) 
__ kBTqiqj V 



S(q)(72 + K{q)qi 



(10.5) 



which allows a completely independent determination of 
rjK and tjb from the approach based strictly on the equal- 
time correlation function described earlier. Consistency 
of the two approaches would demonstrate the validity of 
the exact scaling relation between 77^, rjB, and tja- 

We remind the avid experimentalist that the entire 
above discussion only applies if the SBG phase is stable. 
If the SBG phase does not exist, then our experimental 
predictions are reduced to Eq. l.l3| for the X-ray corrc- 
lati on le ngth, the anomalous elasticity lengths given by 
Eqs.1.3, and the "ghost" of non-analyticity in T in these 
lengths, associated with the remnant of Tp below T^a 
(see Fig. 6), defined by the temperature at which the pure 
smectic order parameter correlation exponent rj passes 4. 

Presumably, there are experiments that could test our 
predictions for the orientational correlation length ^o, 
but we will be unable to make such predictions until we 
have developed a theory of the nematic elastic glass be- 
havior that holds op-ilonger length scales. This is a topic 
for future research .Ea 



There are a number of other promising areas for fu- 
ture research. One is the subject of smectic order ipside 
stretched aerogels, which we are currently studying.E3 In 
that system we have shown that uniaxially stretching or 
compressing the aerogel matrix can stabilize the smectic 
Bragg glass phase. The universality class of the resultant 
SBG now depends upon whether the aerogel is stretched 
or compressed: in the former case, it lies in the universal- 
ity class of the so-called "vortex glass" phase, which has 
recently been Tf^Vfk discussed in the context of random 
field XY-modellS'til and randomly pinned Abrikosov vor- 
tex lattices^ while compression leads to a totally novel 
kind of "m = 1 Lifshitz SBG"rx-|A detailed discussion of 
this problem is in preparation.E3 

We expect a new glass phase similar in many respects 
to the "m=l" glass described above, but in a different 
universality class, to occur when "discotics" (i.e. liquid 
crystals with two solid-like directions and one liquid-like) 
are absorbed into unstressed aerogel. Worlpon this inter- 
esting problem is also currently underway£3 

It is also interesting to consider what hacpens to a 
cholesteric liquid crystal confined in aerogel.u Although 
the symmetry of such a system is identical to that of a 
smectic studied here, the addition of another long length 
scale (namely the cholesteric pitch) may make it possible 
to access a novel type of "random manifold" regime of 
pinned elastic media. 

Probably the most interesting and challenging problem 
is the theory of dynamics of these smectic Bragg glass sys- 
tems, which would enable a detailed understanding of re- 
cent dynamic light scattering experiments. tj If the smec- 
tic Bragg glass phase exists, which we are quite certain it 
does in smectics iitnanisotropic aerogels and discotics in 
isotropic aerogels ,E3 it should have slow glassy dynamics 
similar to those observed by Bellini et altB. Given its 
unusual static properties (e.g. anomalous elasticity), its 
dynamics may be quite unusual, even for glasses. 

As noted in the Introduction, here we have studied 
smectics in aerogel by perturbing around the low tem- 
perature translationally and orientational ordered smec- 
tic phase. However, another complementary approach 
to this problem, more convenient to understanding the 
effects of disorder near and above tiic bulk NA transi- 
tion, is based on de Gennes model.E3 This complemen- 
tary high temperature approach makes it possible to pre- 
dict, e.g., the rounded specific heat and X-ray scatter- 
ing peaks [Observed in many experiments on smectics in 
aerogeLBm 

Finally, we note that the formalism we developed in 
Sec.^ for disordered smectics can be applied to a vari- 
ety of other disordered problems, including gauge glasses, 
randomly pinned Abrikosov vortex lattices and other 
pinned periodic media. These-applications will be dis- 
cussed in a future publication.ca 
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XII. APPENDIX A: DETAILS OF THE D = 3 
SINGLE HARMONIC RENORMALIZATION 
GROUP ANALYSIS 

Here we present the details of the de rivation of t he 
renormalization group flow equations Eqs.3.18a- 3.18b| for 
the three-dimensional disordered smectic model, in the 
harmonic elastic approximation defined by the replicated 
Hamiltonian Eq.|3.8| 



H[Ua\ = i?o[Ua] +ffA[ 



(12.1) 



where the quadratic i/o[ua] and the interaction Hj\[ua\ 
parts of the Hamiltonian are respectively given by 

If " 

(12.2a) 



Ha 



d'^rY, Av cos[go(ua(r) - upir))] . (12.2b) 



In the above and throughout, in the interaction term 
Ha, we have for convenience excluded the a = (3 term 
inside the summation. This exclusion is indicated by a 
prime, and is justified since the a = /3 term is a con- 
stant. Also for simplicity, throughout this appendix we 
have set T = 1; the factors of T can be easily restored 
by the replacement K K/T, B ^ B/T ^h/T^, 
and Ay — > Ay/T^. From Hq, the propagator Gafiifl) 
can be easily obtained. 



Ga/3(q) = GT(q) Sap + Ga^ (q' 



Kql + Bqi {Kql + Bq: 



2\2 



(12.3a) 
(12.3b) 



In three dimensions, a direct perturbative calculation 
in Ay is divergent in the thermodynamic limit even for 



an arbitrarily small Ay. Nevertheless physical observ- 
ables can be computei-utilizing standard methods of the 
renormalization grouped, in which one avoids infrared di- 
vergences by integrating out degrees of freedom a mo- 
mentum "shell" at a time. In this procedure the goal 
is to establish how the effective Hamiltonian functional 
changes after the renormalization group transformation, 
which can be neatly summarized by the renormalization 
group flow equations for the effective parameters. 

We focus on Z" as the physical quantity to keep fixed 
under the infinitesimal renormalization group transfor- 
mation. Following the standard RG procedure we write 
the Fourier-transformed displacement field as Ua{q) = 

(q) + (q) , where (q) is the long wavelength set 
of modes, nonzero for < |qj^| < Ae~^, with q^ unre- 
stricted, and u^(q) are the short wavevector degrees of 
freedom that have support only within a thin cylindrical 
momentum shell 



Ae < |q_L|< A , 
— oo < qz < CO , 



(12.4a) 
(12.4b) 



We first integrate, perturbatively in Ay, the high 
wavevector part out of Z": 



Z"= [du<][du>]e 



<1U„>1 „--H"o[<+«>]-ffA['"^+«S 



<1 g~Hol<]~SH[u<] 



(12.5a) 
(12.5b) 



where SH[Ua] is obtained by integrating out the short 
wavevector degrees of freedom 



-SHlu<] ^ J j^^>j ^-Holu>]~H^lu<+u^ 



1-{Ha[u<+uI]). 



(12.6a) 
(12.6b) 
(12.6c) 



2 (i?A[< + 



The perturbative correction to the non-constant part of 
the effective Hamiltonian can therefore be expressed in 
terms of a cumulant expansion 



1 



5H[u<] = {Ha[u< + ul])l - ± {Ha[u< + u>S)l 



where 



(12.7) 



{Ha[u< + u>]% - {Ha[u< + u>f)^ - {Ha[u< + u>]) 

(12.8; 

To first order in the interaction Ay we obtain 
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— —Ay / d^r ^ cos[go(i 



A^y/ E T:J^'(r-r'), (12. 



15) 



Ay / d^r E cos[go(w< - 



where = J d^r, the prime on the sum now excludes 
a — [3 and/or a' = /?' terms, and 

T:f{v-v')^ 



where 



GL(r = 0)-G> (r=.0) 



in G'^^(O) there is no impHed sum over a and 



J Ac 

Using 



(27r)2 2^ 



^92 



Gai3{q)e 



»;<'(i-5o,3) 

7 

(12.9) 
(12.10) 

(12.11) 



(cos[9o(ua(r) - u/3(r))] cos[go(-Ua'(r') - u/j' (r'))])' 



1 ^ ^iq(u<(r)~u<(v))+iq'(ul,(v')-u<(v')) „ j.') 

(12.16) 



In the above expression the q and g' variables indepen- 
dently range over two values ±go, and S^^f (Sr) (with 
i5r = r — r') is the result of the cumulant average over u> 
fields given by 



Kql + Bq] 



:{l~Sa.0), (12.12) 



which is obviously independent of the tilt-disorder vari- 
an ce A^, and Fourier transfor ming t he last equality in 
Eq.12.12, with 77 as given in Eq.3.19a, we obtain: 



(12.13) 



_ ^^^qiu> (r)-u> {r))+iq'{ul, {v')-u>, (r')) 

-■3-3'(G>^,(<5r)+G>^,(5r)-G>^,(5r)-G>^,(5r)) _ ^ 
-giifcff+fcTJ') e-9<;'G:^(i5r)((5„^/+<53^'-i5„^/-<5^„/) _ ^ 

(12.17) 



Hence, for an infinitesimal S£ the "graphical" correction 
to Ay to first order in Ay is: 



SA 



(1) 



-?7Av 



(12.14) 



We now proceed to compute the corrections to SH[u^] 
to second order in Ay, which we anticipate will result 
in higher order corrections to the random field disorder 
(5 Ay, as well as a correction to tilt disorder variance A^^. 
We compute 



where to obtain the last equation above we used Eq. 12.3a . 
The real space propagator of the short wavelength modes 
((5r) is by definition a partial Fouri er transform of the 
propagator in momentum space, Eqs. 12.3a 12.3b| , which 
for infinitesimal rcnormalization group transformation is 
integrated only over an infinitesimal shell of wavevectors 
Ae~^^ < q± < A. This implies that the exponential in 
the last equation can be expanded in powers of an in- 
finitesimal (Sr) proportional to 5£ 



4qq'G^{dr)Sa,a' + ^{G^{Sr)Y{Saa> + 5pp, - S^^p, - Spa.,f + ... 

4:qq'G^{Sr)Saa' + 2qQ{G^{Sr))^ {Saa' + Saa'Spp' - 2Saf3'Sf3f3') + 



(12.18a) 
(12.18b) 



where to obtain the first terms in the respective square 
brackets of Eqs. 12.18a 12.18b, above, we anticipated that 

a'p' 



a,a',/3,/3'. Inserting this last expression Eq, 12.18b into 



T^J^ ((5r), Eq.12.16, and performing the sum over q,q' 
we obtain 



will be summed over q and 



which allowed 



will uc cuiiiiiivju uvci q ana q 
us to make simplifying changes in summation variables 
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where 



(12.19) 



^±"J'(r, r') = (cos[go«(r) - ^^|(r) - u<{v') + ^(r'))] ± cos[go«(r) - u<{v) + u<{v') - ^(r'))]) (12.20) 



Naively, it appears that the 5af}, Sa'jS' and —^Sap'Sppi 



ious terms with the final results, of cours e, unch anged 



terms above contribute to the renormalization of Ay to Inserting the resulting expression Eq. 12.19 into Eq.12.15 
second order in Ay. However, upon summation over we obtain 
q;, /9, a', /?', these terms, in fact, vanish, because in the 
sum over replica indices we exclude the diagonal terms 
a — (3 and a' — 13' . If we had kept these diagonal con- 
tributions in our definition of H/^, a (somewhat less ob- 
vious) cancellation would have taken place between var- 



h 



{Hl)^=h+h + h, (12.21) 
where we have defined three contributions to {H^^^, 
2G>{5r)+q^,{G>iSr))^] cos[go«(r) - u< {r) ~ u< {v' ) + u< {v' ))] (12.22a) 

~ 2G> (Sr) + ql (G> {5v)) cos[q^{u< (r) - u< (r) + u< (r') - u< (r'))] , (12.22b) 
h^^l [ ^g4^-2,^G>(0)(G>(5r))'x 

(cos[go«(r) -?.|(r) -<(r') +^i|(r'))] +cos[go«(r) -^.<(r) +<(r') -"|(r'))]) , (12.22c) 



r,6r 



E 



ql e-29o^«?(0) 



We now carefully analyze each of the above terms. Al- 
though naively Ii and I2 appear as 3-replica terms, in ~ 2-^ 
fact, as we will see below, Ii renormalizes both Ay and 
A/i, and I2 is irrelevant. 

Ii naturally splits up into a. [3 — (3' term and every- 
thing else: 



cos[go«(r) - u<(r) - <(r + 6y) + u<,{v + 5y))] , 

I 

^ J2 cos[go(u|(r)~u<(r)-K<5r- V«(r)-u<(r)))] , 



where 



ql e-2«oG?(0) 



r.dr 



2G> (Sr) + ql (G> (<5r)) '] (/i, + lu) , (12.23) 



~ E (" ~ ^) cos[(7o(M^(r) — li^, (r))] + irrelevant terms , 

(12.25) 

where the "irrelevant terms" are of the form of a prod- 
uct of cos[(7o(ug (r) — Ug,{r))] and derivatives of (r) — 



/ia= J2 cos[qo«(r) -?/<(r)-u<(r')+u|,(r'))] , 

(12.24a) 



Up, (r) , and are clearly less important at long length scales 
than the term that we have kept. Similarly for /if,, we 
obtain. 



hb = cosbo(w<(r) - u<{r) - u<{r') + u<{r'))] . hb = J2 cos[go(w<(r) - u<{r) - M<(r + 6r) + u|(r + Sr))] 



(12.24b) 



Since the function G^ (Sr) multiplying the above expres- 
sions is a short-ranged function of Sr (with range on the 
order A~^, the inverse of the ultra-violet cutoff A), we 
can safely expand above expressions in powers of Sr. 



J2cos[qoiSr-Viu<ir)^u<{r)))], 



El' 



|Jr.V«(r)-«<(r))p 
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+irrelevant terms , 



(12.26) 



Performing si milar loc al expansion in powers of Sr of 
the term I2, Eq.l2.22t, it is easy to see that I2 con- 
tributes terms of the form of a product cos[qo(2u'^{r) — 
Up{r) — u^,(r))] and derivatives of u^{r) — u^,{r), which 
are all irrelevant at long scales, near the glass transition 
temperature, and we therefore drop them. 

Finally, the analysis of I3 is also very similar. The sec- 
ond cosine in Eq.l2.22c generates a higher harmonic and 



is therefore irrelevant as in the analysis of I2, while the 
first cosine gives 

cos[qo(w<(r) - u|(r) - it<(r + 5r) + w|(r + Sr))] 

J r,5r ^ _LQ 



l_2£|5r.V«(r)-^<(r))|' 



irrelevant terms 



(12.27) 



Putting all this together into Eq.12.21 and dropping ir- 
relevant and constant terms (which renormalize the con- 
stant part of the free energy), we obtain 



2 \= 
> 



-A2.e-29oG>(0) ^ 



^ [A(2-n)cosbo«(r)-u<(r))] 

ic,,9.(u<(r) - u<{r))d,[u<{v) ~ u<(v)) 

(12.28) 

where we have defined 



-I- 



A 



C, 



Sr 



2G>{Sr) + ql{G>{Sr))' 
G>iSr)+q',{G>{Sr))' 



Sri Sri 



(12.29a) 
(12.29b) 



The first term in A and Cij is just a Fourier trans- 
form of G'^(r) evaluated at q = 0. Because the propa- 
gator of the high wavevector modes, G'^(q) by definition 
only has support near the lattice cutoff q± — A, with 
Gy(q = 0) = 0, the first terms in A and vanish. 
With this simplification A and Cij can be evaluated via 
contour integration. 

For A we easily find. 



qt {G>{Sr))' 

JSr 



4 



fq± 
(27r)2 



(12.30a) 
(12.30b) 



4 



(Pq± 



dqz 



1 



(12.30c) 



Ae-« (2^)2 2n {Kql + BqD^ 

Si (12.30d) 

(12.30e) 



4 



AiSi . 



The evaluation of Ca is a similar but a little more in- 



1I 



« 



volved. Because of the anisotropic scaling, q^ 
g^, enforce d b y the quadratic part of the elastic Hamil- 
tonian, Eq.3.8, it is clear that we are only interested in 
Cij with i,j taking on values x, y, with 



C^^2g6 / {G>{Sr))'SrlSrj- 

JSr 



(12.31a) 
(12.31b) 



where the second equation follows by rotational invari- 



ance in the xy-plane, which enforces Cxx = Cy 
and Cxy = Cyx = 0, with 



C^ = ql I d\dzrl{G>{r^,z)y 
-.St , 



cell 



= A^Si . 



C± 

(12.32a) 

(12.32b) 
(12.32c) 



In the above c is a dimensionless constant of order 1. 
Now comparing the expressi on fo r {H'^') , Eq. 12.28, with 
the "bare" Hamiltonian, Eq.3.8, we fincf 



SA,, 



n)AiAlrS£ 



-AiAlSl , 



A^AirSi 



(12.33a) 



(12.33b) 



We note that the graphical correction to the random 
field disorder Ay is stabilizing (negative), leading to a 
stable Cardy-Ostlund-like glassy fixed line (Fig. 6), only 
for n ^ (n < 2), which can be physically understood 
from the discussion given in Sec. III. 

After the above integration of the short wavelength 
modes u^, the effective Hamiltonian is a functional of 
the long wavelength modes , which have an ultra-violet 
cutoff Ae^^ that is different from the cutoff A of the orig- 
inal theory. Therefore in order to identify the new ef- 
fective coupling constants Ay {t) , Ah{t) , ■ ■ ■ we need to 
perform the second part of the renormalization group 
transformation, that involves the rescaling of variables, 
which restores the cutoff in the effective theory back to 
A. The rescaling that accomplishes this is 



u<{v) ^ e^'u^{v') 



(12.34) 
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Because the random-field nonlinearity is a periodic func- 
tion, it is convenient (but not necessary) to take the 
arbitrary field dimension = 0, thereby preserv- 
ing the period p25r/(7o under the renormalization group 
transformation^ Under this transformation the result- 
ing effective Hamilto nian functional can be restored into 
its original form, Eq.3.S with the effective ^-dependent 
couplings, that satisfy differential recursion relations. 



Eqs.3.18 



XIII. APPENDIX B: DETAILS OF THE 3 < L» < 7 
FUNCTIONAL RENORMALIZATION GROUP 
ANALYSIS, INCLUDING ANHARMONIC 
ELASTICITY 

In this appendix we present the details of the func- 
tional renormalization group (FRG) calculation for the 
effecti ve re plicated Hamiltonian, H = Hq + iJint, given 
by Eq. B.IO, with the quadratic part Hq 



1 /" " 



(13.1) 



and the anharmonic part of the effective Hamiltonian 

-ffint = i?anh + -f^A givCU by 



1 /■ " 



B 



(13.2a) 



f " 1 



Ay(Ua - U/3)j 



(13.2b) 



which incorporates both the elastic anharmonicities, pre- 
viously studied in Sec. VI, and the random field and tilt 
disorders of a randomly pinned smectic liquid crystal. 
The result of the analysis presente d below is the se t of th e 
FRG flow equations given in Eqs.|]T^|l3| of Sec.[vni. 

Our RG analysis will closely follow the approach pre- 
sented in the previous RG sections, in particular paral- 
leling the calculation and notation for a single harmonic 
disorder in three dimensions, presented in Appendix A. 
The difference here is that we are interested in the ef- 
fect of disorder in dimensio ns hig her than 3. As a result, 
as discussed in detail in Sec. VIII, in these higher dimen- 
sions all harmonics of the random pinning potential are 
equally relevant and need to be kept track of by studying 
the evolution of arbitrary periodic functions Ay (u^ — up) 
and A/i(uq — up) under the RG transformation. Further- 
more, in contrast to Appendix A, here, in addition, we 
are including anharmonic nonlinearities, which lead to 



the weak logarithmic anomalous elasticity discussed in 
Sec. VIII for 5 < d < 7, and to the power-law anomalous 



elasticity for d < 5 discussed in Sec. VI 



A. Assuming tilt disorder is irrelevant for d > 5 



We have argued in Sees. HI and IV, that, although in 
a full model of a randomly pinned smectic both random 
tilt and field disorders are present, for 3 < d < 5, the 
random tilt disorder always dominates over the random 
field (periodic) disorder. However, as the power-counting 
suggests, for d > 5, the random tilt disorder becomes ir- 
relevant and smectic correlations are dominated by the 
random field disorder. In the second subsection of this 
appendix we will explicitly demonstrate that this con- 
clusion remains valid even when potentially singular di- 
agrammatic corrections to A/i(wq — up) are taken into 
account. Hence, for the remainder of this subsection, con- 
fining our discussion to dimensions d > 5, where random 
tilt disorder is irrelevant, we leave it out of our analysis. 

The logic of the renormalization group analysis here 
is the same as that used in Appendix A: we separate 
the fields Ua into high and low wavevector components 
and , respectively, and perturbatively integrate the 
out of the replicated partition function. At leading 
order, the correction to the quadratic Hamiltonian is 
formally given by a cumulant expansion 



5H[u<] = (Hi„t[u< 



(13.3) 



where as in Appendix A, we have set T — 1, but in con- 
trast to the calculations there have not excluded diagonal 
terms in the replica sum. 

We first focus on the contributions to the random field 
disorder Ay(u). We begin by noting that the renormal- 
ization of Av(u) gets no contributions from the elastic 
anharmonicities, dzUa (Vj_Ma)^ and (V_lUq) , i.e. from 



the vertices in ifanh, Eq. 13.2a. This is because the graphs 
that look like they might renormalize lS.y(ua — up) have 
g's on the external legs and hence renormalize only the 
tilt disorder (and other less relevant operators), which, 
as we discussed above, is irrelevant for d > 5, which is 
the case we are considering here. 

To first order in the disorder Ay{ua — up), keeping 
only the relevant, non-constant terms, we obtain 



^^^A^ = -E / {^v{Uo.-Up))^, 



< _ n.<\/(,.> _ „>\2 



- Un)((U„ - U 



-d£ 



a,P 



J21K{<~u<), (13.4) 



which when compared to the definition of H^^ gives 
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2VBK 



with the primes indicating a partial derivative with re- 
spect to u. To obtain the final result above, we used the 
Gaussian part of the total effective Hamiltonian 



to calculate 



5H 



(2) 



4 ^ 

ai,;3i,a2,/32 



(13.5) ((«> - u>Y)^ = 2 / [(l^.(q)P) - («„(q)uM-q))] , 



(13.7) 



using the corresponding propagator 



(Ma(q)u/3(-q)) = 



Bql + Kqi [Bql + ifgjl^ 



(13.8) 



(13.6) obtained from -ffcauss by manipulations virtually iden- 
tical to those described earlier for the tilt-only model, 
Seclvi. 



The contribution to second order in Av(u) is given by 
A'^«(ri) - (ri))A^«(r2) - u<^{v^))l2fpl{v,,v^) , 

(13.9) 



where 



5H 



(2) 



«(ri)-<(ri))'«(r2)-^>(r2))^ 



(13.10a) 



:<(rl)-<(rl))«(r2)-^.>(r2)) 



(13.10b) 



G^(ri -rs) 



(13.10c) 



5t = Ti — r2, we used the Wick decomposition theorem 
for Gaussian random variables (a, b, c, d) 

{abed) = {ab){cd) + {ac){bd) + (ad) (be) , (13.11) 



to go from Eq.13.10a to Eq.l3.10b|, and used the symme- 

■ (21 ^fl T] 

try of the summ and in S H\ , Eq. |13.9| under a, (3 inter- 
change to get Eq. 13.10d. 

Substituting the last expression for I^l^^i Eq. 13.10d 



)(G>(5r))' 



into SH^^' above and using the short-range property of 
G^{6r) to expand ri as ri = r2 -t- 6r, keeping only the 
(13.10d) most relevant terms, we obtain 

-SiG2 IJ2 [K i< (r) - u< {r)f - 2 A'^ « (r) - u< (r)) A^ (0) + J] K i< (r) - u< (r)) A^ « (r) - u< (r))" 

■ Q,/3 7 

(13.12) 



where the constant G2 is defined by 
G2S£= f {G>{Sr)Y , 

dqz 



i5r 

A ^-l, 

Ae-« (2t 
A 



Ae-« (27r 



\d-l 
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1 



27r (iCg^ + Bq 



the two-replica terms and can therefore be neglected. 
Co mparin g the resulting expression for 5H)^ ' with Ha , 
(13.13a) Eq.^3^wefind 



(13.13b) SAPiu) « 6e ^^'-^' (AUuf - 2A'^(^.)A'^(0) 



sVk^b 



(13.13c) 



(13.14) 



Combining the first and second order contributions to 
(13.13d) Av(u), Eqs.p^, [13.14 with the length and field rescal- 



ings, Eq. 12.34, necessary to bring the ultra-violet cutoff 
It is easy to see by power- count ing that the three- rephca back to A we obtains the recursion relation for Ay('u) 



(2,) I ri 

(last) term in SH}^\ Eq. 13.12 is irrelevant relative to given by 



d,Av{u,e) = {d+l)A{u,£)+'^A'{,{u,e) 



-{AUu,e)f-A'{.{uJ)A'{,{Q,e) 



Cd-iA 



d-7 



(13.15) 
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where, for simplicity, in the above we have set the field 
rescaling (f> — and used to = 2. 

We now turn to the renormalization of the elastic mod- 
uli B and K. As discussed in Sec. [11, statistical sym- 



metry forbids renormalization of these by disorder alone 
(The same argument applies to the XY spin stiffness in 
the random field XY model). However, anharmonic elas- 
ticity, special to smectic liquid crystals, conspires with 
the random field Ay(u) to renormalize B and K. 

Calculations very similar to those for the tilt-only 
model treated in SecJVl| show that the graphical correc- 
tions to B and K found in that section can be taken over 
for the random field disorder, with the identification 



A'f = -2A'^(0)/A2 



(13.16) 



In addition to a detailed perturbative RG calculation, 
this can also be easily seen by comparing 

Av{u^ - up) « const. + iA'^(0)(u„ - upf , (13.17) 

to the disorder part of the the tilt-only Hamiltonian. 
Using this mapping we immediately find 

1/2 



dK 



12(i + 23 ( B 



dB 
d£ 



8(rf2 - 1) 
d-3 + l(^)KiO)Cd^iA'^' 



K 
(13.18a) 



B 



(13.18b) 

Now identifying the dimensionless measure of the random 
field disorder 



Ay (m, i) = ^d-lA'^ ^ ^ 



(13.19) 



and using it inside Eqs. |l3.15| |l3.18a| a nd |l3.18bi we ob- 
tain the FRG flow Eqs j8.13| , Is.lll and ^.12| given in the 

Seclvn^. 



B. Proving that tilt disorder is irrelevant for d > 5 

In the previous subsection we have ignored the tilt dis- 
order for d > 5, arguing for its irrelevance based on 
power-counting in these higher dimensions. However, 
the (supernaturally) alert reader might be concerned (as 
we were at first) that the singular behavior of Ay(u), 
Eq. |8.23 (in particular, the divergence of its fourth deriva- 
tive with respect to u, Ky"\u,l) as £ — > oo) could po- 
tentially invalidate the simple power-counting argument 
that tilt disorder is irrelevant for d > 5. Here, we demon- 
strate that this does not happen, by deriving a. Junctional 
renormalization group recursion relation for Ah{u,l) to 
linear order in Ah{u,t) itself at the (FRG) fixed point 



A*y{u), Eq.3.23, found in SccVHI, We obtain 



d,Ah{u,i) = (d- 1) Ah{u,l) + + D{u)^Al{u,i) + c ^K%i)Bie)A'-'/Cd- 



(13.20) 



r 



where we have defined the "position-dependent diffusiv- 
ity" 



combination 

g{u,i) EE A^(u,^) (;§|y) ^ C.^iA'-' , (13.22) 

and c and C2 arc dimensionless constants of order unity. whose recursion relation readi ly follows from th ose for 
As in the tilt-only model, the important quantity K{£), B{£), and Afi{u,£), Eqs. 13. 18a , 13.18b , and 13.20 
proves not to be Ah{u,£), but rather the dimensionless respectively: 



D{u) ^ C2 [a'Uu) - A'{,m) , 



(13.21) 



degiu,£) 



A2(£) 



XH£) 



(13.23) 



where cg = 3/4+ (d^ - 12d + 23)/(2(d2 - 1)). 

The solution to this inhomogeneous linear partial dif- 
ferential equation can, as usual, be written as a sum of 
the solution to the homogeneous equation obtained by 
dropping the last term, plus any solution of the inhomo- 
geneous equation. 

The solution to the homogeneous equation can be writ- 
ten as an expansion in a complete basis 



ghom{u,£) 



oo 

E 

n=0 



gni£)(f>niu) , 



(13.24) 



where the set </>«(«) are the eigenfunctions of the operator 
D{u)d'^ defined and discussed in Sec. VIII. 

The e xpansi on coefficients gn{£) can be readily found 
from Eqjl3.23| to be 
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exp 



A2(F) 



(13.25) 



with 



A„ = rf-5 + |r„| 



(13.26) 



which is positive for all n for d > 5, and we have 
dropped the r](£) term, which vanishes exponentiall y fast 



13.25 



as i —> oo. Since the second exponential in Eq 
is a monotonically decreasing function of £, eve ry gn{ £), 
and, hence, the entire homogeneous solution Eq. 13.2^ for 
g{u,£), vanishes as ^ oo. 

Thus, the only possible problem with neglecting tilt 
disorder (for d > 5) must come from the inhomogeneous 
solution for g(u,£). We will now show that this vanishes 
as well, thereby proving that tilt disorder is irrelevant for 
d > 5. 

We do this by finding an asymptotic (large £) solution 
gi{u,£) to the full, inhomogeneous equation of the form 



gi{u,£)^Y. 

n=l 



fn{u) 



(13.27) 



which clearly vanishes as £ — > oo. 



Inserting Eq. 13.27 into Eq. 13.22, and keeping in mind 



that l/X'^{£) is itself cx 1/^, as £ ^ oo (see Eq.^.l7D, we 



find an equation for /i {u) by equating coeffic ients o f the 
leading order l/£ terms on both sides of Eq. 13.23, The 
resulting equation reads 



Diu)f'^iu)-{d-5)h{u) 



cA(d)|A- 



(0)1 



(13.28) 



where we ha ve de fined A{d) = (d^ - l)/{d'^ + 6d - 13) 
and used Eq j8.17| for A(^). Once fi{u) is known, the re- 
maining /„'s for n > I can be determined by equating 
coefficients of 1/^ " on bo th sides of Eq. 13.23 after insert- 
ing the ansatz Eq. 13.27. This gives 



I)(u)/:(u) -id- b)Uu) - (c3A(d) + 1 - n)/„-i(u) . 

(13.29) 

Now, to complete the proof that gi{u,£) vanishes as 
oo, w e need merely show that all /„'s obtained from 



Eq. |l 3.281 and Eq. |13.29| are finite for all n. If /i(u) is 
finite, then clearly f2{u) is also finite, since no deriva- 
tives of /i(u) (which c ould p otentially diverge if fi{u) 
is cusped) appear in Eq.13.29. This finiteness argument 
then obviously recursively carries through via Eq jl3.29| 
to all of other /„(m)'s. 

Although fi{u) can be found explicitly applying stan- 
dard integration methods to Eq. 13.28| , our only goal here 
is to demonstrate that /i(u) is finite. To do this we sim- 
ply expand /i (u) in eigenfunctions (j)n (u) of the operator 
D{u)dl: 



(13.30) 



Inserting this expansion into Eq. 13.28, multiplying both 
sides by 4im{u), integrating from u — to w = a, and 
using the eigenval ue eq uation, Eq. 8.46 and the orthogo- 
nality relation Eq. 3.47, we obtain 



zA{d) 



\KM\ 



(13.31) 



where the source 



IynD{u) 



du 



(13.32) 



was show n to be finite for all m (including m = 0) in 

Seclvm. 



Hence, fi{u) will be finite for all u, provided that the 
sum 



E°° Sn4> n{u) 



I A", 



d-5 + |r„| 



cA{d) 



-AH 



(13.33) 



converges for all u. 

To prove that it does, let us follow the obvious con- 
vention of ordering the eigenmodes n such that |r„| is a 
monotonically increasing function of n. Then, for large 
n, we can use the WKB solution for the eigenfunctions 

(f>n{u) 



and eigenvalues r„ 

r„ = 



■1/2 



nqoD 



-Dn\l 



du' 



D^^iu') 



(13.34) 



(13.35) 



where we have defined the suitably averaged diffusion 
constant 



D = 



du' 



a Jo D^/^{i 



(13.36) 



to prove the convergence of the sum in Eq. |l3.33 
Eq. 13.34 shows that (pniu) is bounded from above for 
all u (by [maxu£'(u)]^/* ), which in turn implies that the 
Sn's are also bounded above (since A'y^{u) is fini te for 
all u, and so is (/)«(«)) as n ^ oo. Then using Eq. 13.35 
for the r„, we see that the large n behavior of the sum- 
mand in Eq. 13. 33] is bounded above by const. /n^. Hence, 
the sum converges, and so fi{u) is in fact finite for all 
u; consequently all /„(u)'s are also finite. This therefore 
demonstrates that g(u, £) vanishes as ^ — > oo, proving our 
assertion that tilt disorder is irrelevant for d > 5. 
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XIV. APPENDIX C: DERIVATION OF SMECTIC 
"COULOMB GAS" DISLOCATION THEORY, 
WITHOUT EULER-LAGRANGE EQUATION 

In this appendix we derive the Coulomb gas description 
of smectic dislocatio n loo ps characterized by an effec- 
tive Hamiltonian, Eg. 5. 15 . However, here, in contrast to 
the derivation of the main text, we accomplish this with- 
out minimizing the Hamiltonian with respect to smooth 
smectic deformations. 

Our starting point is the elastic Hamiltonian for a tilt- 



only model of a randomly pinned smectic, given in Eq.5.1 
'B 



H[u] 



K 



(Viii)" + h(r) • V^u 



(14.1) 



We include dislocations in above description by allowing 
the layer displacement u(r) to be a multi-valued function, 
such that. 



V X Vu = m 



(14.2) 



We then decompose Vu into a sum of a singular, purely 
transverse part Vd satisfying 



V X Vd = m 



(14.3) 



and a nonsingular, purely longitudinal part Vup, accord- 
ing to 



(14.4) 



The difference between this approach and our earlier 
derivation in Sec.^ is in our arbitrary choice of the sepa- 
ration of Vm into phonon and dislocation parts. In Sec.^, 
we chose to minimize H given m; here, we simply 
choose it to be purely transverse. We will now demon- 
strate that this difference between these two arbitrary 
choices has no effect on the final answer, as it should 
not. 

Substituting this decomposition into H[u\, Eg, 14.1 , we 
obtain 

Htot[Up,^fd] Hc\[u.p\ +-ffv[vd] +i?mt[Wp,Vd] , (14.5) 

where (dropping the subscript d on v), 



and 

Hint 



B 



K 



(Vx-vx)2 + h(r).vi 



d^Y 



Bv^d^u + K{V± ■ v_L)Viw 



BvzdzU — KdzVz^\u 



(14.6a) 
(14.6b) 

(14.7a) 
(14.7b) 



In go ing fr om Eq. 14.6a to Eq. l4.6b| and from Eq. 14.7a 
to Eq. 14.7b, we used the purely transverse property of v, 
V • V = 0, or equivalently. 



V_L • V_L + dzVz = 



(14.i 



to eliminate v_l in favor of Vz ■ 

We now integrate over the single-valued phonon de- 
gree of freedom Up. After a simple Gaussian integration, 
some algebra, and Fourier transformation, we obtain an 
effective Hamiltonian 



H, = 



2f; 



|«,(q)P+b(q).h(-q) 



where, Tq — + ^^"^ t)(q) is given by 



^^(q) = 



^ 9 



w_Li(q) 



with, 



±ij 



(q) = St 



1l 
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(14.9) 
(14.10) 

(14.11a) 
(14.11b) 



the transverse and longitudinal projection operators in 
the space (_L) perpendicular to z. 



Now Fourier transforming Eq.14.3, and solving it for 
v(q), keeping in mind that v is purely transverse, we find 



Wi(,qj 



ieijkqjrrik/q^ 



(14.12) 



Using this solution to eliminate v inside H^, Eq.14.9 



above, in favor of the dislocation density m, we obtain 
the final expression for the dislocation loop Hamiltonian, 
Eq.5.15, quoted in the main text. 



XV. APPENDIX D: ANALYSIS OF 
FLUCTUATIONS IN THE DUAL MODEL OF 
DISORDERED SMECTICS IN TYPE I LIMIT 

In this appendix we analyze the type I regime of the 
dual disordered smectic liquid crystal model derived and 
studied in Sec.^. In this regime, we will compute ex- 
actly the effective free energy and as a by-product obtain 
from it the fluctuation corrections to the reduced dual 



transition temperature, given in Eq. p.37 . 

Our starti ng po int is the replicated dual "action" Sr 
given in Eq.5.33, together with the quenched "gauge 
field" a variance, given by Eq. 5.34b. We build on the 

liF 



a variance, 

ideas developed in Refs. |15|, generalizing them here to 
disordered systems. As discussed in these references and 



in Sec.^ in the type I regime the order parameter ip is 
significantly stiffer than the gauge field A and can there- 
fore be accurately treated in a mean-field approximation. 
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This amounts to taking ■0(r) to be constant in space. The 1 , -pr |.^^ , 

constant V' is then calculated by minimizing the resultant ~ Na J 'i-L^ ^ 

free energy, which can now be computed exactly, without 

further approximations. As is clear from the discussion 

in Sec.0, to do so, we must calculate with 



•S„y:[V,AQ,a] 



Smf[lp,Aa,a] = Sr[lp,Aa,a] + ^ ^ ^ % „2 -^j (q)^^ (~q) ' 



^ a a 



(15.1) 



(15.2) 



(15.3) 



The prime on the integral in Eq.15.1 indicates that it 
is constrained to be taken only over the transverse parts 
of Aa and a, V is the volume of the system, Na is the 
normalization factor coming from the probability distri- 
bution of a, and will be chosen such that Z" ^ 1 as 
n ^ 0, and 



(15.4a) 



In Z« 



n(t|^P+.|V.n + ^^^^ / ln(G-(q) 



(15.4b) 



can be read off from Eqs. |5.33| and 5.341 , taking the ap- 
propriate long wavelength limit in the latter. 

The great virtue of the effective "action" Smf is that it 
is quadratic in Aq and a and therefore allows exact eval- 



uation of the functional integrals over them in Eq.15.1 
Taking these Gaussian integrals, we obtain 



(d-2) 



ln(G-i(q)-n|Vp) -IniV, , (15.5) 



where the factors of d — 2 correspond to the fact that 
there are d — 2 transverse "gauge-field" components of 
Aa and a, due to the transversality constraint on the 
functional integral discussed above. 

We observe that, as expected, the normalization factor 
Na for the probability distribution of a, given by 



lniV„ 



guarantees that In Z" ^ as n 0, 



In 



-Vn 



(d-2) 



ln(G:,Hq) + IV^I') + 



Incorporating Na into Eq.15.5 gives 
{d~2 



2n 



In l-nGa(q)|^|^ 



(15.6) 



(15.7) 



Expanding this in n, we find: 



InZ" = -Vn 



d-2 



Ga(q) +0{n^) 



(15.8) 



Now using the standard replica expression for the av- by 
erage free energy, we finally obtain 



lim — 

Fmf = ?i ^ -ksT 

n 



1 



lim 

n ^ 



-kuT 



keTV 



n 

|2 I „,|„;,|4 



tnm'+u\ip\' 



d-2 



(15.9a) 



tR^t 



d-2 



G-^(q) , 



(15.10) 



(15 9b) ^hi*^^! using Eq. [l5.4b| , is precisely the result used in 
Eq.5.37 of Sec^ to assess the stabihty of the disordered 



(15.9c) 



smectic to dislocation unbinding. Thus, the result of our 
more rigorous and systematic perturbation theory is re- 
covered by this fluctuation corrected mean field theory. 
The annealed term in Eq. 15.9c is analogous to the term 



where the renormalized reduced temperature ta is given that leads to a fluctuation-driven first-order transition in 
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type I superconductors. Here, however, its effects are 
innocuous: expanding the integrand in ji/jp leads only to 
a finite renormalization of t, plus a non-analytic |'(/;|'^"'"^ 
piece, which only leads to a finite renormalization of u in 
d = 3. 
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